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a  100(usec^  CW  pulse  centered  at  80  kHz  at  the  wind-roughened  surface  of  a 
freshwater  lake  at  a  10.5  degree  grazing  angle.  The  returns  were  sampled 
digitally  for  a  10  msec  interval  during  which  surface  reverberation  was 
predominant,  and  were  formed  into  ensembles  at  each  sample  timej_>  Various 
univariate  moments  and  the  tipe-difference  component  of  the  covariance  were 
computed  experimentally ^utilizing  ensemble  averages.  The  statistical  validity 
of  the  ensembles  was  determined  by  testing  each  ensemble  for  randomness 
and  homogeneity.  The  ensembles  were  also  tested  for  normality  and  found  to 
be  non-Gaussian,  in  contrast  to  previous  studies. 

The  experimental  results  were  compared  to  a  theoretical  model  developed 
by  D.  Middleton  and  others,  based  on  point  scattering  and  Poisson  statistics. 
The  reverberation  is  modeled  as  weak  scattering  from  random  point  sources 
representing  inhomogeneities  at  the  air-water  interface  of  an  otherwise 
homogeneous  medium.  In  the  present  study  the  point  source  was  assumed  to  be 
a  perfect  point  reflector  distributed  uniformly  over  the  surface.  The  imple¬ 
mentation  of  the  model  took  into  account  most  of  the  relevant  geometrical 
parameters  (spatial  distribution  of  sensors,  grazing  angle,  range)  and  acoustic 
parameters  (frequency,  pulse  length,  directionality,  aperture  response,  band¬ 
width,  signal  spectrum)  but  did  rot  take  into  account  environmental  parameters 
(surface  wave  height,  wave  spectrum,  wind  direction). 

'  Significant  differences  were  observed  between  the  covariances  of  the 
vertical  and  horizontal  arrays.  The  envelope  of  the  covariance  between  the 
vertical  receivers  maintained  a  significant  level  at  much  larger  separations 
than  the  horizontal  receivers.  The  phase  of  the  vertical  covariance  was  found 
to  change  linearly  with  time,  resulting  in  a  slow  oscillation  of  the  covariance 
with  time,  while  the  phase  of  the  horizontal  covariance  was  constant  (and 
non-zero).  The  vertical  covariance  was  also  shown  to  depend  on  the  location  of 
the  elements  on  the  array. 

The  theory,  which  had  been  applied  previously  only  to  horizontal  arrays, 
was  extended  to  a  vertical  array.  The  theory  correctly  predicted  the  dependence 
of  the  horizontal  and  vertical  covariance  on  the  time  of  observation  and  on  the 
time-difference  of  the  observations.  In  addition,  it  predicted  the  dependence 
of  the  envelope  of  the  covariance  on  horizontal  separation  of  the  receivers. 

The  change  of  the  phase  of  the  vertical  covariance  with  time  was  also  predicted 
by  theory.  However,  the  theory  failed  to  predict  the  dependence  of  the  envelope 
of  the  covariance  on  vertical  separation  and  the  non-zero  phase  of  the  hori¬ 
zontal  covariance.  The  two  failures  of  the  theory  were  shown  to  most  likely 
be  the  result  of  the  failure  to  include  the  environmental  parameters  into  the 
implementation  of  the  model.  Ways  to  include  these  environmental  parameters 
were  suggested. 
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INTRODUCTION 


The  scattering  of  waves,  either  electromagnetic  or  acoustic, 
from  a  boundary  or  a  volume  has  been  the  subject  of  study  in  many  varied 
fields.  The  eommuni cat  ions  field  studies  scattering  of  electromagnetic 
wav>.s  from  the  ionosphere  to  provide  over-the-horizon  communications. 
Acoustics  has  recently  provided  improved  diagnostic  capabilities  to  the 
medical  community  by  considering  the  scattering  of  ultrasound  from 
organs  and  tissues  of  the  body.  Sound  .scattered  from  volume  scatterers 
can  provide  information  on  the  lumber  of  pollution  particles  in  the 
atmosphere  or  the  number  of  fish  in  a  school.  Other  diverse  fields  such 
as  seismic  exploration,  non-destructive  testing,  and  architectural 
engineering  all  claim  interest  in  the  scattering  of  waves. 

The  present  study  is  most  directly  related  to  the  field  of 
sonar.  Techniques  have  been  developed  in  this  field  to  enable  active 
sonars  to  better  detect  signals  in  the  presence  of  noise  due  to  surface 
tackscattering  (reverberation).  However  many  of  these  techniques  require 
a  knowledge  of  the  statistical  properties  of  surface  reverberation,  in 
particular  the  coherence  of  the  reverberation.  This  requirement  has 
led  to  interest  in  the  fundamental  problem  of  the  scattering  of  sound 
from  the  surface.  The  present  study  is  one  response  to  this  interest. 

The  focus  of  this  study  is  on  the  coherence  of  the  scattered  sound. 

I . 1  Previous  Experimental  Studies 

Much  of  the  previous  experimental  work  which  has  been  performed 
in  this  area  has  investigated  the  coherence  of  forward-scattered  sound. 


1 


o 


1-3  1|_7 

The  sound  source  was  usually  either  continuous  or  explosive;  in 

Q  1  i 

some  cases  a  pulsed  source  was  used.  In  almost  every  case  tie 

received  signal  was  assumed  to  be  stationary  and  ergodic.  Thus  the 

coherence  was  calculated  using  time  averages. 

Other  experimental  work  has  dealt  with  the  coherence  of  back- 

scattered  sound,  or  reverberation.  As  opposed  to  the  experiments  with 

12-10  20-22 

forward-scattered  sound,  the  source  was  pulsed  or  explosive 

rather  than  continuous.  In  addition,  the  coherence  was  almost  always 

computed  using  ensemble  averages  rather  than  time  averages.  Thus  the 

results  were  not  constrained  by  the  assumptions  of  stationarity  and 

ergodicity.  The  present  study  is  limited  to  the  coherence  of  back- 

scattered  sound  from  a  pulsed  sound  source.  Ensemble  averages  are 

used  to  experimentally  estimate  the  reverberation  coherence. 

Most  of  the  previous  experimental  studies  of  the  coherence  of 

reverberation  from  a  pulsed  sound  source  utilize  a  single  receiver  and 

12  1 6  id 

so  only  consider  temporal  coherence.  ’  Two  studies  used  a  multi¬ 

sensor  array  and  examined  the  dependence  of  the  coherence  on  the  spatial 

separation  of  the  observation  points,  as  well  as  the  time  of  observa- 
17  19 

tion.  ’  However,  these  studies  considered  the  spatial  dependence  of 
the  coherence  only  for  a  linear  horizontal  array.  The  present  study 
employs  both  a  horizontal  and  a  vertical  array  to  examine  the  spatial 
dependence  of  the  coherence  for  both  orientations  simultaneously. 

The  spatial  coherence  of  forward-scattered  sound  has  been  exam¬ 
ined  for  both  horizontal  and  vertical  arrays.  However,  since  forward- 
scattered  sound  is  predominantly  specular,  significant  differences  in 


the  spatial  coherence  of  forward-scattered  sound  and  back-scattered 
sound  are  to  be  expected.  In  fact,  a  comparison  of  forward-scattered 

t,  22 

data  by  Wille'  and  back-scattered  data  by  Urick,  both  using  explosive 
sources,  indicates  that  the  back-scattered  coherence  is  significantly 
less.  Thus  it  is  to  be  expected  that  the  present  study  of  back-scattered 
sound  from  a  pulsed  source  will  demonstrate  results  which  differ  signi¬ 
ficantly  from  similar  forward-scattered  data. 

The  spatial  coherence  of  back-scattered  sound  has  also 
been  examined  for  both  horizontal  and  vertical  arrays  utilizing  an 
explosive  sound  source  and  performing  time  averages  to  estimate  the 

op  oo 

coherence."  v  '  The  results  indicated  significant  differences  between 

the  coherence  of  horizontal  and  vertical  arrays.  Thus  it  is  to  be 

expected  that  the  present  study  will  also  demonstrate  significant 

differences  between  the  coherence  of  horizontal  and  vertical  arrays. 

j _ 2  Theoretical  Studies  in  Surface  Scattering 

Theoretical  models  of  underwater  acoustic  scattering  at  an 

air-water  interface  can  be  divided  generally  into  two  .ategories,  de- 

23  2h 

pending  upon  whether  the  boundary  is  treated  as  periodic  ’  or 
25-30 

random.  The  present  study  considers  the  surface  to  be  random, 

that  is,  the  height  of  the  surface  waves  at  any  point  on  the  surface 
is  treated  as  a  random  variable. 

There  are  two  basic  approaches  to  scattering  from  a  randomly 


rough  surface.  One  approach  introduces  a  scalar  wave  equation  and 
expresses  the  scattered  field  in  terms  of  the  Helmholtz  integral  over 
elementary  sources  at  the  surface.  The  boundary  conditions  are  usually 


applied  by  assuming  a  pressure  release  surface  and  utilizing  the 

25 

Kirchhoff  approximation.  The  other  approach,  developed  most  generally 

by  Middleton,  treats  the  surface  as  a  random  distribution  of  point 

scetterers  representing  inhomogeneities  in  an  otherwise  homogeneous 
27-30 

medium.  The  interaction  of  the  incident  acoustic  wave  with  the 

point  scatterers  is  expressed  in  terms  of  an  impulse  response  function, 

eliminating  the  need  for  solving  the  scalar  wave  equation  for  complex 

boundary  conditions.  Middleton's  model  has  been  described  as  the  most 

complete  theoretical  model,  but  one  which  lacks  experimental  verifica- 
31 

tion.  It  is  the  model  which  will  be  used  in  the  present  study. 

1.3  The  Present  Study 

The  present  study  has  four  primary  objectives.  The  first  is 
to  simultaneously  measure  the  dependence  of  the  reverberation  co- 
variance  on  the  spatial  separation  of  the  receivers  in  both  the  hori¬ 
zontal  and  vertical  orientations.  In  the  process,  the  dependence  of 
the  covariance  on  other  parameters  will  also  be  measured.  A  thorough 
analysis  of  the  horizontal  and  vertical  covariance  will  be  performed. 

The  second  objective  is  to  develop  a  theoretical  model  of  surface  rever¬ 
beration  and  compare  it  to  the  experimental  results.  The  model  being 
used  has  previously  been  applied  to  the  covariance  of  a  horizontal  array 
with  good  success.  Thus  the  present  study  seeks  to  extend  the  model  to 
the  covariance  of  a  vertical  array.  The  third  objective  of  this  study 
is  to  experimentally  examine  the  first  four  univariate  moments  of  the 
reverberation,  and  apply  univariate  tests  for  normality.  This  will 
supply  some  information  concerning  the  distribution  of  the  reverberation 
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process.  The  fourth  objective  involves  refinement  of  the  techniques  for 
measurement  and  validation  of  reverberation  from  multiple  receivers. 

Care  must  be  taken  in  the  measurement  process  in  order  to  form  valid 
sample  ensembles  of  the  reverberation.  In  addition,  it  is  important  to 
apply  statistical  hypothesis  testing  techniques  to  validate  the  sample 
ensembles  for  randomness  and  homogeneity  before  utilizing  these  ensembles 
to  compute  the  various  moments  and  covariance. 

The  present  study  is  organized  as  follows.  The  experimental 
equipment  and  measurement  techniques  are  described  in  Chapter  II.  The 
reverberation  data  were  collected  by  scattering  100  ps  CW  pulses  at  a 
center  frequency  of  80  kHz  from  the  wind-roughened  surface  of  a  fresh¬ 
water  lake  at  a  grazing  angle  of  approximately  10.5  degrees.  Approxi¬ 
mately  1000  reverberation  returns  were  recorded  in  a  short  period  of 
time  from  nine  vertical  receivers  and  four  horizontal  receivers. 

The  theoretical  model  used  in  the  present  study  is  developed  in 
Chapter  III.  A  summary  of  Middleton's  theory  is  given  as  it  relates  to 
the  present  study.  Certain  assumptions  are  then  introduced  which  allow 
the  model  to  be  simplified  to  the  form  used  by  the  present  study.  The 
essential  assumptions  are: 

1.  The  scatterer  is  modeled  as  a  perfect  point  reflector, 
i.e.,  scattering  is  omnidirectional,  frequency  insensitive 
(within  the  bandwidth  of  the  signal),  location  insensitive, 
and  introduces  no  time  delay. 

2.  The  signal  is  narrowband. 


These  and  other  assumptions  are  then  applied  to  derive  the  form  of  the 
model  used  in  the  present  study. 


Chapter  IV  discusses  the  statistical  testing  of  the  experimental 
data.  Assumptions  inherent  in  the  sampling  process  which  are  necessary 
to  transform  the  reverberation  returns  into  ensembles  are  discussed. 

Then  various  statistical  tests  for  randomness,  homogeneity,  and  normal¬ 
ity  are  described  and  the  results  of  the  tests  as  applied  to  each  data 
ensemble  are  given. 

The  experimental  estimates  of  the  covariance  are  given  and 
compared  to  the  results  from  the  theoretical  model  in  Chapter  V.  The 
procedure  for  computing  the  moments  and  covariance  of  the  experimental 
data  is  described  first.  Then  the  moments  are  presented.  The  cc variance 
is  then  analyzed  as  a  function  of  time,  time  delay,  and  spatial  separa¬ 
tion  for  the  various  vertical  and  horizontal  receivers. 

Chapter  VI  then  summarizes  the  study  and  draws  some  conclusions. 


i  T .  EXPERIMENT 


Chapter  JI  is  devoted  to  a  discussion  of  the  experimental  equip¬ 
ment  used  to  generate  and  receive  the  reverberation  data  and  to  a 
description  of  the  recording  and  digitization  techniques  necessary  for 
this  type  of  study.  First,  a  simple  analysis  of  the  sources  of  measure¬ 
ment  error  is  performed  to  help  establish  certain  requirements  for 
the  measurement  system.  Then  the  hardware  configuration  and  measurement 
techniques  are  described  which  allow  these  measurement  errors  to  be 
reduced  to  an  acceptable  level.  A  description  of  the  transducers  and 
the  deployment  of  the  experiment  at  the  ARL:UT  Lake  Travis  Test  Station 
is  then  given.  Finally,  the  conditions  under  which  the  data  were 
collected  are  described. 

II. 1  Sources  of  Measurement  Error 

Several  considerations  are  necessary  in  order  to  accurately 
measure  reverberation  coherence  between  spatially  separated  hydro¬ 
phones.  A  reasonable  reverberation  to  background  ratio  must  be  main¬ 
tained  where  the  background  noise  is  due  to  both  the  ambient  acoustic 
noise  and  the  system  electronic  noise.  Any  feedover  of  the  signal  from 
one  nydrophone  to  the  other  signals  will  also  limit  the  accuracy  of  the 
coherence  measurement.  Thirdly,  since  the  reverberation  coherence  is 
calculated  experimentally  from  an  ensemble  of  signals,  the  sample  ensem¬ 
ble  size  limits  the  accuracy  of  the  coherence  calculations.  Finally, 
it  is  necessary  to  maintain  the  phase  integrity  of  the  narrowband 
acoustic  signals  at  the  face  of  the  hydrophones.  Actually,  it  is 
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necessary  only  to  maintain  the  phase  difference  of  the  acoustic  sign-t..; 
between  the  hydrophone  channels  at  the  face  of  the  hydrophones;  that 
is,  the  addition  of  an  arbitrary  constant  phase  shift  to  all  signals 
does  not  affect  the  coherence  of  the  signals.  The  following  discussion 
will  analyze  these  various  sources  of  inaccuracies  and  describe  the 
measurement  equipment  and  techniques  used  to  reduce  the  inaccuracies 
to  an  acceptable  level.  Although  the  analysis  makes  use  of  several 
assumptions  which,  strictly  speaking,  may  not  be  valid  for  the  data 
used  in  this  study,  nevertheless  it  provides  a  useful  insight  into  the 
relative  importance  of  the  various  sources  of  measurement  error. 

To  get  an  estimate  of  the  effect  of  reverberation  to  noise 
ratio  on  the  measurement  of  coherence  between  two  hydrophones,  let  the 
two  reverberation  processes  be  represented  by  and  X^,  which  are 

identically  distributed  zero  mean  gaussian  random  processes  with  a 
2 

variance  of  o  .  Likewise,  let  Xnl  and  X^0  represent  two  additive 

white  noise  processes  which  are  identically  distributed  zero  mean 

2 

gaussian  random  processes  with  a  variance  of  .  The  reverberai ion 

plus  noise  processes  are  then  X,  =  X,  +  X  ,  and  X^  =  X0  +  X  , . 

1+n'  1  nl  2+n  2  n2 

Since  it  is  reasonable  to  assume  that 


s(xAi>  *  e(x2x„2»  ■ 0 

then  Xg+n  are  identically  distributed  zero  mean  gaussian  random 

2  2 

processes  with  variance  of  a  +  o^  .  Since  it  is  also  reasonable  to 


assume  that 


E(XaXnl! 


*  E(X1X„2: 


"  E(KnlX„2)  '  0 


then 
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The  eorrejation  coefficient  between  X, ,  and  is  then 

l+n  J+n 


lOv  (X*  X-  ) 

Py  Y  - _ UlLZ  illS 

Xl+nA2+n  o 
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The  rms  reverberation  to  rms  noise  ratio  is  R/U 


(R)2 

'  N ' 


Xi+nX2+n 


♦  <§,*  ‘  V, 


— .  Thus , 
a 
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This  gives  an  expression  for  the  measured  reverberation  plus  noise 
correlation  coefficient  in  terms  of  the  actual  reverberation-only 
correlation  coefficient  and  the  reverberation  to  noise  ratio.  The 
measured  correlation  coefficient  is  always  less  than  or  equal  to  the 
actual  correlation  coefficient,  and  the  difference  decreases  as  the 
correlation  coefficient,  decreases.  For  example,  if  R,'N  =  10,  the 
measured  correlation  coefficient  is  .99  when  the  actual  correlation 
coefficient  is  1.0.  Thus,  for  a  reverberation  to  noise  ratio  of  20  dB , 
the  maximum  difference  in  the  correlation  coefficients  occurs  for  maximum 
correlation  and  reduces  the  measured  correlation  coefficient  by  only  .01. 

Now  consider  the  effects  of  feedover  between  channels  on  the 
accuracy  of  the  correlation  coefficient.  Let  X^  and  X0  be  as  above. 

Let  Y^  =  aX^,  and  =  aX^ ,  where  a  is  a  constant.  Y^  and  Y^  are  iden¬ 
tically  distributed  zero  mean  gaussian  random  processes  with  variance  of 

p 

( cio )  .  The  measured  reverberation  signal  is  then  the  sum  of  the  actual 


reverberation  signal  and  the  feedover  from  the  other  channel,  i . 
Z1  =  X1  + 

Z2  =  X2  +  Y1 

where  and  are  the  measured  reverberation  processes. 


E(Z)  =  E( X. )  +  E(Yj  =  0 
1  ±  £ 

E(zx2)  =  E{(x1  +  y.)2}  a e  (X,-')  +  v,)  +  e(yJ' 


=  o  +  2a  E(X1  X2)  +  (a  o)‘ 


(1  +  2ap  +  aw)  a2 

X1  X2 


Likewise,  Z^  is  a  zero  mean  gaussian  random  process  with  the  same  var¬ 
iance.  The  correlation  coefficient  between  the  measured  reverberation 

processes  is  p  : 

Z1Z2 


Z1  V 


E(ZX  Z2) 


E{(X1  +  aX2)  (Xp  +  aX  11 


nv  o 
"l  z'o 


(1  +  a2)  E(X  X2)  +  2a  o2 


(1  +  a  )  p  +  2a 

X1  X2 


1  +  a  +  2ap 


X1  X2 


This  gives  a  relation  between  the  correlation  coefficient  of  the 
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measured  reverberation  r roeess  and  the  actual  reverberat ion  process  as 
a  function  of  a,  the  feedover  separation.  In  this  case  the  measured 
correlation  coefficient  is  greater  than  or  equal  to  the  actual  corre¬ 
lation  coefficient,  and  the  difference  increases  as  the  correlation 
coefficient  decreases.  For  a  feedover  separation  of  UO  dB,  i.e., 
a  =  .01,  the  measured  correlation  coefficient  is  .02  when  the  actual 
correlation  coefficient  is  0.  Thus,  for  a  channel  to  channel  separation 
of  h 0  dB,  the  maximum  difference  occurs  when  the  correlation  is  a  minimum 
and  increases  the  measured  correlation  coefficient  by  .02. 

The  third  source  of  inaccuracy  in  measuring  the  reverberation 
coherence  is  the  limited  size  of  the  sample  ensemble.  An  ensemble  of 
reverberation  returns  in  collected,  and  from  this  ensemble  of  returns 
a  correlation  coefficient  can  be  calculated.  The  standard  deviation 
of  the  correlation  coefficient  as  a  function  of  the  sample  ensemble 


size  N  and  the  correlation  coefficient  p  for  normally  distributed 
17  12 

samples  is: 


The  standard  deviation  thus  increases  with  decreasing  correlation.  For 
a  sample  ensemble  size  of  1000,  the  standard  deviation  is  .03  for  a 
correlation  coefficient  of  0.  Thus,  the  standard  error  introduced 
due  to  the  limited  sample  size  increases  with  decreasing  correlation 
and  is  a  maximum  of  .03. 

If  the  sample  ensemble  size  is  constrained  to  be  1000  samples, 
resulting  in  a  maximum  standard  error  of  .03,  then  a  crude  estimate  of  the 


minimum  reverberation  to  background  ratio  and  the  channel  to  channel 
separation  desired  to  reduce  measurement  errors  to  an  acceptable  level 
can  be  obtained  by  upper  bounding  the  maximum  difference  between  the 
true  and  measured  correlation  coefficients  by  the  maximum  standard 
error  due  to  the  limited  sample  ensemble  sice.  As  has  been  shown,  a 
reverberation  to  bad  ground  ratio  of  20  cLB  and  a  channel  to  channel 
separation  of  Uo  dB  would  maintain  measurement  errors  within  this 
bound. 

As  was  mentioned  earlier,  it  is  also  necessary  to  preserve  the 
relative  phases  of  the  narrowband  acoustic  signals  at  the  face  of  the 
hydrophones  in  order  to  accurately  determine  the  spatial  coherence  . 
Phase  errors  will  occur,  however,  if  the  measurement  system  introduces 
different  phase  shifts  to  the  different  signals.  These  different 
phase  shifts  may  even  be  time-varying.  There  are  many  points  at  which 
the  measurement  system  can  introduce  different  phase  shifts  on  the 
different  signal  channels.  The  hydrophone  elements  themselves  can 
have  different  phase  responses  from  hydrophone  to  hydrophone.  Since 
each  hydrophone  has  its  own  set  of  preamplifiers,  amplifiers,  filters, 
and  other  associated  electronics,  each  set  of  electronics  can  introduce 
a  different  phase  shift  to  that  signal.  When  the  signals  are  recorded 
on  an  analog  tape  recorder,  phase  variations  of  the  tape  recorder  elec¬ 
tronics  and  record  heads  from  channel  to  channel  as  well  as  dynamic 
skew,  static  skew,  flutter  and  tape  stretching  can  all  contribute  to 
phase  errors  when  the  signal  is  reproduced.  Finally,  during  the  analog- 
to-digital  (A/D)  conversion  process,  phase  errors  will  be  introduced  if 
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there  are  any  variations  in  the  A/D  dork  or  A/D  electronics  from 
channel  to  channel.  Thus,  the  measurement  system  must  be  carefully 
designed  so  that  these  phase  errors  can  be  eliminated  or  corrected. 

The  techniques  used  to  maintain  the  phase  integrity  of  the  acoustic 
signals  involve  a  combination  of  electronic  design,  equipment  calibra¬ 
tion,  special  recording  and  digitizing  techniques,  and  processing  of  the 
digitized  signals.  These  techniques  will  be  discussed  in  the  next 
section. 

II . 2  Hardware  Configuration  and  Measurement  Techniques 

The  purpose  of  the  previous  section  was  to  help  clarify  the 
reasoning  that  led  to  the  choice  of  hardware  configuration  and  measure¬ 
ment  techniques  used  to  accurately  measure  the  reverberation  coherence. 
This  section  will  describe  the  measurement  system  used. 

The  receive  system  electronic  equipment  is  divided  into  under¬ 
water  electronics  and  surface  electronics.  The  underwater  electronics 
are  located  within  each  of  the  two  hydrophone  arrays .  Each  element  which 
is  used  in  the  hydrophone  arrays  has  associated  with  it  a  preamplifier 
with  70  dB  of  gain  and  a  differential  line  driver.  The  surface  receive 
electronics  consist  of  a  differential  line  receiver,  gain  and  phase 
adjustable  amplifier,  bandpass  filter,  analog  gate  and  summing  circuit 
for  each  receive  channel.  The  complete  receive  system  block  diagram  is 
shown  in  Figure  II-l. 

The  receive  electronic  equipment  was  designed  to  help  minimize 
the  electronic  noise  level  and  the  feedover  between  chan  '’Is.  Low  noise 
preamplifiers  are  used  which  were  measured  to  have  a  wideband  noise 


15 


level  related  to  the  input  with  the  input  grounded  of  approximately 
-110  dBV.  The  noise  level  in  the  10  kHz  band  about  80  kHz  will  be 
significantly  less  than  this.  A  differential  line  driver/line  receiver 
combination  is  used  with  twisted  pair  cable  to  reduce  the  level  of 

feedover  between  channels  in  the  cable  from  the  arrays  to  the  surface 
electronics.  Each  channel  also  has  a  bandpass  filter  10  kHz  wide  cen¬ 
tered  on  80  kHz  to  reduce  the  noise  level  but  still  provide  sufficient 
bandwidth  for  a  100  us  pulse.  Figure  II-2  shows  a  typical  filter  re¬ 
sponse. 

Measurements  were  made  on  this  system  to  determine  the  electron¬ 
ic  noise  level  and  the  amount  of  feedover  between  channels.  The  elec¬ 
tronic  noise  of  the  receive  electronics  at  the  input  to  the  tape  re¬ 
corder  with  the  surface  amplifier  set  to  0  dB  gain  was  -58  dBV  +  1  dBV. 
The  maximum  signal  level  at  the  tape  recorder  input  is  -6  dBV,  which 
will  give  a  maximum  signal  to  electronic  noise  ratio  of  approximately 
50  dB.  This  is  well  above  the  reverberation  to  electronic  noise  ratio 
needed  to  ensure  that  the  measurement  error  due  to  electronic  noise  is 
below  the  measurement  error  due  to  a  limited  sample  ensemble  size. 

Ambient  acoustic  noise  was  also  measured  immediately  prior  to 
collecting  the  reverberation  data.  The  rms  reverberation  to  rms  ambient 
noise  ratio  was  typically  30  dB  or  better.  Thus  the  presence  of  ambient 
noise  did  not  contribute  significantly  to  measurement  error. 

The  feedover  between  channels  was  measured  and  found  to  be 
primarily  in  the  surface  amplifier.  Channels  which  had  physically 
adjacent  amplifier  circuits  exhibited  approximately  -35  dB  to  -Uo  dB 
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or  feedover,  while  channels  which  were  not  physically  adjacent  had  no 
measureable  feedover. 

The  receive  electronic  equipment  was  also  designed  so  that 
phase  differences  from  channel  to  channel  could  he  corrected  during  the 
digital  processing  of  the  data.  The  topside  amplifier  is  gain  and  phase 
adjustable  so  that,  any  gain  and  phase  differences  up  to  the  inputs  to  the 
analog  tape  recorder  can  be  removed.  To  remove  these  differences,  an 
acoustic  signal  at.  80  kHz  is  placed  in  the  water  in  the  farfield  of  the 
hydrophone  arrays  and  is  picked  up  by  the  arrays  and  measured  at  the  inputs 
to  the  tape  recorder.  After  the  hydrophone  arrays  are  oriented  toward 
tne  source,  the  gain  and  phase  of  each  channel  are  adjusted  so  that  each 
hydrophone  appears  to  be  collinear  with  the  other  hydrophones  in  its  array 
and  has  the  same  amplitude  at  the  inputs  to  the  analog  tape  recorder. 

However,  this  calibration  does  not  correct  for  phase  differences 
due  to  the  tape  recorder  electronics,  reproduction  process,  or  A/D 
process.  To  allow  for  this  correction,  several  signals  are  summed  to¬ 
gether  with  each  channel  at  the  tape  recorder  inputs,  as  can  be  seen 
in  Fig.  TI-1.  A  continuous  sine  wave  at  l60  kHz  is  provided  as  the  A/D 
clock.  The  A/D  is  clocked  at  the  zero  crossings  of  the  sine  wave  to 
provide  sampling  at  320  kHz.  By  summing  the  clock  signal  gith  each 
channel  rather  than  recording  the  clock  signal  on  only  one  channel  or 
not  recording  a  clock  signal  at  all,  any  phase  variations  due  to  skew, 
flutter,  or  tape  stretching  can  be  minimized.  During  the  reproduction 
and  conversion  process,  each  channel  is  filtered  to  separate  the  80  kHz 
signal  and  the  l60  kHz  clock,  and  each  channel  is  sampled  with  i4 s  own 


clock  signal. 


This  technique  alone,  though,  will  not  correct  for  phase  varia¬ 
tions  of  the  tape  recorder  electronics  and  record  heads  from  channel  to 
channel,  since  the  phase  response  of  the  tape  recorder  electronics  and 
recorder  heads  is  not  the  same  at  80  kHz  and  160  kHz  and  the  difference 
in  phase  response  changes  from  channel  to  channel.  Thus  an  80  kHz 
pulsed  sinusoid  is  summed  with  each  channel  to  provide  a  reference 
pulse  for  each  ping.  The  reference  pulse  is  digitized  as  a  part  of 
the  reverberation  data  for  each  channel  and  ping.  The  phase  of  the 
reference  pulse  is  then  computed  digitally  and  compared  to  the  phase  of 
the  reference  pulse  from  an  arbitrary  reference  channel  and  ping.  The 
digitized  reverberation  signal  is  then  digitally  phase  shifted  by  an 
amount  which  is  the  difference  between  the  phase  of  these  two  reference 
pulses.  This  is  performed  for  every  channel  and  ping. 

However,  this  technique  is  unambiguous  only  when  the  phase 
difference  between  two  channels  is  less  then  l80°.  This  will  not 

always  be  the  case.  To  compensate  for  large  phase  differences,  an 
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80  kHz  phase-encoded  pulse  (Barker  pulse)  is  also  summed  with  each 

channel  on  each  ping  and  digitized  with  the  reverberation  data.  The 

Barker  pulse  has  a  particular  phase-encoding  that  produces  a  very 

narrow  cross-correlation  peak.  Thus,  the  digitized  Barker  pulse  is 

digitally  cross-correlated  for  each  channel  and  ping  with  the  Barker 

pulse  of  the  arbitrary  reference  channel  and  ping  for  several  different 
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time  shifts.  The  digitized  reverberation  signal  and  reference  pulse 
are  then  shifted  by  the  number  of  samples  necessary  to  produce  the 
largest  cross-correlation  value.  This  is  performed  for  every  channel  and 

Since  the  data  are  digitized  at.  four  times 
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ping  prior  to  phase  shifting. 


the  center  frequency,  each  sample  represents  90°  of  phase  shift.  Thus, 
the  cross-correlation  of  the  Barker  pulses  allows  all  the  reverberation 
data  and  reference  pulses  to  be  aligned  to  within  +  90°.  The  reference 
pulses  will  then  provide  unambiguous  measurements  of  the  phase  differ¬ 
ences.  It  should  be  noted  that  by  recording  the  reference  and  Barker 
pulses  on  every  ping  as  well  as  every  channel,  phase  variations  in  time 
are  also  corrected. 

To  minimize  phase  variation  between  channels  due  to  the  A/D 
equipment,  the  channels  are  digitized  separately  so  that  each  channel 
will  be  processed  through  the  same  set  of  equipment.  A  2h0  kHz  pulsed 
sinusoid  is  summed  with  the  reverberation  signal  and  is  filtered  upon 
reproduction  to  separate  the  2h0  kHz  pulse  from  the  reverberation 
signal.  This  pulse  is  used  as  a  sync  pulse  to  start  the  A/D  on  each 
ping. 

The  receive  electronics  also  contains  an  analog  gate  on  each 
channel  which  gates  off  the  reverberation  signal  while  the  sync,  Barker, 
and  reference  pulses  occur  to  facilitate  separating  them  from  the  rever¬ 
beration  signal. 

II. 3  Transducer  Assembly 

The  transducers  consist  of  one  projector  and  two  line  receiving 
arrays.  Each  receiving  array  has  h2  elements.  The  two  receiving 
arrays  are  oriented  at  right  angles  to  each  other.  Figure  II-3  shows 
the  relative  orientation  of  the  projector  and  receiving  arrays,  and  also 
the  spacing  of  the  elements  within  the  receiving  arrays.  Nine  elements 


are  used  in  the  vertical  array  and  four  in  the  horizontal  array.  Each 
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FIGURE  11-3 

TRANSDUCER  ARRAY  DIMENSIONS 
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element,  is  a  planar  reet  angul  nr  element  sc’.  t  y 

ments  of  the  vertical  array  have  their  f  ar  d  'mens !  -i: 

while  the  elements  of  the  horizontal  array  have  the;  r  ' 

vertically.  The  faces  of  the  elements  in  each  array  •  i  . 

the  elements  in  that  array .  Directivity  ra*  t  orm  r  si:.--  •• 

element  are  shown  in  Firs.  T 1-4  and  TT-r'.  The  :  at '  <•>•:..•  -  v-  :  •’ 

75,  80 ,  and  85  kHz,  and  are  typical  of  the  ra*  :•  rt...  r  .  .  •.  i rtv  • 

The  receive  sensitivity  of  each  receiver  channel  was  r.eac.r*  :  s 
and  found  to  be  approximately  -107  dB  re  1  V/yFa  with  #  in  cad.  it.  • 
topside  receiver  electronics.. 

The  projector  is  a  rectangular  planar  array  7.0.?  cm  by  •?•.. 

Its  lirectivit.y  patterns  are  shown  in  Figs.  TI-o  and  11-7  t--  tree,  .enclos 
of  7%  80,  and  85  kHz.  The  project i ny  response  was  mens-i.ro  :  a*  k H  • 
and  found  to  be  approximately  171  dr  re  1  viTa/V  at  ’■  :::. 

The  projector  and.  receiving  arrays  are  mounted  in  a  -id  ■' -v>  - 
work  and  buoyed  by  flotation.  The  framework  is  at  t  ached  to  moors 
which  allow  it.  to  be  tilted  and  rotated  remotely,  fynohro  t  r.ansmi  i  ten¬ 
on  the  tilt  and  rotate  motors  allow  remote  monitoring  of  the  oriental  im 
of  the  arrays.  Figure  TT-8  is  a  photograph  of  the  assembly  mounted  ‘ o 
a  tower  section  in  air.  The  entire  assembly  was.  deployed  on  ;u-,  und.er- 
wat.er  tower  in  Lake  Travis  near  the  ARI.tUT  Lake  Travis.  Test  .’tat  i. mi 
(Fig.  TT-9) .  ^  The  tower  is  located  in  approximately  '*0  J;  if  vat  ••>!•.  "*h 

center  of  the  transducer  assembly  was  approximately  !  1.  '  r:  below  t  ho 


surface . 
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FIGURE  11-4 

DIRECTIVITY  PATTERNS  FOR  A  SINGLE  RECEIVER  ELEMENT 
AT  75,  80,  AND  85  kHz,  MEASURED  ABOUT  THE  LONG  DIMENSION 
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FIGURE  11-5 

DIRECTIVITY  PATTERNS  FOR  A  SINGLE  RECEIVER  ELEMENT 
AT  75,  80,  AND  85  kHz,  MEASURED  ABOUT  THE  SHORT  DIMENSION 
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FIGURE  11-6 

PROJECTOR  VERTICAL  DIRECTIVITY  PATTERNS 


.1  I  .  > : 


i  onditlons  for  Data  Collection 


A  little  over  1000  reverberation  returns  were  recorded  in  an]  rox 
imately  lU  minutes.  A  pulsed  CW  signal  was  transmitted  at  a  center 
frequency  of  80.0  kHz.  The  pulse  length  was  100  ms.  The  entire  array 
was  tilted  to  produce  a  grazing  angle  with  the  surface  .f  approximately 
10.5  degrees. 

The  water  was  isothermal  to  a  depth  cf  at  least  . n .  n;  a 
temperature  of  13.8°C.  The  sound  speed  in  the  water  was  determined  t  c> 
be  approximately  1U63  m/s.  The  temperature  of  t  h<-  air  was  lh.V'O  and 
winds  were  from  the  north  at  35-50  kph.  The  wave  heights  were  estimated 
to  be  .  1  to  .6  m.  The  array  was  pointed  approximate  ly  1»T°  east  of  north 


_  J 


id.  THEORETICAL  MWEI. 


There  are  two  fundamental,  approaches  to  the  theory  of 

backscattering  from  a  randomly  rough  surface.  The  "physical"  model  seeks 

the  solution  of  a  scalar  wave  equation  for  statistically  random  boundary 

conditions.  It  considers  scattering  from  a  rough  surface  in  terms  of  the 

statistical  distribution  of  surface  wave  heights.  The  other  approach  is 

the  "quasiphenomenological”  model,  developed  most  generally  by  Middleton, 

which  treats  reverberation  as  weak  scattering  from  random  point  sources 

representing  inhomogeneities  in  an  otherwise  homogeneous  medium.  It 

thus  avoids  the  probLem  of  complicated  boundary  conditions,  but  instead 

requires  a  knowledge  of  the  interaction  of  the  incident  sound  with  the 

point  sources  in  term:;  i.f  .in  impulse  response  function  fo'  the  point 

sources.  It  is  generally  agreed  that  Middleton's  "quasiphenomenological" 

model  is  the  most  complete  model  of  reverberation^1  but  suffers  from 

the  difficulty  of  having  to  specify  a  scattering  function.  Fortunately, 

it  has  been  shown  that,  at  least  in  some  cases,  a  very  simple  scattering 

17 

function  works  well. 

The  remainder  of  this  chapter  is  organized  as  follows.  First 
the  development  of  Middleton's  model  as  given  in  Ref.  27  is  outlined. 
Additional  details  of  this  development  are  given  in  Appendix  A.  Special 
attention  is  given  to  the  assumptions  inherent  in  the  model  ,  and  the 
results  for  the  covariance  function  are  shown.  Next  the  simplifications 
which  lead  to  the  specific  form  of  the  model  used  in  the  present  study 
are  discussed.*  This  includes  specifying  the  scattering  function.  This 

*This  same  form  of  the  model  was  used  by  Frazer  (Ref.  17). 
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1 


scattering  function  is  used  in  Appendix  B  to  develop  an  expression  for 


what  Middle-ton  refers  to  as  the  elementary  scattered  waveform.  An  appro 
mation  of  the  covariance  function  which  is  valid  for  bandlimitel  signals 
is  then  developed,  and  some  additional  approximations  which  facilitate  t 
numerical  computation  of  the  covariance  function  are  discussed.  Details 
of  the  computations  of  the  directivity  patterns  are  pi veil  in  Appendix  2. 
i T I . 1  The  General  Model  as  Developed  by  Middleton 

The  following  is  a  summary  of  the  main  features  of  Middleton's 

model : 

1.  The  total  scattering  process  is  a  linear  superposition  of 
scattering  from  independent  (in  space  and  tine)  random  point 
sources  representing  inhomogeneities  in  an  otherwise  homo¬ 
geneous  medium.  This  is  sufficient  to  specify  the  scatter¬ 
ing  distribution  as  Poisson. 

2.  The  interaction  of  sound  with  the  scattc-rers  is  represented 
by  an  impulse  response  function. 

i.  General  geometries  and  transmit  and  receive  apertures  are 
allowed. 

k.  General  illuminating  signals  are  allowed. 

Reverberation  from  multiple  receivers  is  considered. 

b.  Shadowing  effects  for  rough  surfaces  are  included. 

7.  Doppler  effects  due  to  motion  of  the  source,  receiver,  and 
scatterers  are  included. 

In  Ref.  27  a  characteristic  function  for  the  reverberation  pro¬ 
cess  is  derived  based  on  Poisson  statistics.  From  the  character! si i c 
function  the  various  moments  can  be  determined  by  taking  the  appropriate 
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derivatives.  An  abbreviated  explanation  of  this  derivation  is  given 
in  Appendix  A.  The  characteristic  function  depends  upon  a  superposition 
of  elementary  scattered  waveforms,  each  representing  the  scattered 
signal  from  a  point  source.  These  waveforms  are  determined  from  wave 
and  ray  theory  in  a  homogeneous  medium  and  include  the  scattering 
impulse  response  function.  The  following  assumptions  are  made  in  this 
development : 

1.  Scattering  from  the  inhomogeneities  is  sufficiently  weak 
that  only  primary  scattering  is  considered,  i.e.,  multiple 
scattering  is  ignored. 

The  probability  that  exactly  N  scattering  events  occur  in 
a  sufficiently  small  interval  is  given  by  the  Poisson 
probability  1  aw. 

3.  The  elementary  scattered  waveform  for  each  scatt.erer  is  not 
a  function  of  the  specific  scatterer;  that  is,  the  func¬ 
tional  form  of  elementary  scattered  waveforms  are  identical 
for  all  scatterers.  Also  any  random  parameters  associated 
with  each  scattering  event  have  the  same  distribution. 

This  implies  that  only  one  impulse  response  function  is 
sufficient  to  describe  the  scattering  process. 

The  scattering  zone  is  in  the  farfield  of  the  projector  and 
receivers . 

With  these  assumptions  the  Joint  (in  time)  characteristic  func¬ 
tion  of  the  reverberation  process  is  shown  in  Ref.  27  to  be  given  by: 


i 


F„(1h-  tli 


At,  ,  t  t',  A) 
n  n 1 


=  exp/^p(X-,t')<exp  i  { ^  !  C  e  U'(t„  ;X  ,d  I  }-i.)0d> 


whore 

t,  ...  .  ,t  are  the  times  of  observation; 

1  n 

t'  represents  the  temporal  dependence  of  the  .-.’at  t  or  density; 

A  is  the  scattering  region; 

X  represents  a  point  in  the  scattering  region; 

P  is  the  density  of  scatterers  at  point.  \  and  time  t. '; 

U  is  the  elementary  scattered  waveform  at  time  t.,  due  to 
scatterers  at  point  X; 

0  represents  a  set  of  random  parameters  associated  with  the 
scattering  process; 

and 

<  >  represents  the  exnoctei  value  with  re. -feet  t  o  tic  rani  m: 

par a.met ei-s  0. 

The  various  moments  of  the  reverberation  process  are  given  by  the  appro¬ 
priate  derivatives  of  the  characteristic  function.  i-'or  instance,  if 
X(t,t")  represents  the  reverberation  process  ;<t  time  t.  due  ton  rattier 
density  at  time  t  then  the  moments  are  river  tv.  tin  f  1  lowing; 

<X(t,tl)  =  /Ap(X;t  •)<U(t;\,0)>o  dX 

<X2(t,tO)  =  /Ap(X;t')<U2(t;X,0)>od\+<v(t  ,t  M>- 

<X(t1,t')X(t?,t')>  =  /AptX;t')(n(ti;\,ol'nt  d\ 

+  <x(t3  ,t'))-(x(tp,f)> 


Si 


The  I'V'var ianee  of’  rho  proct ss  is  then: 

Kyi.  trt2|tO  x(t1,f)x(t.0,t')>-<'x(t]  ,t')><X(t2,t')> 

=  /Ap(A;t.  )<!/{  tj  ;X,e)V(tn;\^))Qd\ 

For  spatially  joint  processes ,  for  exampLe  a  single  projector 
and  two  receivers,  the  joint  characteristic  function  of  the  two  random 
processes  and  X^“^(t,,t')  is  a  modification  of  the  charac¬ 

teristic  function  given  above: 

FJDx(2)  =  eX*'{V  ^.t'KexpFi^Up+iC^l-l)  ^dXJ  . 

Some  of  the  relevant  moments  are  given  by: 

K1  ,tp|t  >  -  (  X-j  ( t- 1  ,  t ' )  X]  ( t ,  t.  ‘  )> 

-=  /AP1(A,t')<ir1(t);\,e1)Ui(t2;A,e1)>edA 

K(,?)(t1,t,,|t')  =  (x,(trt'iyt2,ti> 

=  Pp(  A  S  -, (  f  i  >  A  ) U^(  t0  ;  A  ,0O  )^pdA 

K(ir>)(tltt0|t*)  =  <X1(t1  ,t')X0(t?,t')>-(X]  (t]  ,t')><X(0(tp,t')> 

XII. 2  The  Theoretical  Model  as  Used  in  the  Present  Study 

A  snecial  case  of  this  model  is  now  considered.  Two  essential 
assumptions  are  invoked.  rt.  will  be  shown  later  that  the  covariance  can 
be  written  as  the  sum  of  a  time  difference  component  (depending  on  the 
time  difference  tp-t^ )  and  a  time  sum  component  (depending  on  the  time 
sum  t^  +  t  ) .  It  is  assumed  that  the  time  sum  component,  is  negligible  in 
comparison  to  the  time  difference  component.  Thus  only  the  expression 


for  the  time  difference  component,  of  the  covariance  is  developed.  This, 


assumption  is  reasonable  when  the  transmitted  signal  is  narrowbanu, 

i.e.,  the  center  frequency  of  the  transmitted  signal  is  much  larger  than 
i4-  bandwidth.  Secondly,  the  scatterer  is  modeled  as  a  perfect  point 
reflector,  i.e.,  scattering  is  omnidirectional,  frequency  insensitive 
(within  the  bandwidth  of  the  signal),  location  insensitive,  and  introduce 
no  time  delay.  Since  much  of  the  scattering  is  from  bubbles  at  vr  near 
the  surface  which  have  diameters  much  smaller  than  the  wavelenct of  the 
transmitted  sound,  this  assumption  is  a  reasonable  one.  Other  assump¬ 
tions  are  as  follows: 

1.  No  random  parameters  are  associated  with  the  scattering 
event,  i.e.,  there  is  no  dependence  on  0. 

2.  The  directivity  patterns  of  the  projector  and  receivers 
are  frequency  insensitive  within  the  baniwidth  of  the 
signal  (narrowband  signal). 

3.  p(X,t')dX  =  o  df,  where  o  is  the  surface  scatterer  ienri’v 

s  s 

and  is  constant,  independent  of  position  or  time,  and 
where  the  region  of  integration  is  over  the  surface  of  the 

water. 

a.  There  is  no  doppler  of  the  projector,  receivers,  s' cat - 

terers . 

5 .  There  are  no  shadowing  effects . 

Using  these  assumptions,  it  is  shown  in  Appendix  B  that  the 
elementary  scattered  waveform  can  be  written  as 


wh°  re 


M  )  ■=  /  ‘'I'  (  i  )t  ,(  r)r; .  ( rk1" 1  rdf 

I*  ~ "oo  ri  X  111 


and 


t  -  t  - 
r  i' 


T,Jf)  is  the  frequency  response  of  the  receiver; 
K 


T  ( f )  is  the  frequency  response  of  the  projector; 


(f)  is  the  frequency  spectrum  of  the  input  electrical  signal 


Rp  is  the  distance  from  the  receiver  to  the  point  scatterer  at 
the  coordinate  A; 


R,p  is  the  distance  from  t.he  transmitter  to  the  point  scatterer 


at.  the  con  ruinate  A; 

c  is  the  speed  of  sound; 

Ap(R^)  is  the  directivity  function  of  the  receiver  in  the 
direction  R^; 

A^R^)  is  t.he  directivity  function  of  the  projector  in  the 


direction  R^; 


Y  represents  the  response  function  of  a  scatterer  which  intro¬ 
duces  no  frequency  shift  or  time  delay  to  the  scattered 
waveform. 

Using  tills  form  of  the  elementary  scattered  waveform,  the  co- 
variance  of  a  joint-  process  can  now  be  written  as 


K(1P)  (1-,  'A  =  «SJA  VVX)l,2(t2»X)dr’ 


YlYf°s  r  AR1(RR1)AR2^R2^V|V_ 
j  Rt  -R  « 


(lmT  "12  “Rl  RR2“T 


X  v^t^Jv^U^Jds 


where 


t  =  t  - 
rl  1 


t,-P  ~  to 


+R__ 
Rl  T 


RR2+RT 


This  surface  integral  can  now  be  evaluated  numerically  to  compute  *he 
covariance . 

This  expression  for  the  covariance  can  be  expanded  further  by- 

considering  the  product  V^(t  ^)Vn(t^).  Representing  V(t  '  ir.  a 

ture  formulation  as 

V(t  )  =  X  cosoj  t  +  Y  sinw  1 
r  or  or 

where  u)q  is  the  angular  center  frequency  of  the  projected  signal,  allov 

the  product  V, (t  , )V„(t  to  be  written  as 
1  rl  2  r2 


x  (X2(trt>)c0“‘otr-’+!,-'(tr  0)sinw^tr-J 


Multiplying  these  two  expressions  together  and  utilising  the  trigono¬ 
metric  identities  for  the  product  of  sine  and  cosine  functions  gives 
V,V  =  ^[X  X0+Y  Yjcosw  (  t-AP)  +  ^[X  Y  -Y  X.^lsinu)  (  r-AF’' 

+  !5tXlX2"YlY2]c0sa,o  (Vl+tr2) 

+  yXlVYlX2]sin“o  (trl+tr2) 


where 


AR 


RR2-R 


Rl 


c 

and  the  dependence  of  the  quadrature  components  on  t  ,  t 

rl  r2 

dropped  for  compactness  of  notation.  Applying  trigonometric 


has  been 
ident i ties 


for  the  cosine  and  sine  of  a  sum  to  this  express!  .n  and  rearranging 
terms  Rives 


V1V.,  =  hi  [X-jX^+Y^Y^lcosu)  .AR  -  [  X  Y  -Y  .  X j  s  i  nu._  AR  jeesu 

+  y[X.Y.-Y  Xjconui  AR  +  fX.X,+Y  Y  jsinu  AR'isin* 

+  I*  [ X,  X  -Y  Y  ,  jcosiij  (t  +t  J 
1  .  1  i-1  c  rl  !• 


+  V  x  Y  t  Y  X  .  |  e  i  s;j  (I  .  +t. 

J  .  !  .  ft  )•.' 


Substituting  this  expression  for  v  V  into  the-  expression  for  the  a 

results  in  a  time  difference  component  (terms  containing  cosw  t  ana 

si  nor  t  )  and  n  time  sum  component  (terms  containin'  cos...  ( t  .  +*  .  1  a: 
o  rl  r. 

sinu  (t  , +t  j).  tor  a  narrowband  si  anal  the  time  sum  terms  are  rat 
o  rl  re 


■  :  .11  v 


varying  functions  over  the  area  of  integration  for  the  covariance,  since 

t.  ,  +t  „  contain  the  variables  R  ,  and  R„  ,  which  vary  -ver  the  area 
rl  rP.  T’  Rl  R. 

of  integral  ion.  Thus  t  h-»  time  sum  component  will  average  <  .  approx  ir.at  el; 

•1  h 

zero  in  comparison  to  the  time  differer.ee  component  of  the  eovar  . 

Therefore  in  the  expression  for  the  covariance  it  is  roa.  na: ■:<  ‘  ..  •  •  rr<  . 

the  product  of  and  V,  as 


where 


V .  ( t  )V_(t  )  =  Xcosio  i+Ysinm  i 

1  rl  2  r.2  o 


Rr  '"Rd - 

X  =  yfX.,(t  .)X0(t  )  +Y  ( t  ,)  ,tt  )]cnru„ 

1  vi  ■  r.  1  rl  .  n  r 

-  [  Xi ( t  r  i  )  Y;J  ( t  )  -  Y 1 1  tr a  )  Xp ( t rp ) ]  s  i  nu.  Rl  )  ) 


Vr.lsL) 


Y  =  -if  [  X^  ( )  Y0(  t rp)-^  j  ( ^  r  l  )Xp(  trp )  lcosii)^(  f. 

R  -R 

+  [X,(t  )X0(t  )+y  (t.  ) Y  ( t.  )  ]sinu)  (  --  ^  )  1 

1  rl  ■  r.  l  rl  2  r.  o  c 


oubsti  tutiiiK  thin  expression  t'or  the  product  of  V  arid  Vo  into  the 

( li' ) 

equation  ‘'nr  the  covariance  K  '  (t  ,t  )  Kives  the  difference  com- 

J.  L 

ponent 

r  / .  )  _  rYlY2°s  r  ARl^Rl^AR?(RRP)AT(RT^ 

Ki^it  ,i)  ~  i  ~T  J  .  - o -  XdM  costo  i 

#i2  RrA24 

,  ,VA  ,  VVVV>T(V  . 

+  1  '  r  J.  - ~ -  YdSJsinoi  T 

(WT  12  R  R  R2  0 

R1  R2  T 


who  re 


X12(tl’T) 

YlVs 

~  - T 

( hn  ) 

Y,Yo0 

Y12(tl,T' 

12s 
~  ,  K  . 

( Urr ) 

X12(tl’T) 

and  Y  0(  t 
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\iVt 


difference  component  of  the  covariance.  As  such  the  difference  com¬ 
ponent  of  the  covariance  can  be  expressed  in  terms  of  an  envelope 
function  and  a  phase  function: 

Rl°(  tl ,T  )  =  E12^tq>T^C0S^w0T+<J’]_2^tl,T^  ’ 


where 


El2(tl’x)  =VXi2(tl’T)  +  Y12(tl’T) 


VVt)  =  tan  (  X12(t1;x) 
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Implementation  of  thin  model  on  i  computer  is  facilitated  by 
several  additional  approximations: 

1.  The  function  V(t  )  is  approximated  from  experimental 
measurement s . 

The  surface  integral  is  evaluated  numerically,  using  the 
Simpson's  integration  technique. 

i.  The  directivity  functions  of  the  projector  and  receivers 
are  computed  from  the  theoretical  response  for  planar 
arrays . 

The  function  v(t.  )  represents  the  input  electrical  signal  after 
it  has  been  fi  li.ered  by  the  projector  and  receiver.  It  can  be  measured 
exneri mental ]y  for  the  particular  projector,  receiver,  and  signal  used 
by  projecting  at  an  omnidirectional ,  frequency  insensitive  (within  the 
narrowband  assumpi i  1)  target  and  recording  the  received  echo.  This 
echo  is  used  to  renresent  V(t.  ),  within  an  amplitude  scale  factor.  The 
onset  of  the  echo,  which  would  represent  the  time  t  =0,  cannot  be  pre¬ 
cisely  defined,  since  the  exact  distance  between  the  projector  and 
target  is  generally  difficult  to  determine  to  the  required  degree  of 
accuracy.  Thus  a  reasonable  point  on  the  received  echo  is  picked  and 
defined  as  the  point  t  =0.  Quadrature  samples  of  the  received  echo  are 

j  n 

obtained  by  directly  sampling  V(t  ),  ’  and  the  quadrature  components  of 
V(t  )  arc  computed  for  any  required  time  t^  by  linear  interpolation  on 
the  quadrature  samples  of  the  received  echo.  For  t.  less  than  zero  or 
greater  than  1  h<-  extent  of  the  echo,  V(t  )  is  defined  as  zero. 


The  waveform  V(t  )  was  measured  experimentally  by  transmitting 
a  100  us  CW  pulse  centered  at  80  kHz  at  a  fluid-filled  spherical  lens 
with  a  diameter  of  15*2  cm  and  recording  the  echo  received  by  the 
arrays  from  the  sphere.  The  sphere  was  at  the  same  depth  as  the  arrays 
and  at  a  range  of  approximately  100  m.  The  envelope  and  phase  of  the 
received  echo  on  each  receiving  hydrophone  is  given  in  Figs,  I1I-L  and 
III-2 . 

The  surface  integrals  are  evaluated  numerically  using  the 

Simpson's  integration  technique.  This  integration  technique  evaluates 

an  integral  by  summing  over  fixed  step  sizes.  The  integral  of  a 

function  f(x)  with  respect  to  the  variable  x  between  the  limits  a  and  b 

is  approximated  by 

b  .  N+l 

/  f(x)dx  =  -y-  mf ( a+( n-1 ) Ax) 

a  n-1 

where 

Ax  is  the  chosen  step  size; 

N  is  the  integer  closest  to  7—^ — ; 

Ax  ’ 

and 

!1,  if  n  =  1  or  N  +  1 
2,  if  n  is  even 
** ,  if  n  is  odd. 

The  step  size  Ax  is  chosen  to  provide  the  desired  accuracy  for  the  inte¬ 
gration  . 

To  evaluate  the  surface  integral  required  for  the  covariance 


computation,  the  integration  is  performed  over  the  variables  t. '  and  , 
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FIGURE  III -2 

THE  AMPLITUDES  a;  AND  PHASES  4>-y  OF  THE  WAVEFORMS  Vj(tf) 
FOR  THE  HORIZONTAL  ARRAY 
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with  the  transformation 


dS  ►  vc  tTdtTd4> 

where  $  is  the  azimuthal  angle  to  a  point  scatterer,  and  t  is  the  time 
of  the  scattering  event  relative  to  the  beginning  of  the  transmission. 
The  transformation  factor  w  is  a  function  of  both  P  and  t_ ,  and  is 
given  in  explicit  form  in  Appendix  C.  For  the  parameters  of  this 
study,  w  was  very  close  to  1.  Si  nee  =  ct  ,  tic  n 


dF  =  f  -  dt  .Lp 
ZT  1 


The  Simpson's  integration  technique  is  used  to  evaluate  both  integrals, 
with  the  t^  integration  performed  first  since  the  limits  of  the  tT 
integration  turn  out  to  be  <J>  dependent. 

The  quadrature  components  of  the  covariance  are  then 
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where  the  area  of  common  scat.terers  A^0  determines  the  limits  on  the  $ 
and  t^  integrations.  Figure  IIT-3  shows  the  common  scattering  region 
for  two  horizontal  receivers  and  a  projector.  Only  scat.terers  whies  were 
illuminated  between  the  times  t  and  t  contribute  to  the  waveform 

-L  Ij  L  U 

of  receiver  1  at  time  t^  ,  where  t.  and  t  are  given  by 


[■ 


i  J 
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where  T  is  the  transmit. ted  pulse  length.  Likewise,  only  the  scatterers 
which  were  illuminated  between  the  times  tp^  and  tp^  contribute  to  the 
waveform  of  receiver  2  at  time  1 =  t  +  i,  where  t  and  tOJJ  are  given 


hv 


*2L  "  V-Jh— 


t2U  =  V1- 


RR2(  t'2U) 


Thus  the  scatterers  which  contribute  to  both  receivers  at  both  times  t^ 
and  t0  are  those  which  were  illuminated  between  the  times  t  and  t  , 

c.  L  U 

where  t  and  t  are  given  by 
t  =  maxlt^  ,  t,?IJ 


tu  =  min[tU!,  t?[JJ 


The  above  equations  for  t  ,  t  ,  t  ,  and  t  are  solved  iteratively 

i  1  j  ill  c  L  rll 

i.n  order  to  determine  the  limits  for  the  t  integration. 

The  limits  on  t.he  <p  integration  are  much  easier  t.o  determine, 
and  in  fact  depend  only  on  the  horizontal  directivity  pattern  of  the 
projector.  The  limits  are  chosen  so  that  most  of  the  projected  sound 
is  within  the  chosen  limits.  If  the  sidelobes  of  the  projector  are 


sufficiently  low,  it  can  be  assumed  that  to  a  good  approximation  the 
scatterers  are  illuminated  only  by  t.he  main  lpbe  of  the  projector.  The 
limits  are  then  typically  chosen  to  include  only  the  main  lobe. 


uc 


The  directivity  functions  used  in  the  covariance  calculation 
are  computed  from  the  theoretical  response  for  planar,  shaded  elements. 
The  projector  is  bizonally  shaded  in  both  dimensions,  while  the  receiving 
elements  are  bizonally  shaded  in  the  long  dimension  and  unshaded  in  the 
short  dimension.  The  projector's  directivity  function  is  given  by 


a(rt) 


sinX+'^sin^) 

5X 

X 


sinY+Jssin( 

5Y 

X 


where 


X  =  2c  *TVsinaTCOs6T  ’ 


Y  =  2c  £THsinaTsineT  ’ 

R^,  is  a  unit  vector  from  the  projector  to  the  point  scatterer 
with  respect  to  a  coordinate  system  centered  on  the 
projector,  as  illustrated  in  Fig.  1 1 1— U  ; 
is  the  vertical  length  of  the  projector; 
is  the  horizontal  length  of  the  projector;  and 

in 

aT  and  B^  are  spherical  angular  coordinates,  represented  by  the 
angles  a  and  8  in  Fig.  III-1+. 

Likewise,  the  directivity  function  for  the  horizontal  receivers  is 
given  by 


*(eR1)  • 


slpX-Hisln(f)  .  ^ 

5X  *  Y 

It 


and  for  the  vertical  receivers  is  given  by 


A("Ri> ' 
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FIGURE  111-4 

COORDINATE  SYSTEM  CENTERED  AT  AN  ELEMENT  OF  THE  ARRAY 
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where 


X  =  —  i.  .  sina  . cos8„ . 
2c  RiV  Ri  Ri 


Y  =  SiH3inaRisinBRi 


.th 


Rp^  is  a  unit  vector  from  the  i  receiver  to  the  point  scatterer 
with  respect  to  a  coordinate  system  centered  on  the  receiver; 

9,  is  the  vertical  length  of  the  i^  receiver; 

Rl  V 

is  the  horizontal  length  of  the  i'th  receiver;  and 

HlH 

a  .  and  6_.  are  spherical  angular  coordinates  represented  by  the 

nl  nl 

angles  a  and  B  in  Fig.  III-U. 

Details  of  the  computation  of  the  directivity  functions  in  terms  of  the 
variables  of  integration  t^  and  <p  are  given  in  Appendix  C. 

A  comparison  of  the  computed  and  measured  directivity  patterns 
for  the  projector  is  given  in  Figs.  111-5  arid  1 1  [-6  for  a  frequency 
of  80  kHz.  The  corresponding  directivity  patterns  for  a  typical  re¬ 
ceiver  element  are  given  in  Figs.  III-7  and  III-8. 

The  steps  sizes  for  the  tj,  and  $  integrations  were  chosen  to 
achieve  an  accuracy  of  at  least  one  percent.  A  step  size  of  2.5  x  10  ^  s 
was  found  to  be  adequate  for  the  t^,  integration.  The  required  step 
size  for  the  <p  integration  was  found  to  depend  upon  whether  the  elements 
were  separated  vertically  or  horizontally.  The  integrand  of  the  $  inte¬ 
gration  for  vertically  separated  elements  changed  only  slowly  over  the 
range  of  <f>.  Thus  a  step  size  of  .?  degrees  was  sufficient.  However, 
for  horizontal  separations  the  int.egrand  changed  more  rapidly.  For 
separations  out  to  5  cm,  a  step  size  of  1  degree  was  sufficient,  but 
for  a  separation  of  7-5  cm,  a  0.5  degree  step  size  was  necessary. 
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FIGURE  111-5 

COMPARISON  OF  COMPUTED  AND  MEASURED 
PROJECTOR  VERTICAL  BEAM  PATTERN 

- MEASURED  PATTERN 

-  COMPUTED  PATTERN 

ARL  UT 
ASW613 
GRW  GA 
1  23  80 
REV3  3-81 


FIGURE  111-6 

COMPARISON  OF  COMPUTED  AND  MEASURED 
PROJECTOR  HORIZONTAL  BEAM  PATTERN 


MEASURED  PATTERN 
COMPUTED  PATTERN 
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FIGURE  III  7 

COMPARISON  OF  COMPUTED  AND  MEASURED  RECEIVER 
BEAM  PATTERN  IN  THE  LONG  DIMENSION  AT  80  kHz 

- - MEASURED  PATTERN 

-  COMPUTED  PATTERN 
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Further  spacings  required  such  smal  i  step  sizes  that  in  view  of  the 
negligible  levels  of  coherence  between  these  spacings,  the  theoretical 
computations  were  not  performed.  The  b  integration  was  performed  over 
the  limits  from  -15  to  15  degrees. 


IV. 


STATISTICAL  TESTING  OF  SAMPLE  DATA 


The  reverberation  field  is  campled  to  form  experimental  ensem¬ 
bles,  and  various  statistical  tests  are  performed  on  the  ensembles  to 
validate  the  sample  data.  This  chapter  discusses  the  sampling  procedure, 
including  the  assumptions  inherent  in  the  sampling  process.  The  statis¬ 
tical  tests  are  also  discussed,  and  the  results  from  applying  these 
tests  to  the  sample  data  are  presented. 

IV. 1  Formation  of  Sample  Ensembles 

The  electrical  signals  produced  at  the  output  of  each  receiving 
element  by  the  acoustic  reverberation  are  treated  in  this  study  as 
continuous  time  random  processes,  denoted  by  ( t  ,u, ) ,  where  i  indexes 

the  particular  random  process  (l<i<L3),  t  represents  the  observation 
time  relative  to  the  beginning  of  the  pulsed  transmission,  and  w  indexes 
the  particular  realization  of  the  random  process  which  was  generated 
from  the  parent  population  fi  of  all  possible  realizations  of  the  random 
process.  For  this  study,  1<j><1081.  For  a  fixed  m,  v.  (t,w)  is  a  partic¬ 
ular  realization,  or  sample  function,  of  the  i*'  random  process. 

The  continuous  time  random  process  v  (t,ui)  is  sampled  in  time 
during  the  analog-to-digital  (A/D)  processing  for  each  realization.  The 

sampled  function  is  denoted  bv  v.(t  where  t  denotes  one  of  the 

*  "in  n 

discrete  times  at  which  the  sampling  occurred.  For  this  study,  l<n<^200, 
where  the  sampling  occurred  at  four  times  the  center  frequency  of  the 
transmission,  or  320  kHz.  Thus  the  reverberation  process  was  sampled 
for  a  total  time  of  10  ms. 

SU 
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For  a  fixed  t  and  receiver  element  i,  the  collection  of  all 
n 

v^(t^,to),  l<u<108U  represents  a  finite  sample  set  of  a  random  variable, 

and  is  referred  to  as  a  sample  ensemble.  To  be  a  valid  sample  ensemble, 

the  sampling  assumptions  must  be  met:  each  sample  of  the  ensemble  must 

be  independent  and  identically  distributed,  i.e.,  the  samples  are  a 

random  collection  of  realizations  of  the  random  variable  v.(t  ,00)  for 

1  n 

a  fixed  t  and  i .  The  results  of  testing  the  ensembles  for  randomness 
are  presented  in  a  later  part  of  this  chapter. 

The  A/D  conversion  not  only  reduces  the  continuous  time  process 
to  a  discrete  time  process,  but  it  also  quantizes  the  value  of  the 
signal  at  each  time  into  one  of  a  finite  number  of  discrete  values.  A 
12-bit  A/D  converter  was  used  for  the  digitization.  The  instantaneous 
value  of  the  electrical  signal  was  thus  converted  to  an  integer  in  the 
range  of  -20^8  to  +20^7,  giving  a  dynamic  range  of  approximately  66  cLB. 
Thus  the  samples  of  the  random  process  are  no  longer  continuously 
valued,  but  instead  take  on  discrete  values.  This  aspect  of  the  sampling 
has  a  significant  effect  upon  the  statistical  testing  of  the  ensembles. 
Most  of  the  statistical  tests  assume  that  sampling  is  from  a  continuous 
distribution  and  thus  no  two  or  more  samples  will  have  precisely  the 
sane  value.  This  assumption  cannot  always  be  met  when  the  samples  are 
quantized  into  a  finite  number  of  values.  Thus  the  statistical  tests 
must,  be  adjusted  to  account  for  ties  in  the  sample  data. 

The  previous  discussion  has  served  to  point  out  the  assumptions 
inherent  in  the  technique  used  by  this  study  to  form  sample  ensembles. 

The  data  upon  which  the  computations  in  this  study  are  made  are  limited 
ensembles  of  quantized,  fin.it. e  range,  discrete  time,  random  processes. 
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It  is  assumed  that  these  data  adequately  represent  continuously  valued, 
continuous  time  random  processes.  Generally,  the  statistical  tests 

account  for  the  fact  that  the  data  is  a  limited  ensemble  of  samples 
from  a  continuous  distribution  and  also  generally  account  for  quantiza¬ 
tion  effects  which  lead  to  ties  in  the  data.  Thus  finite  ensemble 
size  and  quantization  are  not  significant  limitations.  The  limited 
dynamic  range  of  the  sampled  data  must  be  considered  during  data  collec¬ 
tion  and  digitization  to  ensure  that  virtually  the  entire  distribution  of 
values  of  the  reverberation  process  falls  within  the  available  dynamic 
range  of  the  data  collection  and  digitization  system.  It  must  be 
assumed  that  the  extreme  tails  of  the  distribution  of  the  random  process 
which  were  not  included  in  the  sampled  data  due  to  limited  dynamic  range 
contribute  negligibly  to  the  results  of  this  study.  Observation  of  the 
sampled  data  indicated  that  only  occasionally  was  the  value  of  the 
reverberation  process  outside  the  dynamic  range  of  the  A/D  converter, 
so  that  it  is  felt  that  the  dynamic  range  is  also  not  a  significant 
limitation . 

IV. 2  Statistical  Testing  of  Ensembles 

Since  the  reverberation  is  treated  as  a  random  process  and  the 
reverberation  ensemble  as  a  random  sampling  of  the  parent  population  of 
the  random  process,  then  it  is  assumed  that  the  ensembles  consist  of 
random  variables  which  are  independent  and  identically  distributed.  In 
an  attempt  to  verify  that  the  collection  of  reverberation  events  forms 
a  valid  ensemble,  each  ensemble  was  tested  for  randomness  or  independence. 
Since  the  elements  of  a  valid  ensemble  must  be  not  only  independent  but 
also  identically  distributed,  further  tests  for  homogeneity  were  also 


iS  V/ 

applied  to  each  ensemble.  1  An  ensemble  was  accept. e i  as  a  valid 
sample  ensemble  it'  the  hypotheses  of  randomness  and  homogeneity  could 

not  be  rejected  by  these  tests  at  a  reasonable  level  of  s igni f i canoe . 

After  the  sample  ensembles  were  verified,  they  were  then  tested 
for  normality.  That  is,  each  ensemble  was  tested  to  determine  if  it. 
represented  a  sampling  from  a  normal  parent  population. 

IV.?. 1  General  Remarks  on  Hypothesis  Testing 

One  of  the  requirements  of  hypothesis  testing  is  to 
establish  a  criterion  for  accepting  or  rejecting  an  hypothesis  for  an 
ensemhLe  or  group  of  ensembles.  Before  engaging  in  a  discussion  of  the 
details  of  statistical  hypothesis  testing,  a  few  qualitative  remarks 
may  help  to  clarify  the  criterion  chosen  in  this  study.  The  basic 
"output"  of  an  hypothesis  test  upon  which  this  criterion  is  usually 
based  is  the  probability  of  making  a  mistake  by  rejecting  the  test  hypo¬ 
thesis.  That  is,  an  hypothesis  test,  is  performed  on  an  ensemble  and 
the  basic  result  of  the  test  is  this  probability.  Typically  if  this 
probability  is  small,  that  is,  the  probability  of  erring  by  rejecting 
the  hypothesis  is  small,  then  the  hypothec i can  b<  rejected  for  t.iie 
ensemble  with  some  degree  of  confidence,  while  keeping  in  mind  that,  there 
is  some  small  probability  of  being  wrong. 

Jn  the  present  study  many  ensembles  are  obtained,  one  .-V 
each  sample  time  for  each  channel.  Instead  of  accepting  or  rejecting  Indi¬ 
vidual  ensembles,  as  is  typically  done,  it  is  desirable  to  either  accept 
or  reject  the  entire  set  of  ensembles  for  each  channel.  That  is,  it  is 
desirable  to  either  accept  or  reject  the  hypothesis  being  tested  for 
each  channel  over  the  entire  10  ms  of  the  recorded  reverberation  returns. 
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Of  rourie,  this  acceptance  technique  assumes  that  the  property  of  the 

data  which  is  being  tested,  such  as  normality  or  homogeneity,  does  not 

change  during  the  10  ms  reverberation  return.  This  may  not  always  be  the 
1 5  >'  i 

case,  5  and  thus  it  may  be  necessary  to  accept  one  segment  of  the  data 
while  rejecting  another  segment. 

In  any  case  it  is  necessary  to  establish  a  procedure  which 
will  allow  the  acceptance  or  rejection  of  the  entire  set,  or  some  segment  of 
the  entire  set,  of  -Tf-v-ables  for  each  channel.  One  such  procedure  is  to 
examine  the  percentage  of  ensembles  which  have  a  small  probability  of 
erring.  If  a  disproportionately  large  percentage  of  the  ensembles  have 
small  probabilities  of  erring  by  rejecting  the  hypothesis,  then  it  is 
reasonable  to  reject  the  hypothesis  for  the  entire  set  of  ensembles.  On 
the  other  hand,  if  the  percentage  of  ensembles  with  small  probabilities 
is  commensurate  with  the  level  of  probabilities,  then  the  hypothesis  is 
accepted  for  the  entire  set  of  ensembles.  For  example,  if  approximately 
10  percent  of  the  ensembles  have  a  probability  of  making  a  mistake  less 
than  or  t  qual  to  0.1,  then  the  hypothesis  can  be  accepted.  If,  however, 
the  percentage  of  ensembles  is  significantly  different  from  10  percent, 
then  the  hypothesis  is  rejected  at  this  level  of  0.1.  The  probability 
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at  which  an  hypothesis  is  tested  is  referred  to  as  tin'  level  of  signif¬ 
icance.  .n  this  study,  the  level  of  significance  was  chosen  to  be  0.1 
for  all  the  hypothesis  tests. 

Having  qualitatively  discussed  the  results  of  an  hypo¬ 
thesis  lost,  and  established  a  criterion  for  acceptin'*  op  rejecting  an 
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hypothesis,  a  more  detailed  discussion  m'  l.he  general  proceaures  for 
applying  these  tests  is  appropriate.  F'ach  of  the  statistical  tests 
applied  in  this  study  tests  a  single  null  hypothesis  against  a  single 
alternative  hypothesis  H^ .  No  testing  is  performed  against  multiple 
hypotheses.  For  example,  the  test  for  randomness  tests  the  null  hypoth¬ 
esis 

H  :  the  sample  ensemble  is  random, 
against  the  alternative  hypothesis 

:  the  sample  ensemble  is  not  random. 

To  accept  the  alternative  hypothesis  is  to  reject  the  null  hypothesis 
H^.  No  attempt  is  made  to  discriminate  between  various  types  of  non¬ 
randomness  in  the  test  itself.  However,  the  particular  test  chosen  may 
be  more  sensitive  to  one  class  of  nonrandomness  than  another.  Thus  the 
test  should  be  chosen  which  is  most  sensitive  to  the  type  of  nonrandom¬ 
ness  expected  in  the  data  in  order  to  produce  results  which  can  be  con- 
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fidently  relied  upon. 

Hypothesis  testing  is  performed  by  first  computing  a 

test  statistic  which  is  a  function  of  the  sample  ensemble  being  tested. 

For  example,  denoting  the  random  variables  of  the  sample  ensemble  by 

v.(t  ,o>)  and  the  test  statistic  by  T.  (L  ),  then 
in  in 

T. ( t  )  =  H ( v . ( t  ,w) ) 
i  n  i  n 

where  F  is  some  function  which  defines  the  particular  test,  being  used. 
The  test  statistic  is  itself  tnerefore  a  random  variable  which  has, 
(perhaps  unknown)  probabi 1 i ty  distributions  under  H^  and  H  .  That  is, 
if  the  sample  ensemble  meets,  the  criteria  of  the  null  hypothesis,  then 
the  test  statistic  will,  have  a  particular  probability  distribution.  If 
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the  sample  ensemble  meets  the  criteria  of  the  alternative  hypothesis, 
then  the  test  statistic  will  have  a  (hopefully)  different  probability 
distribution. 

The  probability  distribution  of  the  test  statistic  under 
the  null  hypothesis  is  generally  known  or  can  be  approximated.  From  this 
distribution,  the  probability  of  obtaining  or  exceeding  the  value  of  the 
statistic  which  was  calculated  from  the  sample  ensemble,  assuming  the 
sample  ensemble  met  the  criteria  of  the  null  hypothesis,  can  be  calcu¬ 
lated.  If  this  probability  is  small,  then  the  assumption  of  the  null 
hypothesis  is  rejected  and  the  alternative  hypothesis  is  accepted.  This 
probability  was  referred  to  previously  as  the  probability  of  making  a 
Type  I  error,  and  is  designated  a.  Thus  a  is  the  probability  of  announc¬ 
ing  when  HQ  is  true,  written  as 
a  =  P(H1|Hq) 

The  probability  a  provides  a  measure  of  the  confidence 
with  which  the  alternative  hypothesis  can  be  accepted.  If  a  is  very 
small,  then  there  is  little  chance  of  making  an  error  in  announcing  , 

and  thus  one  can  be  confident  about  announcing  H  .  However,  if  a  is 
not  small,  then  cannot  be  reliably  accepted  and  instead  is 
accepted.  Unfortunately,  the  value  of  a  provides  no  indication  of  the 
confidence  with  which  can  be  accepted.  In  order  to  confidently  accept 
Hq,  the  probability  of  announcing  H^when  is  true,  written  as 
P(Hq|H1),  must  be  known  to  be  small.  This  probability  is  known  as  the 
probability  of  making  a  Type  II  error  and  is  designated  by  8,  i.e., 

8  =  P(H0|H1) 


1-8  is  sometimes  referred  to  as  the  power  of  a  test.  'J'o  calculate  this 
probability,  the  distribution  of  the  test  statistic  under  the  alterna¬ 
tive  hypothesis  must  be  known.  This  distribution  is  generally  not 
known  for  the  types  of  tests  used  in  this  study.  Thus,  the  best  that 
can  be  said  is  that  if  H  cannot  be  accepted,  then  cannot  be  reliably 
rejected. 

In  an  attempt  to  minimize  this  deficiency,  several 
criteria  are  applied  to  hypothesis  tests  in  an  effort  to  define  a  test 
which  will  most  likely  be  the  most  powerful  test  to  use.  One  of  the 
most  common  criteria  is  asymptotic  relative  efficiency  (ARE).  It  is 
a  comparison  of  the  efficiency  of  one  test  to  another  test,  and  is 
defined  as  follows.  Let  be  the  number  of  samples  that  test  A 
requires  to  achieve  a  power  1-8  at  a  given  significance  level  a,  and 
let  N  be  the  number  of  sample j  that  test  B  requires  to  achieve  the 

D 

same  power  at  the  same  significance  level.  Tests  A  and  B  test  the  same 
null  hypothesis  against  the  same  alternative  hypothesis  and  each 
has  test  statistics  which  are  asymptotically  normally  distributed.  The 
ARE  is  the  ratio  of  to  H0  in  the  limit  as  NA  and  Ng  go  > 

infinity  and  as  goes  to  H^. 

The  ARE  is  useful  primarily  as  a  comparison  of  several 

different  tests  to  the  same  reference  test.  While  the  magnitude  of 

the  ARE  has  only  a  limited  meaningfulness,  a  comparison  of  the  ARE  of 

several  comparable  tests  is  generally  a  good  indication  of  the  ranking 

in 

of  the  power  of  the  tests.  Thus  the  test  with  the  largest  ARE  is 
usually  the  most  powerful  test.  If  the  tests  are  nonparametric ,  i.e., 
no  assumption  is  made  as  to  the  form  of  the  probability  distribution 
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of  the  sample  ensemble  (except  possibly  for  continuity  or  some  other 
rather  general  assumptions),  then  usually  a  most  powerful  parametric  test 
is  used  as  the  reference  test  from  which  the  ARE  of  the  nonparametric 
tests  are  determined.  In  this  study,  the  tests  for  randomness,  indepen¬ 
dence,  homogeneity,  and  one  test  for  normality  are  all  nonparametric  tests. 

Alternatively,  an  empirical  method  is  often  used  with 
parametric  tests  to  estimate  the  efficiency  of  various  tests  when  the 
distribution  of  the  test  statistic  under  the  alternative  hypothesis  is 
unknown.  The  sample  ensemble  is  assumed  to  have  a  certain  distribution 
under  the  alternative  hypothesis,  and  from  this  distribution  the  distri¬ 
bution  of  the  test  statistic  under  the  alternative  hypothesis  is  empiri¬ 
cally  calculated  for  various  sample  sizes.  It  is  then  possible  to  compute 
the  power  of  the  test  from  this  empirical  distribution.  This  procedure 
is  usually  repeated  for  a  number  of  different  sample  ensemble  distribu¬ 
tions.  In  this  way  the  power  of  several  tests  can  be  compared  for  a 
number  of  different  distributions  to  determine  the  test  which  generally 
has  the  largest  power  over  the  widest  class  of  distributions.  Alterna¬ 
tively,  if  the  sample  data  is  known  a  priori  to  come  from  one  of  a 
restricted  class  of  distributions,  then  a  test  can  be  chosen  which  has 
the  largest  power  over  this  restricted  class  of  distributions.  In  this 
study,  three  of  the  tests  for  normality  are  parametric  tests. 

The  tests  used  in  this  study  were  generally  selected  to 
provide  the  most  power  for  the  type  of  data  being  tested  and  the  likely 
alternatives.  The  runs  up  and  down  test  was  used  as  a  test  for  random¬ 
ness,  while  the  Kolmogorov-Smirnov  two-sample  test  and  the  Wilcoxon  rank- 
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sum  test  wer«  us.-d  tor  t  <  sting  hr-nogeneit .  Pin-  four  1  "S*  s  f.  r  norma  i  i  * ;/ 
were  Pearson  *  s  test  of  skewness,  Pearson  ’  s  test  of  kurtosis,  e '  Agost  i  no  ’  s 
test,  and  the  Kolmogorov-Smi rnov  one-sample  test.  The  (two-tailed)  level 
of  significance  was  chosen  to  be  0.1  for  ail  tests. 

The  tests  were  performed  on  each  of  the  l.-'OO  ensembles 
from  each  of  the  13  channels.  For  each  test,  the  probability  of  a  Type  I 
error  was  computed  and  compared  to  the  level  of  significance.  The  per¬ 
centage  of  the  3200  ensembles  from  each  channel  that  failed  a  particula" 
test,  i.e.,  that  had  a  probability  of  a  Type  I  error  telow  the  level  of 
significance,  was  also  computed.  In  order  to  determine  if  the  hypothesis 
being  tested  should  be  accepted  for  all  3200  ensembles,  the  results  of  a 
particular  test  were  treated  ns  outcomes  from  independent,  repeated 
Bernoulli  trials.1^  That  is,  in  this  scheme  two  outcomes  from  a  test  arc 
possible:  the  probability  of  a  Type  T  error  n  is  greater  than  or  equal 

to  the  significance  level,  or  it.  in  less  than  tic.-  significance  level.  The 
probability  tha*  a  is  less  than  the  signi  fi  e.me<-  level  under  the  null 
hypothesis  is  .lust,  the  value  of  the  signi  ficance  1  eve-1  .  That  is,  for  a 


si  gni  finance  level  of  0.1  t.  tie  re  is  a  probability  of  o.l  that  a  will  be 
less  than  0.1  if  the  null  hypothesis  is  true.  If  tins  probability  is 
denoted  by  p,  then  for  all  the  tests  considered  in  this  study, 

V  =  l'robf  u<0 . 1  1  =  0. 1 

According  to  the  binomial  law,  the  probability  that  ,  <>ut.  ,  •'  n  inlepenb  nl 
outcomes  each  with  a  probability  j.  of  having  a  value  of  -  less  firm  ‘In 
significance  level,  k  or  more  of  the  outcomes  will  actually  have  values 
of  a  less  than  the  significance  level  is.  given  by 
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B(k;n,p)  =  £  (i  )  P^l-P)0-1 
i=k  '  ' 

This  binomial  distribution  has  a  mean  of  np  and  a  standard  deviation  of 
[np(l-p)]  .  For  a  particular  test  the  number  k  of  outcomes  which  had 

values  of  a<0.1  were  determined  and  the  resulting  binomial  probability 
was  computed  for  each  channel.  If  this  probability  was  small,  less  than 
0.05,  indicating  that  too  many  outcomes  were  below  the  significance  level, 
or  if  it  was  large,  greater  than  .95,  indicating  too  few  outcomes,  then 
the  null  hypothesis  was  rejected  in  favor  of  the  alternative  hypothesis 
for  that  channel.  If  the  probability  was  greater  than  0.05,  then  the 
null  hypothesis  could  not  be  confidently  rejected. 

This  technique  assumes  that  the  outcomes  are  independent. 

It  will  be  shown  in  the  next  chapter  that  the  reverberation  used  for  this 
study  has  a  correlation  time  of  approximately  250  us.  Thus  tests  per¬ 
formed  on  ensembles  representing  sample  times  which  are  less  than  250  us 
apart  may  not  have  independent  outcomes,  although  it  is  difficult  to  deter¬ 
mine  the  precise  relationship  between  the  correlation  of  the  reverbera¬ 
tion  process  and  the  dependence  of  the  outcomes  of  a  test.  The  approach 
taken  in  this  study  to  attempt  to  extract  independent  test  outcomes  was 
to  select  the  outcome  of  the  test  from  every  80th  ensemble,  corresponding 
to  a  time  separation  of  250  us,  from  the  set  of  3200  ensembles  from  each 
channel.  These  1+0  values  of  a  were  then  examined  to  determine  the  num¬ 
ber  k  of  values  less  than  0.1,  and  this  number  was  used  to  compute  the 
binomial  probability. ^  Thus  n=Uo  in  the  computations  of  the  binomial 
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In  the  present  study,  the  reverberation  samples 
cannot  be  so  easily  divided  into  two  distinct  classes  since  the  samples 
are  continuously  valued  and  thus  take  on  a  range  of  values.  One 

possibility  is  to  test  the  sign  of  the  difference  between  the  samples 

hi 

and  the  sample  mean.  However,  if  the  distribution  of  the  samples 
is  symmetric  about  the  mean,  which  is  very  likely  the  case  for  rever¬ 
beration,  then  this  test  will  be  insensitive  to  any  changes  in  the 
variance  of  the  distribution  which  might  have  occurred  during  the 
collection  of  the  sample  ensemble.39  That  is,  the  sample  ensemble  may 
not  be  random  due  to  a  lack  of  homogeneity,  but  the  test  will  fail  to 
detect  it.  A  test  of  this  sort  is  a  test  only  of  the  randomness  of 
the  signs  and  not  a  test  of  the  randomness  of  the  samples. 

An  alternative  pair  of  events  which  is  more 
suitable  for  an  ensemble  of  continuously  valued  samples  is  an  unbroken 
sequence  of  increasing  or  decreasing  observations.  Although  the  proba¬ 
bility  distribution  associated  with  the  total  number  of  runs  of  this 
type  is  different  from  the  probability  distribution  associated  with 
the  previously  described  runs,  it  can  still  be  calculated.  The  total 
number  of  runs  r  in  a  sample  ensemble  of  size  n  is  asymptotically 
normally  distributed  with  mean  and  variance  given  by:^*9 


2n-l 


Thus  the  statistic  r  is  transformed  to  a  normal  deviate 

r-u 

a 


■*  ■>■*«? 
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and  is  referenced  to  the  normal  distribution  function.  Since  the  exact 

distribution  of  r  is  discrete  and  the  normal  approximation  is  continuous, 

a  continuity  correction  is  applied  by  reducing  the  absolute  value  of 

r-n  by  The  test  statistic  is  then: 

z  =  sgn(r-p) • [ |r-n | — ■ 5 ] 
a 

where  sgn(*)  is  the  sign  function. 

To  guard  against  nonrandomness  due  to  either 
too  few  runs  (Z  significantly  less  than  zero)  or  too  many  runs  (Z  sig¬ 
nificantly  greater  than  zero),  the  absolute  value  of  Z  is  used  as  the 
test  statistic.  The  test  is  then  referred  to  as  a  two-tailed  test. 

The  distribution  of  |z|  under  the  null  hypothesis  can  be  easily  deter¬ 
mined  from  the  distribution  of  Z,  and  the  significance  level  of  )z|  can 

h2 

be  calculated  from  the  significance  level  of  Z.  Denoting  the 
distribution  function  of  Z  by  F,(z)  and  the  distribution  function  of 
|z|  by  F|z|(z),  then  since  Z  has  a  symmetric  density  function  about 
zero , 

F|z|(z)  =  2Fz(z)  -  1  ,  z  >  0 

Since 

a  =  P(H1|Hq)  =  P( | Z | > z )  =  1  -  F| z | ( z ) 

then 


a  =  2(1-Fz(z))  =  2P(Z>z)  ,  z  >  0 

Thus  the  calculated  value  z  of  the  test  statistic  jz|  can  be  related  to 
the  normal  distribution  function  and  the  resulting  probability  of 
exceeding  that  value  can  be  multiplied  by  2  to  determine  the  signifi¬ 
cance  level  for  | Z | . 


The  runs  up  and  down  test  assumes  that  there 


a re  no  tied  observations.  If  ties  do  occur,  only  adjacent  tied  values 
could  possibly  affect  the  number  of  runs.  In  order  to  account  for  ties, 
the  ties  are  handled  so  as  to  produce  the  minimum  possible  number  of 
runs  and  then  the  maximum  possible  number  of  runs.  The  probabilities 
associated  with  the  minimum  and  maximum  number  of  runs  are  computed 
to  provide  upper  and  lower  bounds  on  the  level  of  significance. 

The  ARE  of  this  test  was  found  to  be  zero 

1*0 

with  respect  to  several  distribution- free  and  parametric  tests. 

However,  even  though  its  efficiency  is  low,  the  power  will  be  large  due 
to  the  large  number  of  samples  employed  in  this  study.  The  test  was 
chosen  for  use  because  of  its  simplicity  and  economy  as  compared  to  a 
more  efficient  test  such  as  the  Kendall  rank  correlation  test  for  ran¬ 
domness. 

The  results  from  the  runs  up  and  down  test  are 
displayed  graphically  in  Figs.  IV-1  through  IV-i*.  The  plots  show  the 
probability  of  a  Type  I  error  for  each  of  the  3200  ensembles  of  each 
channel.  The  percentage  of  probabilities  below  the  value  of  0.1  are 
computed  for  each  channel  and  indicated  on  the  plots.  The  percentages 
are  based  on  all  3200  samples  and  not  on  an  independent  sampling  of  the 
ensembles.  The  two  values  represent  probabilities  associated  with  the 
minimum  and  maximum  number  of  runs  and  can  be  considered  as  bounds  on  the 
true  percentage.  It  can  be  seen  from  these  plots  that  there  was  not  sig¬ 
nificantly  more  than  10%  of  these  probabilities  below  0.1.  It  will  be 
shown  later  that  probabilities  based  on  independent  Bernoulli  trials  were 


not  small  enough  to  reject  the  random  hypothesis.  Thus  the  entire  set  of 
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FIGURE  IV-1 

THE  RUNS  UP  AND  DOWN  TEST  FOR  RANDOMNESS,  CHANNELS  1-3 
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FIGURE  IV-2 

THE  RUNS  UP  AND  DOWN  TEST  FOR  RANDOMNESS,  CHANNELS  4-6 
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FIGURE  IV-3 

THE  RUNS  UP  AND  DOWN  TEST  FOR  RANDOMNESS,  CHANNELS  7-9 
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FIGURE  IV-4 

THE  RUNS  UP  AND  DOWN  TEST  FOR  RANDOMNESS,  CHANNELS  11-14 
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ensembles  for  all  channels  was  accepted  as  random,  that  is,  randomness 
could  not  be  reliably  rejected. 

The  tests  for  randomness  were  performed  only 
on  the  last  500  samples  from  each  ensemble,  instead  of  the  entire  1081* 
samples.  The  results  of  the  next  section  will  show  that  the  entire 
108U  samples  were  not  homogeneous,  but  that  the  last  500  samples  of 
each  ensemble  were  more  nearly  homogeneous.  Thus  all  of  the  statistical 
analysis  was  performed  only  on  the  last  500  samples  of  each  ensemble. 

IV. 2. 3  Testing  for  Homogeneity 

The  non-par ametric  tests  for  homogeneity  employed  in 
this  study  divide  the  sample  ensemble  into  two  sub-ensembles  of  length 
n  and  m  and  test  that  the  two  sub-ensembles  are  identically  distributed. 
Specifically ,  they  test  the  null  hypothesis: 

Hq :  each  of  the  possible  combinations  of  sub-ensembles,  one 
of  length  n  and  the  other  of  length  m,  obtainable  from  the 
total  ensemble  of  length  n  +  m  was  equally  likely  to  have 
become  the  one  which  was  actually  obtained; 
against  the  alternative  hypothesis : 

H^:  certain  combinations  of  sub-ensembles  are  more  likely  than 
others . 

For  the  null  hypothesis  to  be  true,  the  two  sub-ensembles  must  be 
identically  distributed.  Since  the  most  likely  source  of  inhomogeneity 
in  the  reverberation  data  is  a  change  in  reverberation  intensity 
between  the  beginning  and  the  end  of  the  recording  of  the  data  (approx¬ 
imately  a  lU  minute  period  of  time)  the  two  sub-ensembles  initially 
consisted  of  the  first  500  samples  and  the  second  500  samples  of  the 
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ensemble.  Two  tests  for  homogeneity  were  employed.  The  first  one 

to  be  described  is  the  Kolmogorov-Smirnov  two-sample  test.  The  second 

is  the  Wilcoxon  rank-sum  test. 

IV. 2. 3.1  The  Kolmogorov-Smirnov  Two-Sample  Test 

The  test  statistic  for  the  Kolmogorov-Smirnov 

two-sample  (KS-2)  test  is  the  maximum  absolute  difference  between  the 

empirical  cumulative  distribution  functions  of  the  two  sub-ensembles, 

denoted  D  .  Although  the  exact  distribution  of  D  can  be  calculated, 
ran  mn 

1|1| 

for  large  sample  sizes  the  Smirnov  approximation  is  more  useful.  If 

D  ^  is  transformed  to  the  statistic  Z,  where 
mn 


=  D 


mn  m+n 

then  Z  has  as  its  limiting  distribution  the  Smirnov  distribution. 

kh 

a  correction  for  continuity,  Z  becomes 


With 


Z  =  D 


V^  + 


mn  m+n 


where 


2 

C  =  for  n  =  km,  k  =  1,2,3 

2 


=  5*'rrT  for  n  #  km 

The  Smirnov  distribution  F(z)  was  approximated  as 
C  0  ,  -“><z<.22 


F(z)  =  < 


.22<_z_<  .80 


1  -  2[exp(-2z2)-exp(-8z2)+exp(-l8z2)], 
J.  ,  z>3.15 

The  probability  of  a  Type  I  error  is 


.SO<z<3.15 


a  =  1  -  F(z) 
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The  probability  of  a  Type  I  error  is  plotted 
in  Figure  IV-5  for  all  3200  ensembles  of  channel  11  using  the  first 
1000  samples  of  each  ensemble.  Slightly  over  30#  of  the  ensembles  had 
a  probability  of  Type  I  error  under  0.1,  indicating  a  lack  of  homo¬ 
geneity  between  the  first  500  and  second  500  samples  of  these  ensembles. 
Figure  IV-6  gives  a  strong  indication  of  the  source  of  the  inhomogeneity. 
100-sample  estimates  of  the  variance  were  computed  throughout  each 
ensemble  of  1081+  samples  and  are  plotted  for  every  32nd  ensemble  of  the 
entire  set  of  3200  ensembles  from  channel  11.  Since  the  variance  is 
proportional  to  the  intensity  of  the  reverberation,  this  plot  is  indic¬ 
ative  of  the  change  in  reverberation  intensity  during  the  approximately 
lU  minutes  necessary  to  record  the  1084  reverberation  returns  and  the 
10  ms  during  which  the  returns  were  digitized.  It  is  obvious  that 
the  reverberation  intensity  was  significantly  lower  during  the  first 
part  of  the  data  recording  period,  after  which  it  increased  to  a  maximum, 
then  decreased  somewhat  and  was  relatively  stable  during  the  last 
part  of  the  lU  minute  recording  period.  The  change  of  the  reverbera¬ 
tion  intensity  for  a  single  ensemble  is  shown  in  Fig.  IV-7  and  further 
demonstrates  the  change  in  reverberation  intensity  which  is  character¬ 
istic  of  all  the  ensembles.  Thus  it  was  felt  that  this  change  in 
intensity  was  the  cause  of  the  lack  of  homogeneity  in  the  data  set. 

Since  the  intensity  over  the  last  500  samples  of  each  ensemble  appears 
to  be  more  stable,  the  last  500  samples  should  be  more  nearly  homo¬ 
geneous.  For  this  reason  it  was  decided  to  restrict  all  further 
analysis  of  the  reverberation  returns  to  the  last  500  samples  of  each 


ensemble . 
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FIGURE  IV-5 

THE  KOLMOGOROV-SMIRNOV  TWO-SAMPLE  TEST  FOR  HOMOGENEITY,  CHANNEL  11 

1000  SAMPLES  PER  ENSEMBLE 
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Figures  IV-8  through  TV-11  display  the  results 
of  the  KD-2  test  for  all  13  channels,  utilizing  only  the  last  500 
samples  of  each  ensemble.  The  percentage  of  samples  below  the  threshold 
of  0.1  varies  generally  between  10  and  20  percent.  These  percentages 
indicate  that  there  is  still  some  inhomogeneity  even  in  the  last  100 
samples,  bur  it  is  not  as  severe  as  in  the  first  1000  samples.  Even 
though  this  inhomogeneity  is  present.,  it.  was  felt  that  it  woula  not 
be  a  severe  limitation  to  the  rest  of  the  analysis. 

It  can  be  observed  that  physically  adjacent 
channels  in  the  vertical  array  fail  the  KS-2  test  in  approximately  the 
same  regions.  For  example,  channels  1-3  (Fig.  IV-8)  all  fail  signifi¬ 
cantly  in  the  region  between  ensembles  2200-2300.  However,  there  are 
some  exceptions  also,  like  the  region  around  ensemble  2000,  where  channel 
3  fails  but  channels  1  and  2  do  not.  In  general,  this  similarity  in 
Ihe  failure  regions  indicates  a  high  degree  of  correlation  between 
these  channels.  Examination  of  the  four  horizontal  channels  (Fig.  IV-11 ) 
shows  little  similarity  between  adjacent  channels,  which  would  indicate 
a  lower  level  of  correlation. 

TV. 2. 3.2  The  Wilcoxon  Rank-Sum  Test 

The  rank  of  a  sample  within  an  ensemble  repre¬ 
sents  its  relative  size  within  the  ensemble.  The  smallest  sample  has 
a  rank  of  1,  the  next  smallest  a  rank  of  2,  and  so  forth.  If  the 
ensemble  is  divided  into  two  sub-ensembles  of  length  n  and  m  (n<m) , 
then  the  Wilcoxon  rank-sum  statistic  W  is  the  sum  of  the  ranks  of  the 
samples  from  the  sub-ensemble  of  length  n  among  the  entire  n+m  samples. 
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FIGURE  IV-10 

THE  KOLMOGOROV-SMIRNOV  TWO-SAMPLE  TEST  FOR  HOMOGENEITY',  CHANNELS  7  9 
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FIGURE  IV  11 

KOLMOGOROV-SMIRNOV  TWO-SAMPLE  TEST  FOR  HOMOGENEITY,  CHANNELS  1 1  14 
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The  exact  null  distribution  of  W  is  known  and  has  a  mean 

W  =  %n(n+m+l) 

? 

and  a  variance  S  : 

w 

„2  _  nm(n+m+l) 

Sw  “  12 - 

and  is  symmetric  about  the  mean.140  The  distribution  is  asymptotically 
normal  so  that  for  large  sample  sizes  the  normal  approximation  with 

continuity  correction  can  be  used,  where  the  N(0,l)  deviate  is 

Z  -  Bgn(W-W)«[|W-W|-.5l 

S 

w 

The  two-tailed  version  of  the  test  was  used,  so  |z|  was  referred  to  the 
normal  probability  distribution.  A  correction  for  the  normal  approxi- 

1*5 

mation  was  also  applied.  The  calculated  value  z  of  the  test  statis¬ 
tic  | Z |  was  referenced  to  the  normal  probability  distribution  to  deter¬ 
mine  the  probability  P^  of  exceeding  that  value.  was  then  reduced 

by  the  factor 
2  2 

n  +m  +nm+n+m  ,  , 

20nm( n+m+1 )  g  2 

where  g(z)  =  27T-Js(z3-3z)exp(-z2/2) 

This  reduced  probability  was  multiplied  by  2  for  the  two-tailed  level 
of  significance. 

The  Wilcoxon  rank-sum  test  assumes  that  there 
are  no  tied  observations.  To  account  for  ties,  the  test  statistic  W 
is  computed  by  handling  ties  so  that  W  is  as  small  as  possible,  and 
then  as  large  as  possible.  The  associated  probabilities  are  assumed 
to  provide  a  bound  on  the  true  probability. 
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The  test  is  particularly  sensitive  to 
inhomogeneities  due  to  unequal  locations,  but  will  also  detect  other 
differences.  For  testing  for  unequal  locations  between  populations  which 
have  identical  shapes  and  variances,  the  Wilcoxon  test  has  an  ARE  which 
is  no  lower  than  .  86U  relative  to  the  two-sample  t  test.^  It  is 
more  powerful  than  the  KS-2  test  for  detecting  changes  in  location  of 
normal  populations.  In  general,  the  Wilcoxon  test  is  more  powerful  or 
almost  as  powerful  as  other  non-parametric  tests  for  detecting  changes 
in  location. 

The  upper  and  lower  limits  of  the  two-tailed 
probability  of  a  Type  I  error  are  plotted  in  Figs.  IV-12  through  IV-15. 
Although  a  few  channels  have  significantly  more  than  10$  failures  (as 
much  as  about  15$),  most  channels  have  less  than  10$  failures.  This 
indicates  that  there  is  not  a  significant  inhomogeneity  due  to  a  change 
in  location.  It  also  indicates  that  the  Wilcoxon  test  is  not  as  sensi¬ 
tive  as  the  KS-2  test  to  the  type  of  inhomogeneities  which  are  present. 

A  comparison  of  the  results  from  the  Wilcoxon  and  KS-2  tests  shows 
that  in  virtually  every  instance  in  which  the  Wilcoxon  test  failed,  the 
KS-2  test  also  failed.  However,  the  opposite  was  not  true.  Thus  the 
Wilcoxon  test  did  not  detect  inhomogeneities  which  were  not  already 
detected  by  the  KS-2  test. 

IV. 2.1+  Summary  of  Results  of  Tests  for  Randomness  and  Homogeneity 
The  purpose  of  the  tests  for  randomness  and  homogeneity 
was  to  determine  if  the  collection  of  ensembles  for  each  channel  could  be 
considered  to  contain  independent,  identically  distributed  random  variables. 
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FIGURE  IV-12 

THE  WILCOXON  RANK-SUM  TEST  FOR  HOMOGENEITY,  CHANNELS  1-3 
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Three  statistical  tests  were  applied,  one  test  for  randomness  and 
two  tests  for  homogeneity.  As  was  described  earlier,  the  probability  of 
a  lype  I  error  for  every  80th  ensemble  was  compared  to  the  level  of  sig¬ 
nificance  of  0.1  and  the  total  number  of  probabilities  below  this  level 
were  tabulated  for  each  channel.  The  probability  that  this  number  or  more 
would  be  below  the  level  of  significance  was  computed  for  each  channel  and 
each  of  the  three  tests.  The  results  are  reproduced  in  Table  IV. 1. 

Although  there  are  some  probabilities  below  .05  or  greater  than  .95  for 
one  test  of  a  particular  channel,  this  condition  did  not  occur  for  the  other 
tests  for  that  channel.  For  example,  the  binomial  probability  resulting 
from  the  runs  test  for  channel  5  has  a  value  of  .99,  but  the  KS-2  test 
and  Wilcoxon  tests  do  not  exhibit  such  extreme  probabilities.  Thus  con¬ 
sidering  the  three  tests  together,  the  assumption  of  a  statistically  valid 
ensemble,  i.e.,  consisting  of  i.i.d.  random  variables,  cannot  be  confi¬ 
dently  rejected  for  this  channel  or  any  of  the  channels.  Therefore  the 
collection  of  reverberation  ensembles  was  accepted  as  valid  ensembles  for 
all  channels. 

IV. 2. 5  Testing  for  Normality 

Each  ensemble  of  each  channel  was  further  tested  for 
univariate  normality  using  four  different  tests:  Pearson's  test  of 
skewness,  Pearson's  test  of  kurtosis,  D'Agostino's  test,  and  the  Kolno- 
gorov-Smirnov  one-sample  test.  Each  one  tests  the  null  hypothesis: 

H  :  the  samples  of  the  ensemble  consist  of  random  variables 


drawn  from  a  normal  parent  population; 


TABLE  TV. 1 

BINOMIAL  PROBABILITY  OF  k  OR  MORE  OUTCOMES  BELOV.' 
THE  LEVEL  OF  SIGNIFICANCE  FOR  THE  STATISTICAL 
TESTS  FOR  RANPOMNr'SC  ALT'  HOMOGENEITY 


Channel  Number 

Statistical  Test 

Runs 

KS-2 

Wilcoxon 

I 

.',8-.  3 

.10 

. in-. i0 

2 

.'f8-.  7  8 

.01 

.10-. 10 

3 

.78-. 78 

.21 

.58-. 73 

1* 

.21-. 21 

.04 

.78-. 73 

5 

.  99- . 99 

.10 

.98-. 7 3 

6 

.73-. 92 

.37 

. 21  - . 37 

T 

O 

-Cr* 

1 

H- < 

O 

.10 

. o4-.nl 

8 

.98-. 93 

.02 

.04-. 10 

9 

•21-. 37 

.10 

.21-.. ML 

11 

.10-. 10 

.  92 

.0O_. qq 

12 

. 37- • 37 

.78 

.?3-. op 

13 

.37-. 37 

.04 

.21 -.21 

14 

.10-. 10 

.21 

.98-. 98 
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against  the  alternative  hypothesis: 

H^:  the  samples  are  not  normal. 

No  a  priori  values  were  assigned  to  the  mean  and  variance  of  the  ensem¬ 
bles;  that  is,  the  tests  for  normality  were  not  testing  to  determine 
if  the  ensembles  were  normally  distributed  with  a  particular  mean  and 
variance.  Rather,  the  mean  and  variance  were  assumed  to  be  the  sample 
mean  and  variance  computed  from  each  ensemble.  Designating  the  ith 

sample  from  an  ensemble  of  size  N  by  X^ ,  then  the  sample  mean  m  and 

2 

the  sample  variance  s  are  given  by 


m 


i  N 

=  I  y 

N 

1=1 


X. 

l 


2 

s 


1 

N-l 


N 

i=l 


IV. 2. 5.1  Pearson's  Test  of  Skewness 

The  test  statistic  for  Pearson's  test  of 
skewness  is  Just  the  sample  skew  /b^"  ,  given  by 

mo 


m_ 


1  N 

where  =  jj  £  (X^m) 

i=l 


N-l  2 
m2  N  s 

h6 

The  skew  is  asymptotically  normally  distributed  with  a  mean  of  zero 


and  a  variance  of 


hi 


6( N-2) 

(N+l ) (N+3) 


Thus  for  large  sample  sizes,  the  sample  skew  can  be  transformed  to 
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N(0,l)  deviate 


J  6( M-g)  * 
f  (N+l)(N+3) 

and  Z  can  be  referenced  to  the  normal  distribution  function  to  determine 

the  probability  of  a  Type  I  error.  For  a  two-tailed  test,  the  test 

statistic  is  |z|  and  the  resulting  probability  is  multipled  by  2,  as 

usual.  An  empirical  study  of  this  and  other  tests  for  normality  over 

a  wide  class  of  alternative  distributions  has  shown  that  it  is  one  of 

the  more  powerful  tests  for  normality  when  used  in  conjunction  with 

18 

Pearson's  test  of  kurtosis,  described  in  the  next  section. 

Results  from  this  test  are  displayed  in  Figs. 
IV-16  through  IV-19.  The  percentage  of  failures  for  each  channel  ranged 
from  approximately  17%  to  29%,  indi-  iting  some  departure  from  normality. 
There  does  not  appear  to  be  any  correlation  between  the  failure  regions 
of  the  skewness  test  and  the  failure  regions  of  the  KS-2  test,  which 
is  some  small  indication  that  the  level  of  inhomogeneitv  present  in 
the  data  was  not  sufficient  to  upset  the  skewness  test. 

IV. 2. 5.2  Pearson's  Test  of  Kurtosis 

The  test  statistic  for  this  test  is  just  the 
sample  kurtosis  b0,  given  by 

b2  “  2 
m2 

w  U 

where  *  jj  ]£  (X^m) 

i=l 


The  distribution  of  the  kurtosis  only  very  siowly  approaches  the  normal 
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FIGURE  IV- 16 

PEARSON'S  TEST  OF  SKEWNESS  FOR  NORMALITY,  CHANNELS  1-3 
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FIGURE  IV-17 

PEARSON'S  TEST  OF  SKEWNESS  FOR  NORMALITY,  CHANNELS  4  6 
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FIGURE  IV-18 

PEARSON'S  TEST  OF  SKEWNESS  FOR  NORMALITY,  CHANNELS  7-9 
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distribution,  which  makes  the  normal  distribution  unsuitable  as  an 
approximate  distribution.  Rather,  the  distribution  is  usually  approx- 

k6 

imated  by  the  Pearson  Type  IV  curve,  although  other  approximations 
1*9 

have  been  used.  Because  of  the  difficulty  of  calculating  these 
distributions,  the  probability  of  a  Type  I  error  was  not  calculated  for 
this  statistic.  Instead,  Figs.  IV-20  through  IV-23  are  plots  of  the 
statistic  b^  itself  with  an  indication  of  the  upper  and  lower  bounds  of 
b^.  Values  of  b^  outside  these  bounds  result  in  a  two-tailed  probabil¬ 
ity  less  than  0.1.  These  upper  and  lower  bounds  were  derived  by 

U6 

Pearson  using  the  Type  IV  curve.  The  percentage  of  points  outside 
these  bounds  is  indicated  on  each  plot. 

The  percentage  of  failures  varied  from  approx¬ 
imately  2h%  to  almost  60% ,  indicating  a  significant  deviation  of  the 


kurtosis  from  its  value  under  the  normal  distribution.  The  kurtosis 
is  generally  larger  than  expected,  which  implies  that  the  tails  of  the 
sample  distribution  are  larger  than  the  normal  tails.  These  larger 
tails  would  be  expected  if  there  were  a  greater  number  of  large  ampli¬ 
tude  returns.  Taken  together,  the  tests  of  skewness  and  kurtosis  indi¬ 
cate  that  the  reverberation  returns  were  significantly  non-normal,  with 
most  of  the  difference  due  to  a  larger  kurtosis. 

IV. 2. 5- 3  D'Agostino’s  Test  for  Normality 

The  test  statistic  D  for  this  test  is 


where  T 
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FIGURE  IV- 22 

PEARSON’S  TEST  OF  KURTOSIS  FOR  NORMALITY,  CHANNELS  7-9 
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FIGURE  IV-23 

PEARSON’S  TEST  OF  KURTOSIS  FOR  NORMALITY,  CHANNELS  11-14 
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The  expected  value  of  1)  and  its  asymptotic  standard  deviation  are 
approximately  given  hy 


E(D) 


_i _  l  1 

2/7  »  V  308-1)" 


1Z8( H-l ) 


and 


asd(D)  - 


0.029986 

v^r 


A1 1, hough  the  distribution  of  P  is  asympt ot- i rally  normal,  it  is  more 

accurate  to  transform  D  to  a  normal  deviate  using  'ronish-Fishor  •>;pan 

sions.  The  statistic  is  first  standardized  by  the  transformation 
v  =  D-E(P) 

asd  f  P*> 

Then  a  Cornish-Fisher  expansion  using  the  first  four  moments  of  P  is 
used  to  transform  to  the  N(0,l)  deviate  : 


where 


Z  =  Y  -  [iTl(Y2-l)  +  ^Y?(Y3-3Y)-~ y^(lY3-7Y)] 

v  ,  ri-Jiil 

1  /r 


y  =  1QT.9 
2  N 


A  two-tailed  test  is  performed,  so  |z|  is  referenced  to  the  normal 
distribution  and  the  resulting  probability  is  multiplied  by  2  to 
determine  the  probability  of  a  Type  I  error.  An  empirical  power  study 
has  been  performed  over  a  wide  class  of  alternative  distributions  and 
indicates  that  this  test  has  good  power  over  most  of  the  distribu¬ 
tions  . 


The  probability  of  a  Type  I  error  is  plotted 
in  Figs.  IV- through  IV-27  for  each  of  the  channels.  A  high  degree 
of  correlation  can  be  observed  between  these  results  and  the  results 


from  the  kurtosis  test.  Both  tests  have  approximately  the  same 
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D'AGOSTINO  S  TEST  FOR  NORMALITY,  CHANNELS  1 1  14 


108 


percentage  of  failures  for  each  channel,  and  the  failures  occur  in  the 
same  regions.  Thus  it  appears  that  D'Agostino's  test  and  the  kurtosis 
test  have  approximately  the  same  power  for  the  type  of  distributions 
encountered  in  these  reverberation  data,  and  both  indicate  a  signifi¬ 
cant  departure  of  the  data  from  normality. 

IV.Z'.5.,;  The  Kolmogorov-Smirnov  One-Sample  Test 

The  test  statistic  for  the  Kolmogorov- 
Smirnov  one-sample  (KS-l)  test  is  the  maximum  absolute  difference 
between  the  empirical  cumulative  distribution  function  of  the  data 
and  the  theoretical  distribution  function.  The  distribution  of 

depends  on  N,  the  number  of  samples,  but  is  independent  of  the  theoreti- 

52 

cal  distribution  function.  As  a  test  for  normality,  the  theoretical 

distribution  is  the  normal  distribution  with  a  mean  and  variance  given 

by  the  sample  mean  and  sample  variance  of  the  ensemble  being  tested. 

5°  53 

The  exact  distribution  of  can  be  calculated,  but  for  large 

sample  sizes  the  Smirnov  approximation  can  be  used.  If  is  transformed 
to  the  statistic  Z,  where 

z  =  /r  dn  . 

then  Z  has  as  its  limiting  distribution  the  Smirnov  distribution.  The 
approximation  to  the  Smirnov  distribution  has  been  discussed  in  an 
earlier  section  of  this  chapter. 

As  a  goodness-of-fit  test,  the  KS-l  test  is 

lo 

generally  superior  to  the  Chi-square  test.  However,  as  a  test  for 

normality,  it  generally  has  less  power  than  the  skewness  and  kurtosis  tests 

1,8 

or  D'Agostino's  test. 


The  results  from  the  KS-l  test,  presented  in 


iU  , 

Figs.  IV-28  through  1V-31,  tend  to  confirm  it..;  poor  performance  ns.  a  test 
for  normality  against  the  type  of  distributions  eneountr-rel  in  the.'  nut  a. 

Much  of  tile  poor  performance  of  the  KS-1  test 
can  be  attributed  to  the  use  of  the  sample  mean  and  variance  instead  of 
the  true  (but  unknown)  mean  and  variance.  The  distribution  .  f  tin-  t-.r 
statistic  for  the  KS-1  test  is  derived  on  the  assumption  that  all  param¬ 
eters  of  the  theoretical  distribution  function  are  known.  Wh>  s  ar  .  rs 
are  unknown  and  must  be  estimated  from  the  data,  the  distribution  d'  the- 
test  statistic  is  generally  different,  and  in  fact  depends  on  the  theo¬ 
retical  distribution  function.  Thus  the  probability  of  a  Type  I  error 
computed  from  the  test  statistic  distribution  derived  under  the  condition 
that  parameters  are  known  will  generally  be  different  from  the  true  prob¬ 
ability  of  a  Type  I  error  when  parameters  are  estimated.  Usually  the 
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test  will  be  more  conservative;  that  is,  the  probability  of  a  Typo  1 
error  computed  in  this  way  will  be  larger  than  its  true  value,  resulting 
in  fewer  rejections  of  the  null  hypothesis  at  a  given  significance  level 
than  if  the  probability  of  a  Type  I  error  had  been  computed  from  the  true 
distribution  function.  The  conservative  nature  of  the  test  when  param¬ 
eters  are  estimated  is  demonstrated  in  these  results.  As  can  be  seen, 
every  channel  had  less  than  1 %  failures ,  compared  to  typically  50/?  fail¬ 
ures  from  the  kurtosis  test  of  D'Agostino's  test.  At  a  level  of  sig¬ 
nificance  of  0.1,  approximately  10$  failures  would  be  expectel  even  if 
the  data  were  normal.  Although  efforts  have  been  made  to  compute  the 

rj  ^  o 

distribution  function  of  the  KS-1  test  when  parameters  are  estimated,’1 
the  results  are  extremely  difficult,  to  apply  and  were  not  a,  tempted  in 


this  study. 
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IV. 2. 6  Summary  of  Results  of  Tests  for  Normalit 


Four  different  tests  for  normality  were  applied  to  all 
13  channels  of  data.  The  probability  of  a  Type  I  error  from  every  Both 
ensemble  was  compared  to  the  level  of  significance  of  0.1  and  the  result¬ 
ing  binomial  probability  was  computed.  On  all  channels,  the  results  from 
the  D-test  and  the  kurtosis  test  indicated  a  near-zero  probability  that 
the  collection  of  ensembles  from  each  channel  was  normally  distributed. 

The  skew  test  resulted  in  binomial  probabilities  less  than  .05  for  all 
channels  except  6,  9,  11,  12,  and  13.  For  these  channels  the  probabilitie 
were  .1,  .1,  .21,  .37.  and  .58,  respectively.  Since  the  D-test  and  the 
kurtosis  test  rejected  these  channels,  it  was  felt  that  the  hypothesis  of 
normality  could  be  confidently  rejected  for  all  13  channels.  The  KS-1 
test  was  not  applied  to  the  determination  of  normality  due  to  the 
extremely  conservative  nature  of  the  test  as  it  was  implemented  with  the 
sample  mean  and  variance. 

IV. 3  Summary  of  the  Statistical  Testing 

All  thirteen  channels  of  the  reverberation  data  have  been 
tested  for  randomness,  homogeneity,  and  normality  for  all  3200  ensem¬ 
bles  formed  during  the  10  ms  reverberation  return.  Only  the  last  500 
samples  of  each  ensemble  were  used,  since  the  reverberation  intensity 
was  more  uniform  during  this  last  half  of  the  observation  period.  The 
runs  up  and  down  test  indicated  that  alJ  the  data  were  random  at  a 
two-tailed  level  of  significance  of  0.1.  The  KS-2  test  for  homogeneity 
indicated  that  there  was  some  lack  of  homogeneity  in  the  data  (10  to 
20  percent  of  the  ensembles  of  each  channel  failed  at  a  10  percent 
significance  level)  although  the  inhomogenei ty  war,  judged  not.  to  be  too 
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severe.  Some  sequences  of  ensembles  were  more  nearly  homogeneous  than 
others.  For  example,  ensembles  1801-2200,  which  will  be  examined  in 
some  detail  in  the  next  chapter,  had  less  than  10  percent  failures 
for  all  channels  except  11,  13,  and  llj .  For  these  channels  the  failure 
rate  was  19.25,  10.5,  and  12.0  percent,  respectively. 

The  data  were  also  tested  for  normality.  The  tests  indicated 
some  skewness  and  a  significant  amount  of  kurtosis.  Pearson's  tost 
of  kurtosis  and  D'Agostino's  test  typically  indicated  a  50  percent 
failure  rate.  The  tests  for  normality  were  apparently  not  adversely 
affected  by  the  overall  lack  of  homogeneity.  There  appeared  to  be  no 
correlation  between  ensembles  which  failed  the  homogeneity  tests  and 
those  which  failed  the  tests  for  normality.  Ensembles  1801-2200,  which 
were  generally  homogeneous,  failed  the  tests  for  normality  at  approxi¬ 
mately  the  same  rate  as  the  rest  of  the  ensembles.  Although  previous 

tests  of  surface  reverberation  for  normality  have  indicated  no  signifi- 

Ul 

cant  departures  from  normality,  it  has  not  yet  been  determined  under 
what  conditions  a  departure  from  normality  is  to  be  expected.  That  is, 
it  is  not  known  what  acoustic,  geometric,  or  environmental  conditions 
most  affect  the  statistical  distribution  of  surface  reverberation.  In 
particular,  it  is  possible  that  the  large  surface  wave  heights  encounter¬ 
ed  in  the  present  study  tended  to  result  in  backscattered  sound  with 
larger  amplitudes  and  thus  distributions  of  this  backscattering  with 
larger  than  normal  tails,  which  would  account  for  the  large  kurtosis 
observed.  Thus  it  is  felt  that  the  lack  of  normality  in  no  way  conflicts 
with  previous  results. 


-  * 
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V.  MOMENTS  AND  COVARIANCE 


The  previous  chapter  presented  the  results  of  various  tests  to 
determine  the  validity  of  the  sample  data,  that  is,  to  determine  if 
valid  sample  ensembles  could  be  constructed.  In  addition,  these  ensem¬ 
bles  were  tested  to  determine  if  their  underlying  distribution  was 
gaussian.  Since  the  ensembles  were  shown  to  be  non-gaussian ,  a  closer 
examination  of  the  sample  moments  provides  insights  into  the  type  of 
underlying  distribution  from  which  the  data  were  sampled.  Examples  of 
techniques  whereby  the  sample  moments  are  used  to  approximate  a  distri¬ 
bution  are  Pearson  curves,  Gram-Charlier  series,  and  Cornish-Fisher 
expansions.'’1  Although  the  present  study  will  not  attempt,  to  estimate 
the  distribution  of  the  data,  the  mean,  variance,  skew,  and  kurtosis 
will  be  examined  in  this  chapter  in  some  detail. 

After  an  examination  of  the  moments,  the  spatial  and  temporal 
covariance  of  the  reverberation  will  be  extensively  analysed.  Many 
statistical  techniques  for  the  detection  of  signals  in  the  presence 
of  reverberation  depend  upon  a  knowledge  of  the  covariance  of  the  rever¬ 
beration.  For  this  reason  the  covariance  is  an  important  statistical 
parameter  to  study.  Of  course,  most  of  these  techniques  assume  that 
the  reverberation  is  zero-mean  gaussian,  so  that  the  covariance  of  the 
reverberation  is  the  only  statistical  parameter  necessary  to  completely 
specify  the  distribution  function  of  the  reverberation.  For  non- 
gaussian  distributions,  the  covariance  may  not  be  the  only  parameter 
necessary  to  specify  the  distribution,  and  the  detection  techniques 
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based  on  non-gaussian  reverberation  may  require  more  than  ,iust  a 
knowledge  of  the  covariance.  Nevertheless,  the  covariance  of  the  rever¬ 
beration  remains  an  important  parameter  in  the  statistical  description 
of  the  process. 

This  chapter  is  organized  as  follows.  The  equations  for  the 
numerical  calculation  of  the  mean,  variance,  skew,  kurtosis,  and  co- 
variance  are  given.  The  sum  and  difference  components  of  the  sample  co- 
variance  are  discussed,  along  with  the  technique  for  normalizing  the 
covariance.  The  mean,  variance,  skew,  and  kurtosis  are  then  presented 
and  analyzed  as  a  function  of  time  for  all  13  channels.  The  rest  of  the 
chapter  is  then  an  analysis  of  the  covariance  as  a  function  of  time, 
time-difference,  and  spatial  separation.  The  experimental  results  are 
compared  to  the  theoretical  model.  It  is  found  that  while  the  theoret¬ 
ical  predictions  of  the  covariance  between  the  horizontal  elements 
agree  quite  well  with  the  experimental  results,  the  predictions  of  the 
vertical  covariance  differ  ;uite  significantly  from  the  experimental 
results.  An  explanation  of  this  difference  is  put  forth,  and  changes 
to  the  model  which  would  possibly  bring  the  vertical  predictions  into 
agreement  with  the  experimental  results  are  proposed. 

V. L  Numerical  Computation  of  the  Moments  and  Covariance 

The  rth  population  moment  about  the  mean  of  v  (t),  the  random 
process  from  the  ith  channel  at  time  t,  is  given  by 

Pr ( t , i )  =  ((v.  (t)-p'(t.,i)  )r^ 

where 

M '(t,i  )  =  ^v_.  ( t, 


The  sample  ensemble  estimate  of  u  (t.,i  '  at  the  time  t  (  is 

rar(tn’i}  =  N  £  (vi(V“°  "  nl(,n4 


where 


u.=l 


ral‘ ( tr.  ’ 1  }  =  T  £  Vi(tn’lj) 


and  is  the  sample  estimate  of  the  mean, 

;j  =  number  of  samples  in  the  ensemble,  ana 

v  (t  ,a})  =  the  sample  of  the  reverberation  return  f 
i  n 


channel  at  time  t.  from  ping  w. 


The  sample  estimate  of  the  skew  is 


v/Mt  ,i 


n"'n’"'  1 1 T  , 


“‘O 


( t  ,i) 


and  the  kurtosis 


:t  ,ii  = 


m,  ( t.  ,  i  ) 


mytn,: 


The  numerical  computation  of  the  covariance  is  more  involved 
The  covariance  is  shown  to  consist  of  sum  and  difference  component s . 
Only  the  difference  component  is  considered  in  this  rtuiy.  The 
covariance  is 

K-jj(t,t+t)  =  ^  (t  )-p'(  t. ,  i  )  )  {  v  { t+x  )-u^‘(  (  +r  ,.i  ) 

Writing  the  reverberation  signals  in  their  uundrni  ur<  renresent  a‘  '  ve 


gives 


Ki  j  ( ^  ,t"t’T  )  ~  ^(Xi(t)coswot+Yi(t)sinu)ot-^(i(t)cosa)ot+Y.  (t  )sinuot^) 


• (X  (t+T  )cosw  (t+T  )+Y  ( t+x  ) sinui  (t+x  ) 

d  °  J  O 

-  ( t+T )coswq( t+T )+Yj i t+T )sinuo ( t+T 

Carrying  out  the  multiplication,  utilizing  identities  for  the  product 

of  cosine  and  sine  functions,  and  rearranging  terms  results  in 

k  .  .  ( t  ,t+T )  =  X,  cosoj  t  +  Y  sinui  t 
lj  dodo 

+  X  cos2ui  (t+T/2)  +  Y  sin2ui  (t+r/2) 
so  so 

where  X  =  h  (x. (t)X (t+T)  +  Y.(t)Y  (t+x)^ 

1  J  i  J 

Xs  =  52(x.(t)XJ(t+T)  -  Y.Cth^t+x)) 

Yd  =  ^(^(tjYjCt+x)  -  Y.(t)X.(t+T)) 

Y  =  %<X  (t)Y  (t+x)  +  Y  (t)X.(t+x)> 
and  Xi  =  Xi  ~  ^Xi^  ’  e^c • 

It  is  assumed  that  the  variables  t  and  t  have  no  random  components. 

Thus  the  covariance  can  be  written  as  the  sum  of  a  difference  component 

and  a  sum  component : 

K  , ( t ,t+T )  =  K  +  K 
i  j  d  s 

where 


k,  =  X,  cosrn  x  +  Y,  sinui  x 
d  d  o  d  o 


=  E,  cos(io  T  +  <j>.) 
d  o  d 


and 


k  =  X  cos2ui  (t+x/2)  +  Y  sin2u  (t+x/2) 
s  s  o  so 


=  E  cos2ui  (t+x/2+cji  ) 
so  s 


. 
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The  difference  component  is  so  called  because  of  its  dependence  ■  ip ■ 
the  difference  in  observation  times,  while  the  sum  component  depend:: 
upon  their  sum.  The  envelope  and  phase  of  the  sum  and  difference 
components  can  be  written  in  terms  of  their  quadrature  components  in 
the  lisual  way . 

The  present  study  is  concerned  only  with  the  difference  compon¬ 
ent  of  the  covariance.  It  can  be  shown  that,  in  the  case  of  the  var¬ 
iance,  the  envelope  of  the  sum  component  is  always  Less  than  or  equal 
to  the  envelope  of  the  difference  component.  Since 


a.nd 


2  2  2 
Ed  =  Xd  +  Yd 

E2  =  X2  +  Y2 
s  s  s 


then  E2  -  E2  =  (  XI  ( t  ))(y2 ( t  f)  -  (x.  (t  )Y.  ( t ))  ^ 

where  for  the  variance,  i  =  J  and  t  =  0.  From  Schwarz's  inequal. ty, 
< X.  ( t ) ? . ( t )> 2  <  ( X j ( t j)(F ( t . 

Thus 


2  ? 

E  -  E  >  0 
d  s  — 


and 


E  >  E 
d  —  s 

Therefore  the  difference  component  envelope  is  always  greater  than  or 
equal  to  the  sum  component  envelope.  It  will  be  shown  later  that  the 
sum  component  of  the  variance  was  actually  much  smaller  than  the 
difference  component  for  the  segment  of  the  data  which  was  analyzed. 

It  can  also  be  argued  that  when  estimates  of  the  covariance  are 


computed  from  sample  ensembles,  the  time  variable  t  will  have  a  random 
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component  e  due  to  time  sampling  errors  during  the  A/D  conversion  or 
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motion  of  the  receivers.  ’  If  this  e  is  such  that  m  e  is  uniformly 

o 

distributed  from  -it  to  it,  then  taking  the  expected  value  over  the  cosine 
and  sine  functions  will  cause  the  sum  component  of  the  sample  covariance 
to  be  zero.  Even  if  e  is  not  so  distributed,  any  random  error  in  t 
at  all  will  decrease  the  contribution  of  the  sum  component  to  the  sample 
covariance.  Since  it  appears  that  the  difference  component  of  the 
covariance  is  the  dominant  component,  the  present  study  is  limited  to 
an  examination  of  only  the  difference  component.  Further  references  to 
the  covariance  are  intended  to  include  only  the  difference  component 
unless  otherwise  stated. 

The  sample  estimates  of  the  quadrature  components  of  the  differ¬ 
ence  component  of  the  covariance  are  given  by 


X.jU.t)  =  4 


w)X.  (t+T  ,Ul)  +  Y  .  (  t  ,w)Y  (t  +  T  ,w  )  ) 
J  1  J 


=  jjjj  ^  (X.  (t,u>)X^(t+T,w)  +  Yi(t,u.)Y1(t+T,u)) 


OJ-J. 


-'=[?£;  xi(t'“)ff£  v1”' 

U)=l  10=  ± 


U))  +  -  >  Y.  {  t  ,Ol) 

(Jj-l 


The  normalised  envelope  is  eivi  n  by 


and  the  phase  is  unaffected  by  normalisation. 

The  covariance  is  referred  '  as  a  cross-covariance  for  :  * 


If  i  =  ,1  ,  then  the  covariance  is  an  auto-covariance. 


An  aul  o-covari  nr.co 
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at  i  =0  is  just  the  variance.  The  norma] i zed  covariance  is  sometimes 
referred  to  as  the  correlation  coefficient. 

The  quadrature  components  of  the  reverberation  returns  were 
obtained  by  sampling  the  received  waveforms  directly,  i.e.. 


and 


X.  ( 1.  ,ui)  =  v.  ( t  ,io ) 

1  1 

i  ( t. ,  u> )  =  v.  (t  +  ,  <->) 


where  f  is  the  center  frequency  of'  th  transmitted  signal.  The 

method  of  obtaining  (t,u>)  directly  from  the  received  waveform  is  an 
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approximation  consistent  wit.n  the  narrowband  assumption. 


V.2  Moments  of  the  Reverberation  Processes 

The  sample  mean,  variance,  skew,  and  kurtosis  are  plotted  for 
200  ensembles,  representing  a  time  duration  of  0.622  ms,  centered  about 
a  time  after  transmit  of  75-25  ms.  Much  of  the  analysis  of  the  covar¬ 
iance  is  also  centered  about  this  time.  The  short  time  duration  of 
the  plots  allows  the  oscillatory  nature  of  these  moments  to  be  observed. 
V.2.1  Sample  Mean 

The  sample  means  for  all  13  channels  are  shown  in  Figs. 
V-l  through  V-4.  Since  an  absolute  calibration  of  the  receivers  was 
not  maintained  throughout  the  recording  and  digitization  process,  the 
magnitudes  of  the  means  have  little  significance.  That  is,  gain  dif¬ 
ferences  which  occurred  between  channels  during  recording  or  digitiza¬ 
tion  will  be  reflected  in  the  computed  means.  Thus  these  plots  repre¬ 
sent  the  sample  mean  only  to  within  an  undetermined  scale  factor  which 
is  different  for  each  channel.  The  same  is  true  of  the  sample  variance 
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which  is  shown  later..  However,  the  ratio  of  the  mean  to  the  square- 
root  of  the  variance  is  independent  of  any  scale  factors.  In  spite  of 
this  limitation,  it  is  still  possible  to  estimate  the  statistical 
significance  of  the  amplitude  of  the  mean.  The  ratio  of  the  mean  to 
the  square  root  of  the  variance  provides  one  measure  of  the  significance 
of  the  mean.  For  the  data  shown  here,  a  typical  value  for  this  ratio 
is  .07.  The  Student's  t-test  for  zero-mean  of  normal  samples  indicates 
that  for  this  ratio  the  mean  can  be  accepted  as  zero  at  a  1 2%  signifi¬ 
cance  level.  Although  the  data  are  not  normal,  nevertheless  this  is 
some  indication  that  the  mean  is  not  statistically  large. 

Another  measure  of  the  statistical  significance  of  the 
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sample  mean  is  its  sampling  variance,  given  by 
var ( |  ng 

The  sampling  variance  is  a  measure  of  the  spread  of  the  sample  estimate 
of  a  statistical  parameter  (in  this  case  the  sample  mean)  about  its 
true  value.  Using  the  sample  estimate  m0  in  place  of  the  population 
moment  it  was  found  that  most  of  the  values  of  the  mean  fell  within 
plus  or  minus  one  standard  deviation  of  zero.  Essentially  all  of  the 
values  were  within  two  standard  deviations  of  zero.  This  also  indicates 
that  the  sample  mean  is  not  statistically  large. 

Even  though  the  sample  mean  is  not  statistically  differ¬ 
ent  from  zero,  observation  of  Figs.  V-l  through  V-U  .'hows  that  the  means 
are  not  completely  random  about  zero  but  rather  oscillate  at  the  trans¬ 
mitted  frequency  u^.  This  is  to  be  expected  from  a  bandlimited  signal 


if  the  random  component  of  t  is  small  since 
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Thus  the  mean  of  a  handlimited  signal  ir.  also  hand]  imi ted  and 
at  the  frequency  w  .  Any  random  component  of  t  will  have  the  of 
reducing  the  amplitude  of  the  mean.  liowc/er,  it  is  uppain  at  in 
ent  data  that  this  random  component  is  sufficiently  small  tint  i 
oscillatory  fact. or  is  preserved.  Therefore,  it  can  r>e  conclude.', 
although  from  a  statistical  point,  of  view  the  neon  is  small,  n.  v 
it  contains  a  periodic  component,  at  the  center  frequency  <•>  . 


V.2.2  flump  le  Variance 

The  sample  variance  for  all  1  -  channels  is  shown 
Figs.  V-5  through  V-8.  The  actual  sample  variance  is  plotted  al 
with  the  difference  component  of  the  sample  variance.  It  can  b» 
that  the  variance  does  contain  a  component  at  a  frequency  of  ho, 
that  it  is  small  compared  to  the  difference  component. .  i'he  same 
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the  sampling  variance  can  be  estimated  as 
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FIGURE  V -5 

COMPARISON  OF  ACTUAL  VARIANCE  AND 
DIFFERENCE  COMPONENT,  CHANNELS  1  -3 
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FIGURE  V -6 

COMPARISON  OF  ACTUAL  VARIANCE  AND 
DIFFERENCE  COMPONENT,  CHANNELS  4-6 
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FIGURE  V-7 

COMPARISON  OF  ACTUAL  VARIANCE  AND 
DIFFERENCE  COMPONENT,  CHANNELS  7-9 
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FIGURE  V-9 

SAMPLE  SKEW,  CHANNELS  1-3 


ARL  UT 
AS-81  '0 
GRW  (3A 
1  9  81 


CHANNEL  6  CHANNEL  5  CHANNEL  4 


s 


ENSEMBLE  NUMBER 


r - ■— 1 

94  75 


TIME  AFTER  TRANSMIT  -  msec 


FIGURE  V  II 

SAMPLE  SKEW,  CHANNELS  7-9 
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V.2.1*  Sample  Kurtonls 


li? 


Plots  of  the  sample  kurtosis  are  shown  In  Figs.  V—  1  ■? 

through  V-l6.  Oscillations  are  evident  and  appear  to  be  at  a  frequency 

of  20)^,  although  they  arp  probably  at  i*u)Q.  A  normal  dl  stribut.  Ion  has 

a  kurtosis  of  3.  These  samples  are  generally  significantly  larger  than 

3.  The  10#  confidence  limits  for  Pearson's  test  of  kurtosis  are  P.67  - 
1*6 

3.37.  Most  of  the  values  of  the  kurtosis  are  above  the  upper  limit,. 
The  values  of  the  kurtosis  are  not  distributed  symmetrically  about, 
the  normal  value  of  3.  Rather,  the  kurtosis  is  on  the  average  larger 
than  3.  Thus  the  periodic  component  of  the  kurtosis  is  not  the  dominant 
factor  in  causing  a  rejection  of  normality,  as  was  the  case  for  the 
skew. 

V. 3  The  Covariance 

This  section  examines  only  the  difference  component  of  the  co- 
variance.  Both  the  sample  estimates  from  the  experimental  data  and 
the  theoretical  results  from  the  point-scatter  model  are  presented. 

This  section  is  developed  as  follows.  The  difference  component  of  the 
variance  over  the  entire  10  ms  digitization  interval  is  discussed  first. 
Then  the  auto-covariance  as  a  function  of  t  is  presented  and  analyzed. 

The  cross-covariance  as  a  function  of  spatial  separation  for  x=0  is 
then  compared  for  the  horizontal  and  vertical  arrays.  The  covariance  for 
r^O  is  presented  and  the  section  concludes  with  a  discussion  of  the 
failure  of  the  theoretical  model  to  adequately  predict  the  spatial 
dependence  of  the  covariance  for  the  vertical  array. 

V.3.1  The  Difference  Component  of  the  Variance 

Figures  V-17  through  V-20  are  plots  of  the  sample 
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difference  component  of  the  variance  over  the  entire  10  ms  of  the 
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SAMPLE  KURTOSIS,  CHANNELS  13 
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FIGURE  V  14 

SAMPLE  KURTOSIS,  CHANNELS  4-6 
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FIGURE  V-15 

SAMPLE  KURTOSIS,  CHANNELS  7-9 
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FIGURE  V-17 

COMPARISONS  OF  EXPERIMENTAL  AND  THEORETICAL  VARIANCES 
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FIGURE  V-18 

COMPARISONS  OF  EXPERIMENTAL  AND  THEORETICAL  VARIANCES 
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FIGURE  V-19 

COMPARISONS  OF  EXPERIMENTAL  AND  THEORETICAL  VARIANCES 
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digitized  reverberation  returns  for  all  13  channels.  Also  plotted  is  the 
theoretical  difference  component  of  the  variance.  Good  agreement  is 
attained  between  the  sample  and  theoretical  estimates.  The  fluctuations 
of  the  sample  variance  can  be  attributed  to  sampling  errors  due  to  the 
finite  ensemble  size. 

Since  the  variance  represents  the  intensity  of  the 
reverberation,  most  of  the  following  analysis  of  the  covariance  was 
performed  at  a  time  when  the  variance  was  the  largest  in  order  to  max¬ 
imize  the  reverberation-to-ambient  noise  ratio.  As  can  be  seen  from 
the  plots,  the  reverberation  intensity  increases  to  a  maximum  during 
this  10  ms  and  then  begins  to  decrease  slightly  by  the  end  of  the  10  ms 
data  segment.  This  change  in  the  intensity  is  due  primarily  to  the 
vertical  directivity  of  the  projector.  The  time  chosen  for  most  of 
the  covariance  analysis  was  75-25  ms  after  transmit,  which  is  approx¬ 
imately  the  maximum  of  the  reverberation  intensity  on  each  channel. 

V.3.2  Auto-covariance 

The  sample  estimates  of  the  normalized  envelopes 
e. . (t,r)  and  phases  4>ii(t ,t )  of  the  auto-covariances  for  the  horizontal 
and  vertical  arrays  are  shown  in  Figs.  V-21  and  V-22.  Observation  of 
the  auto-covariance  permits  examination  of  the  dependence  of  the  co- 
variance  on  t  independently  of  its  dependence  upon  spatial  separation.  The 
normalized  auto  •  covariance  envelope  has,  by  definition,  a  maximum  of 
1  at  t— 0 .  Also  an  examination  of  the  equations  for  the  phase  of  the 
difference  component  shows  that  the  phase  must  be  zero  at  t=0. 

It  can  be  observed  that  the  auto-covariance  decreases  to 
approximately  zero  for  x  =  +250  us.  It  is  usually  stated  that  the 
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FIGURE  V-21 

NORMALIZED  ENVELOPES  AND  PHASES  OF  SAMPLE  ESTIMATES 
OF  THE  AUTOCOVARIANCE  OF  THE  HORIZONTAL  RECEIVERS 
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FIGURE  V-22 

NORMALIZED  ENVELOPES  AND  PHASES  OF  SAMPLE  ESTIMATES 
OF  THE  AUTOCOVARIANCE  OF  THE  VERTICAL  RECEIVERS 
AS  A  FUNCTION  OF  r  FOR  t  -  75.25  miac 
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.  "ivan.).if,i?  decreases  zer. ;  i  ,j  a  i  i me  1  .i*  erval  equal  to  or.**  pulse 
length.  However ,  fne  covavi ance  envelope  has  a  value  of  about  0.5  at 
the  time  interval  of  100  us,  the  pulse  length  use.i  in  this  study.  This 
apparent  discrepancy  can  be  resolved  by  pointing  out  that  the  usual 
statement  is  valid  only  for  extremely  narrowband  signals,  i.e.,  Long 
pulse  lengths.  For  wider  bandwidths,  the  frequency  responses  of  the 
projector  and  receiver  filter  the  transmitted  signal  more  ;  everely  and 
result  in  an  increase  in  the  duration  of  tne  covariance  foejond  the  time 
interval  of  one  pulse  length.  The  transmitted  pulse,  after  being  filtered 
by  the  projector  and  receivers,  is  represented  by  the  waveforms  v^(t  ) 
arid  is  shown  in  Figs.  TIT-1  and  ITT-.’  for  all  thirteen  receivers.  It 
can  be  seen  that  the  envelope  of  the  transmitted  pulse  has  a  duration 
of  approximately  275  PS  after  being  filtered.  Thus  it  appears  that  the 
covariance  decreases  to  zero  in  a  time  interval  approximately  equal  to 
the  length  of  the  transmitted  pulse  after  it  has  been  filtered  by  the 
projector  and  receivers. 

This  effect  has  been  observed  by  others.  ’  Figure 
V-2'j  shows  representative  envelopes  of  the  auto-covariance  and  the 
filtered  transmitted  signal  for  one  receiver  from  Frazer'1”  and  from  the 
present  study.  In  both  cases  the  pulse  length  was  100  ps;  however, 
after  being  filtered,  the  durations  of  the  transmitted  pulses  were 
di f fererit .  The  time  interval  of  the  auto-covariances  are  also  different 
and  in  both  cases  correspond  approximately  to  the  length  of  the  filtered 
transmitted  pulses.  This  agreement  tends  to  confirm  the  conclusion 
that  the  time  interval  of  the  covariance  is  related  to  the  filtered 
transmitted  pulse  length.  It,  also  points  out  that  a  valid  theory  of 


the  covariance  of  surface  reverberation  needs  to  take  into  account-  tin: 
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A  COMPARISON  OF  A  REPRESENTATIVE  WAVEFORM  Vj(tr) 
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bantwiclt  I.  rml  pulse  length  of  the  transmitted  signal  and  the  frequency 
r  'spouses  ■  ■  f  t  he  projector  and  receivers. 

The  symmetry  of  the  auto-covariance  envelopes  and  phases 
in  Figs.  V- 21  and  V-22  is  also  of  note.  The  envelopes  have  approximately 
even  symmetry  about  t=0  while  the  phases  exhibit  an  odd  symmetry.  These 
conditions  are  sufficient  to  specify  the  reverberation  as  being  locally 
stationary,  at  least  at  this  observation  time  t.  A  random  process  is 
locally  stationary  at  the  time  t  if  the  auto-covariance  is  an  even 
function  of  t,  i.e., 

i ( t  ,x  )  =  kfi (t ,-t  ) 

The  local  stationarity  of  reverberation  from  short  pulsed  CW  trans¬ 
missions  has  also  been  previously  demonstrated.1^ 

Figures  V-2U  and  V-25  are  similar  to  Figs.  V-21  and  V-22 
but  in  addition  include  the  theoretical  estimates  of  the  envelopes  and 
phases  of  the  auto-covariances.  As  can  be  seen,  excellent  agreement 
between  the  theoretical  and  experimental  estimates  is  obtained  for  both 
the  envelopes  and  the  phases.  This  is  true  of  both  the  horizontal  and 
vertical  arrays.  In  particular,  the  time  extent  of  the  auto-covariance 
predicted  by  the  model  Is  in  good  agreement  with  the  experimental  re¬ 
sults.  Thus  it  appears  that  at  least  in  this  case  the  model  is  correctly 
taking  into  account  the  bandwidth,  pulse  length,  and  frequency  responses 

of  the  projector  and  receivers.  Similar  excellent  agreement  has  been 

1  7 

obtained  previously  for  a  horizontal  array  using  the  same  model. 

This  confirmation  provides  increased  confidence  in  the  model  for  the 
cases  examined. 

There  appears  to  he  some  similarity  between  the  phases  of 

the  auto-covariance  and  the  phases  of  the  waveforms  V.(t  )  in  Figs.  TII-1 
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FIGURE  V-24 

COMPARISONS  OF  NORMALIZED  ENVELOPES  AND  PHASES  OF  SAMPLE 
ESTIMATES  AND  THEORETICAL  ESTIMATES  OF  THE  AUTOCOVARIANCE 
OF  THE  HORIZONTAL  RECEIVERS  AS  A  FUNCTION  OF  r  FOR  t  =  75.25  msec 
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FIGURE  V-26 

COMPARISONS  OF  NORMALIZED  ENVELOPES  AND  PHASES  OF 
SAMPLE  ESTIMATES  AND  THEORETICAL  ESTIMATES  OF  THE 
AUTOCOVARIANCE  OF  THE  VERTICAL  RECEIVERS 
AS  A  FUNCTION  OF  t  FOR  t  -  75.26  m$ec 
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and  III-2.  The  similarity  is  more  apparent  for  the  theoretical 
estimates  and  to  a  lesser  extent  for  the  experimental  estimates 


The 


actual  magnitudes  of  the  phases  of  V^(t  )  are  not  reflected  in  the  phases 

of  the  auto-covariances,  but  rather  the  change  in  phase  during  the 

extent  of  V.(t  ).  The  shift  of  the  phases  of  V.(t  )  for  the  different 
l  r  1  r 

receivers  of  the  vertical  array  represents  the  slight  path  length 
differences  from  the  sphere  to  the  vertical  receivers;  the  horizontal 
receivers  exhibit  no  such  shifts  since  their  closer  spacing  resulted 
in  insignificant  path  length  differences  when  pointed  directly  at  the 
sphere.  This  shift  in  the  phases  for  the  vertical  array  is  not  observable 

in  the  phases  of  the  auto-covariances  of  the  vertical  array.  However, 

the  changes  in  phase  during  the  extent  of  the  waveforms  V..  (t^ )  can  be 

correlated  with  similar  changes  in  phase  during  the  extent  of  the  auto- 

covariances  of  the  same  channels.  Considering  the  dependence  of  the 

theoretical  model  upon  the  waveforms  V. (t  ),  the  similarity  of  the 

i  r 

phases  of  V^ft^)  and  the  phases  of  the  theoretical  estimate  of  the 
auto-covariances  is  not  too  surprising.  Some  of  the  experimental  esti¬ 
mates  of  the  phase  appear  to  also  have  shapes  similar  to  their  corre¬ 
sponding  phases  for  V^(t  ),  but  some  do  not.  Although  Frazer  has  reported 

a  noticeable  similarity  between  the  experimental  phases  of  the  auto- 

IT 

covariance  and  the  phases  of  iMt^),  the  present  data  do  not  con¬ 
clusively  show  this  similarity. 

To  this  point,  the  auto-covariances  have  been  examined 
as  a  function  of  r  only  at  one  time  t.  Figure  V-26  shows  the  auto¬ 
covariance  for  channel  1  as  a  function  of  both  t  and  x  over  the  entire 
10  ms  of  digitized  reverberation.  Both  the  unnormalized  envelope 
e.  (t,x)  and  the  normalized  envelope  e  ^t,!)  are  shown,  along  with 
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FIGURE  V-26 

UNNORMALIZED  AND  NORMALIZED  ENVELOPES  AND  PHASES 
OF  SAMPLE  ESTIMATES  OF  THE  AUTOCOVARIANCE  OF 
CHANNEL  1  AS  A  FUNCTION  OF  t  AND  r 
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the  phase  <t>,  ^t,!).  The  unnormalized  envelope  demonstrates  the  change 
in  reverberation  intensity  during  the  10  ms  interval.  The  maximum  of 
the  unnormalized  covariance  increases  with  time  for  most  of  the  interval 
before  it  begins  to  decrease  again  toward  the  end.  This  change  has  been 
observed  in  the  variance  shown  previously.  The  change  is  due  to  the 
vertical  directivity  of  the  projector.  That  is,  the  reverberation 
returns  from  the  earlier  times  are  caused  by  scattering  from  the  surface 
which  is  illuminated  from  an  off-axis  angle  of  the  projector's  vertical 
directionality  response.  Thus  the  reverberation  intensity  is  less.  As 
time  increases,  the  surface  which  is  illuminated  by  the  on-axis,  or 
maximum,  vertical  response  of  the  projector  begins  to  contribute  to  the 
reverberation  return,  resulting  in  an  increase  in  reverberation  intensity. 
The  grazing  angle  between  the  projector  and  surface  at  the  earliest  time 
of  69  ms  is  11.2  degrees,  decreasing  to  9.8  degrees  at  79  ms.  The  change 
in  angle  is  l.U  degrees,  which  is  less  than  the  vertical  beamwidth  of 
the  projector.  Thus  the  entire  digitized  reverberation  return  was 
caused  by  scattering  only  from  the  main  lobe  of  the  projector.  It 
should  be  noted  that  spreading  losses  and  the  vertical  directionality 
of  the  receiver  will  also  affect  the  change  in  reverberation  intensity 
with  time.  However,  these  effects  were  insignificant  for  the  period 
of  time  which  is  observed  here. 

This  change  in  the  covariance  with  time  points  out  one 
form  of  non-stationarity  which  reverberation  characteristically  has. 

A  random  process  is  said  to  be  wide-sense  stationa’y  in  an  interval  if 
its  auto-covariance  is  a  function  only  of  t  within  an  interval  of  tine. 


1  ,  r» . 


k..(t,T)  =  k..(x)  ,  1.  •  t  •'  t.. 

11  n  i  - 

It  is  obvious  from  Fig.  V-2 6  that  the  unnomal ized  autc-er.  var i a  nee 
changes  with  time;  thus  it  is  not  wide-sense  stationary.  However,  the 
unnormalized  auto-covariance  does  appear  tc  represent  a  local iy 
stationary  process  at  each  time  due  to  the  symmetry  of  its  envd<  -a  at. 
phase  about  r  at  each  time.  Even  though  the  unnormal i red  covari  an '•? 
is  not  wide-sense  stationary,  it  does  exhibit  one  property  of  vide-:-' 
stationary  process:  the  covariance  tnr  it.s  maximum  at  r=, ,  -w r 
time  t.  On  the  other  haiid,  the  normalized  auto-covariance  does  aprenv 
to  represent  a  wide-sense  stationary  process  over  the  interval 
69  ms  <_  t  <_  79  ms  since  the  normalized  envelope  and  phase  chan,  e  lit4  1  •' 
in  this  interval.  If  the  unnormalized  aut o-covar i once  k  .  ( l , : )  repre¬ 
sents  a  random  process  X^(t),  then  the  normalized  auto-covar i an ,?o 
represents  a  normalised  random  process  where 

i 


x.(t)  = 


X.  (t) 


k  .  .  ( t ,  0 ) 


Thus  this  technique  for  normalizing  the  random  process  X  ( •  '  has,  in 
this  instance,  removed  much  of  its  nori-stationar  i  ry .  This  need  p  d 
be  the  case  in  general  though.  In  particular,  "jg.  V-.' !  i  t*p  i  %:.i  r. 
the  difference  component  of  the  covariance.  The  act.ua!  v.ar' a:.c-  'v 

be  non-  s  tat  i  onary ,  even  when  normalized,  due  to  contributions  •’run 

,,  .  5  6 

the  sum  term. 

V.3.3  The  Covariance  for  1=0 

The  previous  sections  have  analyzed  t  ho  temporal  -h  : 


dence  of  the  variance  and  auto-eovur i ance . 


The  sn.at  i  a  I  : 


icO 


the  covariance  is  addressed  in  this  section,  and  the  joint  spatial- 
temporal  covariance  is  discussed  in  the  following  section.  The  spatial 
covariance  for  the  horizontal  array  is  presented  first,  followed  by  the 
vertical  spatial  covariance. 

The  envelope  of  the  normalized  covariance  between  the 
horizontal  elements  at  x=0  as  a  function  of  spatial  separation  of  the 
elements  is  shown  in  Fig.  V-27.  The  envelope  of  the  covariance 
decreases  rather  significantly  within  a  spatial  separation  of  just  a 
few  wavelengths.  In  fact,  the  envelope  has  decreased  to  a  value  of  0.1 
in  just  over  four  wavelengths. 

It  is  of  interest  to  quantify  the  error  associated  with 
the  sample  estimate  of  the  covariance  due  to  the  finite  sample  size. 

The  sampling  variance  of  the  normalized  covariance  is  given  by"’’*' 


where 


and 


Urs  =  <(vry{(i))r(v.-u'(j))S) 


Ki,j  N /  U20U02 


The  time  dependence  has  been  suppressed  for  compactness  of  notation, 
either  r  or  s  is  zero,  then  the  bivariate  population  moment  p 
reduces  to  its  corresponding  univariate  moment.  The  usual  procedure 
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FIGURE  V- 27 

COMPARISON  OF  THE  THEORETICAL  AND  EXPERIMENTAL  DEPENDENCE 
OF  THE  ENVELOPE  OF  THE  NORMALIZED  COVARIANCE  ON  SPATIAL 
SEPARATION  FOR  THE  HORIZONTAL  ARRAY 
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when  the  population  moments  are  unknown  is  to  replace  the  population 
moments  with  their  sample  estimates  mrg  in  order  to  evaluate  the  sample 
variance.  In  the  present  study,  the  sample  moments  ^2'  m31’  an(*  m13 
have  not  been  calculated.  However  if  the  Joint  random  processes  are 
Jointly  normal,  then 

v,r(ki,j)  =  if 

which  can  be  approximated  as 

var(k  )  a  i  (l-k2  )2 

i  iJ  N  i ,  J 

The  standard  deviation  of  the  sample  estimate  of  the  normalized  co- 
variance  is  just  the  square  root  of  the  sampling  variance.  Even  though 
the  present  data  are  not  normal,  this  normal  approximation  was  used 
to  quantify  the  sampling  errors. 

The  error  bars  in  Fig.  V-27  represent  one  standard  de¬ 
viation  of  the  sample  estimate  above  and  below  the  measured  envelope  of 
the  covariance.  Since  the  envelopes  of  the  covariance  from  the  four 
largest  separations  axe  all  within  one  standard  deviation  of  the  value 
of  0.1,  no  statistical  significance  can  be  placed  upon  the  changes  in 
the  envelope  at  these  separations.  However,  it  is  significant  to  note 
that  the  envelope  remains  at  0.1  instead  of  going  to  zero  at  these 
separations. 

Theoretical  estimates  of  the  envelope  of  the  covariance 
were  made  at  the  first  three  separations  and  are  also  shown  in  Fig.  V-27. 
It  can  be  seen  that  good  agreement  is  obtained  at  these  separations 
between  the  experimental  and  theoretical  estimates.  Theoretical  esti- 


mates  were  not  made  at  the  larger  separations  due  to  the  significant 
increase  in  computer  cost  required  to  make  the  numerical  integrations 
converge  at  these  separations. 

Similar  agreement  between  the  theoretical  and  experi¬ 
mental  estimates  of  the  horizontal  covariance  has  been  obtained  by 

17 

Frazer  using  the  same  theoretical  model.  However,  Frazer  obtained 
larger  values  of  covariance  than  were  obtained  in  the  present  study. 

The  theoretical  and  experimental  results  from  Frazer  and  the  present 
study  are  compared  in  Fig.  V-28.  At  a  separation  of  0.05  m,  for  example, 
Frazer  obtained  a  covariance  envelope  of  0.6  while  0.3  was  obtained  in 
the  present  study.  It  is  significant  to  note  that  the  theoretical 
model  correctly  predicted  the  envelope  of  the  covariance  in  both  cases, 
even  though  the  values  were  significantly  different.  This  is  an  indi¬ 
cation  that  the  model  is  correctly  taking  into  account  the  experimental 
differences  which  caused  these  significantly  different  results.  It  is 
felt  that  the  wider  horizontal  transmit  beamwidth  used  in  the  present 
study  was  the  primary  factor  which  resulted  in  the  smaller  values  of 
covariance.  The  present  study  employed  a  horizontal  transmit  beamwidth 
of  15  degrees,  compared  to  10  degrees  used  in  Frazer's  experiment.  The 
wider  beamwidth  insonifies  more  "off-axis"  scatterers  on  the  surface, 
which  tends  to  decrease  the  covariance  between  horizontal  receivers. 

Thus  these  results  from  the  present  study  tend  to  confirm  the  results 
obtained  by  Frazer. 

Figure  V-27  showed  the  covariance  as  a  function  of 
spatial  separation  at  only  one  time  t.  It  is  also  of  interest  to  know 
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FIGURE  V-28 

COMPARISON  OF  THE  THEORETICAL  AND  EXPERIMENTAL  DEPENDENCE  OF  THE 
ENVELOPES  OF  THE  NORMALIZED  COVARIANCE  ON  SPATIAL  SEPARATION 
FOR  THE  HORIZONTAL  ARRAYS  OF  WILSON  AND  FRAZER 
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if  there  is  any  significant  change  in  the  spatial  covariance  over  the 
10  ms  observation  interval.  Figure  V-29  shows  the  envelopes  of  the 
normalized  covariance  at  r=0  as  a  function  of  time  for  the  six  separa¬ 
tions  of  the  horizontal  array.  It  can  be  observed  that  the  envelope 
of  the  covariance  does  not  change  significantly  over  the  10  ms  interval 
for  each  separation.  The  random  fluctuations  which  are  present  are 
all  less  than  two  standard  deviations  of  the  mean  value  of  the  envelope 

over  the  interval;  thus  any  apparent  changes  can  be  attributed  to 
sampling  errors  due  to  the  finite  sample  size  of  the  ensembles.  Since 
the  envelopes  do  not  change  significantly  with  time,  the  spatial  depen¬ 
dence  demonstrated  in  Fig.  V-27  for  one  time  will  be  approximately 
the  same  for  all  other  times  in  the  10  ms  interval. 

The  actual  covariance  is  determined  not  only  by  the 
envelope  but  also  by  the  phase.  Figure  V-30  contains  plots  of  the  phase 
of  the  covariance  at  t=0  as  a  function  of  time  for  each  of  the  six 
separations  of  the  horizontal  array.  The  most  obvious  feature  to  note 
is  that  the  separations  which  have  significantly  non-zero  values  for 
the  envelopes  also  have  phases  which  are  constant  with  time,  whereas 
the  separations  with  small  envelope  values  appear  to  have  phases  which 
change  approximately  linearly  with  time.  Specifically,  4^  -^(t.O)  and 
$11  .^(tjO),  representing  separations  of  2.5  and  5.0  cm  respectively, 
are  reasonably  constant  with  time,  while  the  other  phases  change  sig¬ 
nificantly.  It  can  also  be  observed  that  the  four  phases  which  change 
are  all  very  similar,  i.e.,  they  all  begin  at  approximately  180  degrees 
and  decrease  linearly  at  about  the  same  rate. 
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A  non-zero  phase  implies  that  the  effective  scattering 

center  is  off-axis  azimuthally ,  indicating  a  non-uniform  surface  scatter 

density  in  the  azimuthal  angle.  A  change  in  the  phase  of  the  horizontal 

covariance  with  time  is  indicative  of  a  change  in  the  bearing  angle 

to  the  effective  scattering  center  with  range;  thus  the  scatter  density 

changes  with  range  as  well.  However,  it  is  curious  that  this  change  in 

phase  with  range  is  only  noticeable  at  the  larger  separations,  which 

have  smaller  values  of  covariance.  One  possible  explanation  is  that 

some  sort  of  secondary  scattering  is  occurring  which  is  causing  the 

effects  described  here,  but  is  observable  only  when  the  covariance  from 

the  primary  scattering  is  negligible.  However,  at  this  point  no 

mechanism  for  secondary  scattering  which  would  account  for  these  effects 

is  known.  It  should  be  noted  that  this  change  in  phase  with  time  was 
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also  observed  by  Frazer.  Also  of  note  is  the  fact  that  while 
$13  and  12(t,0)  are  constant  with  time,  they  are  also 

non-zero,  indicating  a  non-uniform  surface  scatter  density.  Thus 
for  the  significant  levels  of  covariance,  the  actual  covariance  is 
somewhat  less  than  the  envelope  of  the  covariance  and  does  not  change 
significantly  with  time  over  the  observed  interval. 

Theoretical  estimates  of  the  envelope  and  phase  of  the 
normalized  covariance  as  a  function  of  time  for  t=0  have  also  been  made 
at  three  horizontal  separations ,  and  are  compared  with  the  experimental 
results  in  Figs.  V-31  and  V-32.  The  theoretical  envelope  values  are 
almost  constant  throughout  the  10  ms  interval.  The  envelopes  t,0 ) 

and  ^(t,0)  are  consistently  slightly  larger  than  the  experimental 
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values,  but  the  agreement  is  still  good.  E12  13(t,o)  agrees  very  well 
with  the  experimental  results. 

The  theoretical  phases  are  also  constant  over  the  in  ms 
interval.  In  fact,  the  phases  are  all  very  close  to  zero  for  the  three 
separations.  This  is  to  be  expected  since  the  model  assumes  a  uniform 
scatter  density  at  the  surface.  As  has  already  been  stated,  the  phases 

^11,12^*^  and  <^2  -^(t,*})  are  also  constant  with  time  but  different 
from  zero.  Thus  the  covariance  predicted  by  the  model  will  be  greater 
than  the  experimental  covariance,  even  though  the  theoretical  and  experi¬ 
mental  envelopes  agreed  well,  because  the  magnitudes  of  the  experimental 
phases  were  greater  than  the  theoretical  phases.  It  should  also  be  noted 
that  the  mechanism  causing  <J>12  ^(tjO)  to  change  with  time  is  completely 
unaccounted  for  in  the  model. 

The  previous  discussion  has  examined  only  the  horizontal 
spatial  dependence  of  the  covariance.  In  a  similar  manner,  the  spatial 
dependence  of  the  covariance  for  the  vertical  array  will  now  be  pre¬ 
sented.  Figure  V-33  displays  the  envelope  of  the  normalized  covariance 
at  x=0  as  a  function  of  spatial  separation  for  the  vertical  array. 

Unlike  the  horizontal  array,  some  pairs  of  elements  have  identical 
spacings.  Thus  more  than  one  value  of  the  envelope  of  the  covariance 
occurs  at  some  separations.  It  is  of  interest  to  determine  if  the 
covariance  at  a  specified  separation  is  dependent  upon  the  position  of 
the  elements  on  the  vertical  array.  As  an  indication  of  this  dependence, 
error  bars  representing  one  standard  deviation  of  the  sample  covariance 
above  and  below  the  mean  envelope  value  at  each  separation  are  provided 


on  the  plot  In  Fig.  V-33.  Most  values  of  the  envel  ope  are  Mi  *  ain  *'  * ■ r 
standard  deviation  of  the  mean  value  at  each  separation.  They  are  ail 
certainly  within  two  standard  deviations.  Any  differences  between  the 
various  values  of  the  envelopes  at  each  separation  are  not  statistically 
significant  and  can  be  attributed  to  sampling  errors  due  to  a  finite 
ensemble  size.  Thus  the  envelope  of  the  covariance  is  no4  affected 
by  the  position  of  the  pair  of  elements  on  this  vertical  array, 

The  covariance,  though,  is  determined  not  only  by  its 
envelope,  but  also  by  its  phase.  If  the  phase  of  the  covariance 
changes  with  respect  to  the  position  on  the  array,  then  the  covariane= 
will  also.  Figure  V-3^  shows  the  covariance  as  a  function  of  separation 
on  the  vertical  array.  It  can  be  seen  that  the  covariances  in  general 
differ  significantly  from  their  envelope  values,  indicating  a  non-zero 
phase.  Also,  there  is  a  greater  range  of  values  of  the  covariance  at 
separations  which  occur  at  more  than  one  location  on  the  array.  Tims 
the  covariance  does  change  with  position  on  the  array.  Furthermore, 
this  change  can  be  attributed  to  the  change  of  rhe  phase  of  the  co¬ 
variance  with  position.  The  phase  will  be  examined  in  more  detail  later. 

However,  the  most  striking  observation  about  the  ver* i cal 
covariance  is  that  it  maintains  a  significant  level  at  much  larger 
separations  than  the  horizontal  covariance.  Figure  V-3f>  compares  the 
vertical  and  horizontal  envelopes  of  the  covariance.  The  values  shown 
for  the  vertical  array  represent  the  mean  value  of  the  envelope  at 
each  separation.  The  envelope  of  the  vertical  covariance  decreases 
monotonically  with  increasing  separation,  but  even  the  covariance  at  a 
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FIGURE  V-34 

THE  DEPENDENCE  OF  THE  EXPERIMENTAL  ESTIMATE  OF  THE  NORMALIZED 
COVARIANCE  ON  SPATIAL  SEPARATION  FOR  THE  VERTICAL  ARRAY 
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separation  of  greater  than  50  wavelengths  has  an  envelope  value  greater 
than  0.2.  By  comparison,  the  envelope  of  the  horizontal  covariance 
has  decreased  to  the  same  value  in  approximately  3  wavelengths.  As 
another  example,  the  envelope  of  the  vertical  covariance  has  a  value 
of  0.5  at  approximately  25  wavelengths,  whereas  the  envelope  of  the 
horizontal  covariance  has  this  value  at  approximately  2  wavelengths. 

The  vertical  covariance  was  also  computed  from  the 
theoretical  model.  Figure  V-3 6  compares  the  envelopes  of  the  experi¬ 
mental  and  theoretical  estimates  of  the  normalized  covariance  at  t=0. 

As  can  be  seen  the  theoretical  predictions  are  not  at  all  in  agreement 
with  the  experimental  estimates.  The  model  did  predict  a  covariance 
which  was  larger  for  the  vertical  array  than  the  horizontal  array,  but 
the  covariance  it  predicted  for  the  vertical  array  was  much  larger 
than  what  was  actually  obtained  experimentally.  Considering  how 
well  the  model  predicted  the  dependence  of  the  horizontal  covariance 
on  separation,  and  the  dependence  of  both  the  vertical  and  horizontal 
auto-covariance  on  time  and  time  difference  (t  and  t),  its  failure  for 
the  vertical  covariance  was  quite  unexpected.  However,  a  close  exam¬ 
ination  of  the  simplifying  assumptions  which  were  introduced  into  the  model 
revealed  the  probable  reasons  for  the  failure  of  the  model  in  this  case. 
These  reasons  will  be  discussed  qualitatively  later  in  the  chapter, 
along  with  suggested  improvements  to  the  current  implementation  of  the 
model . 

As  in  the  case  of  the  horizontal  array,  it  is  of  interest 


to  determine  if  the  spatial  covariance  changes  with  time.  To  examine 
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this  dependence  upon  time,  the  covariance  between  channel  1  and  all  the 
other  channels  in  the  vertical  array  was  calculated  at  x=0  for  the  entire 
10  ms  interval.  Although  this  does  not  provide  an  examination  of  all  of 
the  possible  covariances  of  the  vertical  array,  nevertheless  it  should 
demonstrate  the  significant  aspects  of  the  dependence  of  the  vertical 
covariance  upon  time.  The  envelopes  of  the  experimental  covariance 
are  shown  in  Figs.  V-37  and  V-38.  In  like  fashion  to  the  horizontal 
envelopes  of  the  covariance,  the  vertical  envelope  for  each  separation 
does  not  change  significantly  with  time.  The  envelope  of  the  normal¬ 
ized  auto-covariance,  "e  1(t,0),  is  1  by  definition.  As  the  separation 
from  the  first  element  increases,  the  envelope  of  the  covariance 
gradually  decreases.  The  fluctuations  in  the  envelopes  can  be 
attributed  to  the  limited  ensemble  size.  Thus  it  is  expected  that  the 
results  of  Fig.  V-33,  which  show  the  dependence  of  the  envelopes  on 
vertical  spatial  separation  at  one  time  t=75.25  ms,  will  be  closely 
duplicated  at  all  times  within  the  10  ms  interval  examined. 

The  phases  of  the  covariance  for  the  vertical  array 
are  shown  in  Figs.  V-39  and  V-1+0.  Unlike  the  phases  of  the  horizontal 
array,  the  vertical  phases  change  with  time  even  when  the  envelope 
values  are  not  small.  The  phases  change  linearly  with  time,  and  the 
rate  of  change  increases  as  the  separation  increases.  This  effect  can 
be  explained  qualitatively  by  recalling  that  the  covariance  at  increas¬ 
ing  times  is  due  to  scattering  from  the  surface  at  increasing  ranges. 

The  difference  in  path  lengths  from  the  scattering  region  to  two  vertical 
receivers  changes  with  range,  thus  resulting  in  a  change  in  the  phase 
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THE  DEPENDENCE  OF  THE  EXPERIMENTAL  ENVELOPES 
OF  THE  NORMALIZED  COVARIANCE  ON  TIME  AT  r  =  0 
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FIGURE  V-40 

THE  DEPENDENCE  OF  THE  EXPERIMENTAL  PHASES  OF 
THE  NORMALIZED  COVARIANCE  ON  TIME  AT  r  =  0 
FOR  THE  VERTICAL  ARRAY 
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of  the  covariance  with  time.  Receivers  which  are  farther  apart  have 
path  length  differences  to  the  scattering  region  which  change  more 
rapidly  as  the  scattering  region  moves  out  in  range  (time);  thus  the 
phase  of  the  covariance  changes  more  rapidly  with  time  for  these  re¬ 
ceivers.  These  effects  are  not  seen  in  the  horizontal  phases  since 
the  path  length  differences  from  the  scattering  region  to  two  horizontal 


receivers  does  not  change  appreciably  with  range. 

The  significance  of  the  change  in  phase  with  time  is 
that  even  though  the  envelope  of  the  covariance  does  not  change  with 
time,  the  covariance  changes  with  time  in  an  oscillatory  fashion.  For 
t=0,  the  covariance  is  given  by 

k.  At,  0)  =  e.  (t,0)cosU.  ,(t,'0)) 

1  »  J  1  5  J  1  *  J 

Since  the  phase  changes  linearly  with  time,  this  can  be  written  as 

k.  ,(t,0)  =  %.  (t,0)eos(iD.  t) 

1>J  1  >  J 

where  the  term  w.  represents  the  rate  of  change  of  the  phase  with 
i  *  J 

time,  and  is  dependent  upon  the  separation  of  the  two  receivers  (as 
well  as  other  parameters).  Thus  the  vertical  covariance  is  not  constant 
with  time  or  equivalently  range,  but  instead  oscillates  at  a  frequency 
which  depends  upon  various  geometrical  parameters. 

Theoretical  estimates  of  the  vertical  phase  were  made 
at  1  ms  intervals  and  are  compared  to  the  experimental  estimates  in 
Figs.  V-Ul  and  V-l+2.  The  theoretical  estimates  also  indicate  a 
linear  change  in  the  phase  of  the  covariance  with  time,  and  in  fact 
indicate  the  same  rate  of  change  as  was  measured  experimentally.  How¬ 


ever,  the  theoretical  phase  has  a  constant  offset  from  the  experimental 


FIGURE  V-41 
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phase  by  an  amount  which  depends  upon  the  separation. 


One  possible  explanation  for  this  offset  can  be  found 
in  the  orientation  of  the  vertical  array  with  the  surface.  Any  change 
that  the  angle  of  the  normal  to  the  planar  face  of  the  vertical  array 
makes  with  the  surface  will  change  the  path  length  difference  between  the 
scattering  region  on  the  surface  and  two  vertical  receivers,  thus 
changing  the  phase  of  the  covariance.  The  theoretical  estimates  shown 
in  Figs.  V-Ul  and  V-U2  were  made  by  considering  that  the  planar  faces 
of  the  vertical  array  and  projector  were  parallel,  and  that  their 
normals  intersected  the  plane  of  the  surface  at  a  10.5  degree  grazing 
angle.  Theoretical  estimates  were  also  made  with  the  normal  to  the  face 
of  the  vertical  array  intersecting  the  surface  at  a  grazing  angle  of 
9.8  degrees,  while  the  projector  maintained  its  10. 5  degree  grazing 
angle.  The  results  are  shown  in  Figs.  V-h3  and  V-kh.  The  offset  has 
been  removed  and  the  theoretical  phases  are  in  much  better  agreement 
with  the  experimental  results.  Thus  only  a  0.7  degree  grazing  angle 
difference  between  the  projector  and  vertical  receiving  array  can  make 
a  significant  change  in  the  prediction  of  the  phase  of  the  covariance, 
and  thus  the  covariance  itself.  It  must  be  pointed  out  that  it  cannot 
be  said  conclusively  that  this  is  the  reason  for  the  offset  between 
the  experimental  and  theoretical  results.  On  the  other  hand,  the 
mechanical  mounting  of  the  projector  and  vertical  array  in  the  experi¬ 
ment  was  such  that  their  angular  orientations  could  certainly  not  be 
guaranteed  to  be  the  same  to  an  accuracy  of  less  than  one  degree. 

Thus  it  is  feasible  that  this  angular  difference  produced  the  offset  in 
the  theoretical  and  experimental  phases. 
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V.3.1* 


The  Covariance  for  t^O 


Thus  far,  the  variance  as  a  function  of  t,  the  auto¬ 
covariance  as  a  function  of  t  and  t,  and  the  covariance  as  a.  function 
of  spatial  separation  and  t  for  t=0  have  been  examined  both  experimen¬ 
tally  and  theoretically  for  the  horizontal  and  vertical  arrays.  Thi 
section  will  now  present  the  covariance  as  a  function  of  spatial  separ¬ 
ation  and  t.  The  horizontal  covariance  is  presented  first,  followed 
by  the  vertical  covariance. 

Figure  V-U5  presents  the  normalized  envelopes  of  the 
sample  covariance  between  the  various  elements  of  the  horizontal  array 
as  a  function  of  t.  The  diagonal  plots  in  this  "matrix"  of  envelopes 
are  the  auto-covariances  shown  previously.  The  first  row  of  plots 
represents  the  envelopes  of  the  covariance  between  channel  11  and  the 

other  horizontal  receiver  channels,  i.e.,  the  first  plot  is  e  ',t,x- 

J.  ,  1^. 

the  next  one  to  the  right  is  et  _,(t,x),  etc.  Likewise  the  second  r:w 

11  ,  - 

contains  the  envelopes  between  channel  12  and  the  others;  the  first  pio 

is  'e’^2  ^(t,x),  then  e^ 2  ,0(t,x),  and  so  forth.  Because  the  horizontal 

elements  were  -unequally  spaced,  adjacent  plots  do  not  represent  equal 

changes  in  separation.  Thus  e'! .  .  ?(t,x)  is  significantly  less  than 

. ( t , i )  because  the  separation  between  the  elements  changed  from 
-L 

zero  for  e  ..  11(t,x)  (the  auto-covariance )  to  5.25  cm  for  e. ,  .  „(t,:', 
whereas  'e'.  ,  , ,  (t,x)  has  not  decreased  as  much  compared  to  ,  .(t,-’1, 
even  though  they  are  also  adjacent  plots  in  the  matrix,  because  the 
separation  for  et „  , (t,x)  is  only  2.5k  cm. 

-t>  1  1  ^ 

Several  observations  can  be  made  about  the  envelopes 
r:  the  ?o variance .  In  general,  e .  ( t,i;  =  e,  .vt,c  on^y  tor 

i  n1  v'  *  -■ 
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FIGURE  V-45 

NORMALIZED  ENVELOPES  OF  SAMPLE  ESTIMATES  OF 
THE  COVARIANCE  OF  THE  HORIZONTAL  RECEIVERS 
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Thus  symmetry  about  the  diagonal  of  the  matrix  is  not  required.  In 


this  instance,  the  envelopes  are  reasonably  symmetric  about  the  diagonal 
plots,  but  are  not  exactly  symmetric.  If  two  processes  are  jointly 
wide-sense  stationary,  then  it  can  be  shown  that 


ki,j (,)  ' 


This  symmetry  can  be  realized  when 


and 


Thus  some  of  the  symmetry  of  the  envelopes  observed  would  not  be  un¬ 
expected  if  the  reverberation  were  approximately  jointly  wide-sense 
stationary. 


It  can  also  be  observed  that  the  dependence  of  the 
envelope  of  the  covariance  on  r  is  approximately  the  same  for  all 
pairs  of  receivers.  Thus  the  time  extent  of  the  cross-covariance  is 
determined  by  the  duration  of  the  pulse  after  it  has  been  filtered  by 
the  transmit  and  receive  apertures,  as  was  previously  shown  for  the 
auto-covariance . 


The  phases  of  the  sample  covariance  are  shown  in 
Fig.  V-16  for  the  various  elements  of  the  horizontal  array.  The  phase 
is  not  constant  over  the  range  of  x  shown,  indicating  that  the  phase 
will  also  make  a  contribution  to  changing  the  covariance  as  a  function 
of  x.  The  phases  of  the  auto-covariances  are  zero  at  x=0.  and 
t.  ,(t,0)  =  -i> ,  ,(t,0),  as  required  for  the  difference  corr.ror.er.x  . 

t  1  o  v.  5  1 

Theoretical  estimates  of  the  envelope  and  phase  of  the  normal¬ 
ized  covariance  for  the  horizontal  array  were  calculated  for  the  three 


smallest  separations.  Figures  V-U7  and  V-U8  show  the  results  overlaid 
on  the  sample  estimates  of  the  envelope  and  phase.  The  theoretical 
envelopes  show  very  good  agreement  with  the  sample  estimates.  The 
model  is  correctly  predicting  the  shape  of  the  envelopes  as  well  as 
the  magnitude. 

The  theoretical  phases,  shown  in  Fig.  V-hB,  indicate 
good  agreement  between  closely  spaced  elements,  but  poor  agreement  for 
elements  spaced  farther  apart.  However,  since  these  elements  have 
small  envelopes,  the  phase  will  not  make  the  covariance  differ  signif¬ 
icantly  from  zero.  It  is  interesting  to  note  that  the  theoretical 
phases  exhibit  the  symmetry  associated  with  wide-sense  stationar ity . 

It  has  already  been  shown  that  the  covariance  cf  the 
vertical  array  is  much  greater  as  a  function  of  separation  than  the 
horizontal  array.  The  envelopes  of  the  normalized  covariance  between 
all  the  possible  pairs  of  nine  elements  in  the  vertical  array  are  shew 
in  Figs.  V-U9  through  V-52.  The  greater  covariance  ' s  also  obvious 
in  these  plots.  It  can  be  seen  that  all  the  envelopes  have  a  similar 
dependence  upon  t.  That  is,  the  value  of  t  for  which  the  envelopes 
decrease  to  zero  is  the  same  for  all  the  envelopes. 

The  corresponding  phases  are  shown  in  Firs.  Y-55  *hr~. 
7-56.  It  was  shown  earlier  that  the  phase  at  any  specified  separation 
changed  with  time.  In  these  figures  the  change  in  phase  with  serara* i 
at  a  specified  time  is  more  obvious.  The  change  in  chase  as  a  fur.ot  i  * 

cf  t  is  also  evident.  The  phases  dc  not  exhibit  the  symmo*  ry  icr.on- 
strated  by  the  horizontal  phases.  Furthermore,  the  pnase?  ohanre  w ; • 
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FIGURE  V-47 

COMPARISONS  OF  NORMALIZED  ENVELOPES  OF  SAMPLE  ESTIMATES  AND 
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FIGURE  V-50 

ENVELOPES  OF  THE  SAMPLE  ESTIMATE  OF  THE  NORMALIZED  COVARIANCE 
OF  THE  VERTICAL  ARRAY  AS  A  FUNCTION  OF  r  FOR  t  -  75.25  m*«c 
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time.  Thus  the  vertical  covariance  does  not  represent  a  jointly  wi de¬ 
sense  stationary  process. 

Since  the  model  did  not,  correctly  predict  the  dependence 
of  the  envelope  of  the  vertical  covariance  on  spatial  separation,  com¬ 
putation  of  the  theoretical  covariance  for  all  pairs  of  elements  of  the 
vertical  array  would  not  provide  much  additional  information.  However, 
the  covariance  between  channel  1  and  the  other  vertical  elements  was 

computed  and  is  displayed  in  Fig.  V-57 •  Even  though  the  model  predicted 
the  temporal  extent  of  the  envelope  correctly,  the  magnitude  of  the 
envelope  was  overestimated,  as  was  shown  previously.  The  phases  were 
computed  with  the  vertical  array  and  projector  at  the  same  grazing 
angle  with  the  surface.  Thus,  although  the  theoretical  phase  changes 
with  separation  as  does  the  sample  estimate  of  the  phase,  there  i s  an 
offset  between  the  two  estimates. 

V . 3 . 5  Examination  of  the  Failure  of  the  Model 

In  order  to  gain  insights  into  why  the  model  failed  to 
correctly  account  for  the  dependence  of  the  envelope  of  the  covariance 
on  vertical  separation,  it  is  important  to  consider  that  the  model 
provided  good  agreement  with  experiment  in  many  instances.  The  model 
correctly  predicted  the  dependence  of  the  variance  on  time,  indicating 
the  grazing  angle,  directivity  patterns,  and  locations  of  sensors  were 
being  correctly  taken  into  account.  It  also  predicted  the  temporal 
extent  of  the  covariance  by  taking  into  account  the  pulse  length  and 
frequency  response  of  the  projector  and  receivers.  Most  importantly, 
the  dependence  of  the  envelope  of  the  covariance  on  spatial  separation 
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COMPARISON  OF  THE  THEORETICAL  AND  SAMPLE  ESTIMATES  OF  THE  ENVELOPES 
AND  PHASES  OF  THE  NORMALIZED  COVARIANCE  FOR  CHANNEL  1  OF  THE 
VERTICAL  ARRAY  AS  A  FUNCTION  OF  t  FOR  t -  75.26  m»*c 
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in  the  horizontal  orientation  was  correctly  taken  into  account.  Any 
possible  explanations  of  the  failure  of  the  model  to  predict  the  vertical 
covariance,  and  subsequently  any  changes  to  the  model,  must  allow  the 
model  to  continue  correctly  predicting  these  things.  With  this  guideline 
in  mind,  it  is  appropriate  to  examine  the  model  to  determine  why  it 
correctly  predicted  the  horizontal  covariance. 

First  of  all,  it  is  important  to  determine  why  the  model 

correctly  predicted  that  the  horizontal  covariance  would  decrease  rapidly 
with  spatial  separation.  Determination  of  the  predominant  cause  of  this 
decrease  in  covariance  may  indicate  what  additional  factors  should  be 
considered  which  would  decrease  the  covariance  in  the  vertical  case.  The 
basic  assumption  of  the  model  is  that  the  backscatt.ered  return  at  a 
receiver  is  the  result  of  a  linear  superposition  of  returns  from  point 
scatterers  imbedded  in  the  surface: 

(t,Aj)  , 


where  Uj(t,Aj)  is  the  return  from  the  ,1th  scatterer  at  position  A^, 
M(t')  is  the  number  of  contributing  scatterers  at  time  t and  X(t) 
is  the  total  return  at  the  receiver  at  time  t.  The  covariance  between 
two  spatially  separated  receivers  is  due  to  scatterers  which  are 
providing  returns  to  both  receivers.  If  an  elementary  return  is 
different  at  the  two  receivers,  and  if  all  the  elementary  returns  are 
contributing  differently,  then  it  could  be  expected  that  when  all  the 
elementary  returns  are  combined,  the  returns  at  the  two  receivers  cou] d 
be  incoherent.  More  specifically,  the  covariance  between  two  returns 


X1(t1)  and  X2(t2)  has  been  shown  in  Chapter  III  to  be  given  by 
K(l2)ttrt2  |t')  =  /  p12U,t')(ui(t1,X,01) 

•  U2(t2,A,02)^)edX  , 

which  in  the  present  study  was  simplified  through  the  appropriation  of 


several  assumptions  to 


'(t1,t2|t')  =  Os  /  U1(t1,X)U2(t2,A)dS 


where  the  integration  is  over  the  plane  of  the  surface  which  contains 

scatterers  distributed  uniformly  in  the  region  A^2  with  a  density  os. 

l^i ( ti t A )  is  the  elementary  return  at  receiver  1  at  time  t^  due  to  the 

scatterer  at  position  A,  and  U2 C 1 2 , A )  is  the  return  at  receiver  2  at 

time  t2  due  to  the  same  scatterer.  If  the  integration  of  the  product 

of  and  U 2  over  the  scattering  region  is  small,  then  the  covariance 

will  be  small.  This  is  the  case  if  the  product  of  and  U2  is  approxi 

mately  equally  distributed  about  zero  over  the  range  of  integration. 

It  was  shown  in  Chapter  III,  Section  2,  that  this  product  depended  upon 

/RR2-RR11  /RR?-RR1  1 

terras  containing  cosgjI - - - 1  and  sinw  I — — - 1  where  RR^  is  the 


path  length  from  scatterer  to  receiver  1,  and  R_„  is  the  rath  length 

from  the  same  scatterer  to  receiver  2.  Thus  if  this  path  length  differ 

^o  °1T 

ence  R^.-R-,..  changed  significantly  with  respect  to  —  =  —  over  the 

rid  n±  1  c  a 

various  contributing  scatterers,  then  the  product  of  and  would 

-L  <- 

oscillate  and  integrate  to  approximately  zero.  Therefore  it  can  be 
concluded  that  in  a  very  qualitative  sense,  if  the  path  length  differ¬ 
ence  between  a  scatterer  and  the  two  receivers  changes  significantly 


2o8 

over  the  region  of  contributing  scatterers,  then  the  covariance  between 
those  two  receivers  will  be  small.  Conversely,  if  the  path  length 
difference  changes  little,  then  the  covariance  will  be  large. 

Figure  V-58  illustrates  the  path  lengths  between  two 
receivers  separated  vertically  and  horizontally  for  conditions  similar 
to  the  present  study.  The  horizontal  receivers  will  be  considt  *. . 

The  area  that  the  projector  illuminates  on  the  surface  at  any  om  •  i:ne 
has  a  small  downrange  extent  determined  by  the  pulse  length,  and  a  larger 
crossrange,  or  azimuthal,  extent  determined  by  the  horizontal  beamwidth 
of  the  projector.  It  can  be  observed  that  the  path  length  difference 
RR2-RR1  chang«s  over  ^e  region  of  scatterers.  For  a  scatterer  left 
of  center  is  greater  than  R  ,  hut  for  a  scatterer  right  of  center 
the  opposite  is  true.  Thus  the  geometry  allows  for  a  small  covariance, 
depending  upon  how  much  the  path  length  differences  change  over  the 
scattering  region.  This  in  turn  is  determined  by  the  azimuthal  width 
of  the  scattering  region  and  the  horizontal  separation  of  the  receivers. 

For  the  present  study,  most  of  the  energy  from  the  projector  was  confined 
horizontally  to  angles  ±15  degrees  from  center.  An  approximate  calcu¬ 
lation  of  the  change  in  path  length  difference  over  the  30  degree  azimuthal 
width  of  the  scattering  region  for  the  ranges  and  geometries  used  was 
made.  For  two  elements  separated  horizontally  by  2.5l  cm,  the  closest 
horizontal  spacing  used,  the  path  length  differences  changed  by  .7  A, 
where  A  is  the  wavelength  of  the  transmitted  sound.  This  will  result 
in  some  oscillation  of  the  integrand  of  the  covariance  integral  and 
result  in  a  decreased,  although  probably  not  zero,  covariance.  This  is 


(a)  VIEW  FROM  ABOVE  OF  HORIZONTAL  RECEIVERS  AND  PROJECTOR 


(b)  VIEW  FROM  THE  SIDE  OF  VERTICAL  RECEIVERS  AND  PROJECTOR 
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what  was  observed  theoretically  and  experimentally.  For  a  horizontal 
separation  of  15  cm,  the  largest  horizontal  spacing  used,  the  change 
in  path  length  difference  was  U  X,  which  would  produce  significant 
oscillation  of  the  integrand  term  and  result  in  very  low  covariance. 
Again,  this  is  what  was  observed.  This  agreement  indicates  that  the 
qualitative  explanation  of  the  covariance  is  a  reasonable  one.  This 
could  also  explain  the  larger  covariance  measured  by  Frazer,  since  his 
projector  was  confined  horizontally  to  angles  of  only  +10  degrees, 
resulting  in  less  change  of  the  path  length  difference.  It  was  also 
found  that  the  path  length  difference  changed  negligibly  with  the  small 
downrange  extent  of  the  scattering  region. 

The  path  lengths  to  the  vertical  receivers  and  the 
scattering  region  are  indicated  in  Fig.  V-58  from  a  side  view.  The 
path  lengths  are  different,  but  it  can  be  seen  that  the  difference 
changes  little  over  the  region  of  scatterers.  Again,  an  approximate 
calculation  indicated  that  over  the  30  degree  azimuthal  extent,  the 
path  length  difference  for  a  10  cm  vertical  separation  changed  'inly 
.03X.  For  a  vertical  separation  of  1  m,  the  change  in  path  length 
difference  was  only  . 3X .  Thus  it  can  be  expected  from  this  qualitative 
argument  that  the  model  would  predict  a  very  large  covariance  even  for 
large  vertical  separations. 

Since  this  reasoning  has  qualitatively  predicted  the 
performance  which  has  been  obtained  by  the  model,  then  the  same  reasoning 
may  indicate  possible  changes  to  the  model  which  would  bring  the  vertical 
covariance  into  better  agreement  with  the  experimental  results.  In 
particular,  it  can  be  observed  that  in  order  to  change  the  path  length 


differences  for  vertical  receivers  in  the  same  way  in  which  the  horizon¬ 
tal  path  length  differences  were  changed,  it  would  be  necessary  to 
consider  the  scatterers  distributed  vertically  below  the  surface  as 
well  as  on  the  surface  and  integrate  over  the  depth  of  the  scatterers 
as  well  as  on  the  surface.  Thus  the  covariance  could  be  calculated  as 

K(l2)(ti,t2|tO  =  /A  PyUj/t^AjUgU^AteV 

where  the  integration  is  now  over  a  volume  of  scatterers.  The  scatter 
density  is  now  a  volume  scatter  density  and  is  a  function  of  the  depth 
below  the  surface.  This  implementation  would  be  appropriate  if  scatter¬ 
ing  could  be  physically  considered  to  occur  throughout  a  layer  bounded 
by  a  flat  surface.  Although  the  implementation  of  the  model  can  pres¬ 
ently  only  perform  a  surface  integration,  it  was  run  at  several  differ¬ 
ent  depths  and  the  results  manually  integrated  over  the  depth.  It 
was  found  that  the  vertical  covariance  could  be  reasonably  reduced, 
depending  upon  the  particular  dependence  of  the  density  upon  depth  and 
the  depth  of  the  scattering  layer.  For  example,  by  considering  a  con¬ 
stant  scatter  density  as  a  function  of  depth,  the  value  of  the  envelope 
of  the  covariance  at  t=0  between  receivers  1  and  2  could  be  reduced 
from  .99  to  .91  by  integrating  over  a  1.2  m  depth,  and  could  be  further 
reduced  to  as  little  as  .2k  by  an  integration  over  a  depth  of  3.2  m. 

The  measured  value  of  the  envelope  was  .8k. 

Although  this  implementation  is  conceptually  straight¬ 
forward,  it  does  have  some  practical  problems.  Fpecifically ,  it  ir 
not  obvious  how  the  depth  of  the  scattering  layer  and  the  density  of 
scatterers  as  a  function  of  depth  would  be  determined.  Secondly,  it  is 
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not  known  if  scattering  actually  occurs  in  a  layer  below  the  surface 
or  if  it  occurs  Just  at  the  surface.  Thus  another  implementation  could 
consider  scattering  at  the  surface  only  but  allow  the  surface  to  vary 
randomly  in  height.  The  original  expression  for  the  covariance  contained 
an  expected  value  of  the  product  of  and  with  respect  to  random 
parameters  0.  By  retaining  one  random  parameter,  the  surface  wave 
height  z,  and  then  performing  a  surface  integral  over  the  mean  surface 
level,  the  covariance  could  be  expressed  as 
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where  A  is  the  mean  location  of  the  surface  and  <>z  denotes  the  expected 
value  with  respect  to  the  random  surface  waveheight  z.  This  implemen¬ 
tation  has  the  advantage  of  the  direct  inclusion  of  the  statistical 
distribution  of  the  surface  wave  heights,  which  is  potentially  measurable. 

It  is  anticipated  that  neither  of  these  implementations 
would  change  the  results  for  the  horizontal  covariance,  since  the  path 
length  differences  to  horizontal  elements  do  not  change  significantly 
with  depth.  It  is  also  anticipated  that  the  temporal  results  would 
be  unaffected.  Thus  it  appears  that  these  implementations  are  reasonable 
ones  to  correct  for  the  vertical  covariance. 

Finally,  it  is  encouraging  to  observe  that  the  basic 
assumption  of  a  superposition  of  point  scatterers  can  accurately  model 
many  aspects  of  the  covariance  of  backscattering,  and  can  at  least 
qualitatively  explain  those  areas  which  it  cannot  presently  model  well . 
Thus  it  is  felt  that  the  basic  model  itself  is  still  a  reasonable  model 


and  that  its  limitations  are  due  to  the  simplifications  which  were 


made  for  the  current  implementation  of  the  model.  It  is  felt  that  the 
present  study  has  further  demonstrated  the  validity  and  usefulness  of 
the  model,  and  has  provided  a  direction  for  further  improvements  to  its 
implementation . 


VI.  SUMMARY  AND  CONCLUSIONS 


VI .1  Objectives  of  the  Study 

The  statistical  properties  of  surface  reverberation  depend  upon 
many  varied  environmental,  geometrical,  and  sonar  parameters.  A  complete 
understanding  of  surface  reverberation  requires  its  study  under  all 
these  conditions.  Within  this  broader  context,  the  present  study  has 
attempted  to  examine  a  few  of  the  statistical  properties  under  a  limited 
set  of  these  parameters.  Four  objectives  were  identified  for  this  study 

1.  The  primary  objective  was  to  simultaneously  measure  the 
dependence  of  the  covariance  upon  spatial  separation  of  the 
observation  points  in  both  the  horizontal  and  vertical 
orientations.  This  also  involved  determining  the  change  of 
the  spatial  covariance  as  a  function  of  time  of  observation 
(t)  and  the  difference  in  the  observation  times  (x).  Al¬ 
though  the  most  closely  examined  parameter  was  spatial 
separation,  it  was  hoped  that  the  dependence  upon  some  other 
parameters  could  be  inferred  as  well. 

2.  The  second  objective  was  to  compare  the  measured  results  to 
a  theoretical  model  of  reverberation.  In  this  way  it  was 
hoped  that  our  physical  understanding  of  the  acoustic 
scattering  phenomena  could  be  extended.  It  was  of  interest 
to  examine  a  model  which  included  most  of  the  required 
parameters  in  a  general  way,  which  allowed  the  computation 
of  various  statistical  parameters,  and  which  had  not  been 


extensively  validated  experimentally.  Thus  Mi  1  ii don 
point-scatter  model  was  chosen  as  the  particular  mod- 

examine  and  compare  to  the  measured  results.  It  war, 
hoped  that  if  the  experimental  measurements  illustrat 
deficiencies  in  the  model,  extensions  to  the  model  c„, 
identified  which  would  more  accurately  predict  surfae 
reverberation . 

3.  A  third  objective  of  the  study  was  to  experimentally 
the  first  four  univariate  moments  and  perform  univari 
tests  for  normality.  These  measurements  would  rrovii 
indication  of  the  type  of  distribution  of  the  reverse 
Although  no  density  estimates  were  performed,  it  woul 
possible  to  do  so  with  the  data  presented  in  this  stu 

k.  Finally,  the  fourth  objective  of  the  study  was  to  ref 
the  techniques  for  the  measurement  and  validation 
beration  from  multiple  receivers.  A  stable  platform 
desired  to  minimise  possible  effects  due  platform 
Larger  ensemble  sizes  were  desired  ir.  order  t:  reduce 
sample  error.  Statistical  *  ?sts  to  efficiently  val:  : 
these  larger  ensembles  were  also  r.eeied. 


Distinctives  of  the  Stud} 


VI. 2.1  Experimental  Distinctives 


The  measured  sound  was  limited  to  backs 


sound  from  the 


wind-ro 


3.  The  receivers  were  linear  arrays  oriented  horizon¬ 
tally  and  vertically. 

U.  A  stable  platform  consisting  of  a  bottom-mounted 
tower  was  used. 

5.  Ensemble  estimates  were  made  of  the  various  statis¬ 
tical  properties,  as  opposed  to  time-average  esti¬ 
mates  . 

6.  The  possible  sources  of  measurement  error  were 
identified  and  reduced  to  an  acceptable  level 
through  the  use  of  specialized  data  collection/ 
recording/digitization/processing  techniques . 

7-  Extensive  testing  of  the  data  was  performed  to 
validate  the  ensembles  and  to  test  for  normality. 


8.  The  first  four  moments  and  the  covariance  were 


measured. 


VI. 2. 2  Theoretical  Distinctives 


The  model  assumed  a  random  surface  composed  of  a 
collection  of  independent  point  scatterers  in  a 
homogeneous  medium. 

2.  The  model  allowed  for  general  geometries,  including 
the  spatial  distribution  of  the  sensors  and  the 
grazing  angle. 

3.  General  sources  and  receivers  were  allowed:  fre¬ 
quency,  pulse  length,  directionality,  bandwidth. 


signal  spectrum,  and  aperture  response  could  be  taken 
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into  account. 

U.  The  model  was  implemented  for  bandlimited  transmis¬ 
sions  . 

5.  The  scatterers  were  modeled  as  perfect  point  reflec¬ 
tors  distributed  uniformly  over  the  surface. 

VI . 3  Review  of  the  Study 
VI. 3.1  Chapter  I 

Chapter  I  presented  an  introduction  to  the  study, 
citing  several  similar  studies  and  exposing  the  need  for  the  present 
study.  It  was  shown  that  the  coherence  of  forward-scattered  sound 
has  been  examined  for  both  vertical  and  horizontal  arrays,  but  that  it 
was  expected  that  backscattered  sound  would  have  significantly  differ¬ 
ent  statistical  properties.  It  was  also  shown  that  explosive  sources 
have  been  used  to  measure  the  coherence  of  backscattered  sound  for 
both  horizontal  and  vertical  arrays  using  time  averages.  Significant 
differences  were  determined  between  the  coherence  of  horizontal  and 
vertical  arrays;  thus  it  was  expected  that  similar  differences  would 
be  observed  using  a  pulsed  source  and  ensemble  averages. 

VI. 3.2  Chapter  II 

The  second  chapter  described  the  experimental  measure¬ 
ments.  Estimates  of  the  effects  of  reverberation  to  ambient  noise  ratio, 
crossover  between  receiver  channels,  and  finite  sample  size  on  the 
accuracy  of  the  calculation  of  the  coherence  were  made.  It  was  shown 
that  with  the  measurement  system  used,  the  primary  limitation  was  finite 
sample  size.  Since  the  accurate  measurement  of  coherence  between  re- 
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ceivers  also  requires  the  preservation  of  the  phases  of  the  reverberation 
returns  at  the  face  of  the  receivers,  the  various  sources  of  phase 
error  which  could  be  introduced  into  the  system  were  identified.  The 
techniques  for  minimizing  or  correcting  these  phase  errors  were  described. 
The  techniques  required  the  proper  interaction  of  special  data  collection, 
analog  recording,  digitization,  and  digital  processing. 

One  projector  and  two  line  receiving  arrays  were  used 
to  collect  the  data.  The  vertical  array  consisted  of  nine  elements 
with  a  maximum  separation  of  slightly  less  than  1  m.  The  horizontal 
array  contained  four  elements  with  approximately  a  .15  m  maximum  separ¬ 
ation.  A  100  ps  pulsed  CW  signal  at  a  center  frequency  of  80.0  kHz 
was  transmitted  at  the  surface  at  a  10.5  degree  grazing  angle.  The 
center  of  the  array  was  10.2  m  below  the  surface  in  approximately  Uo  m 
of  water.  The  water  was  isothermal  and  wave  heights  in  the  35-50  kph 
wind  were  approximately  .3-. 6  m. 

VI. 3. 3  Chapter  III 

The  theoretical  model  was  described  in  the  third  chapter 
of  this  study.  Reverberation  was  modeled  as  a  linear  superposition  of 
the  scattered  sound  from  independent  point  scatterers  representing  in¬ 
homogeneities  at  the  surface  of  a  homogeneous  medium.  Special  attention 
was  given  to  the  assumptions  inherent  in  the  model,  and  the  additional 
assumptions  made  for  its  current  implementation.  The  general  model 
was  described  first,  and  an  expression  for  the  covariance  function 
was  given.  Several  simplifying  assumptions  were  then  introduced.  Among 
others,  the  transmitted  signal  was  assumed  to  be  bandlimited,  and  the 
scatterer  was  modeled  as  a  perfect  point  reflector.  An  expression 


for  the  elementary  scattered  waveform,  representing  the  waveform  from 
a  single  scatterer,  was  then  developed,  along  with  a  simplified  ex¬ 
pression  for  the  covariance.  A  quadrature  representation  of  the  co- 
variance  was  then  developed  in  order  to  allow  the  presentation  of  the 
covariance  in  terms  of  its  envelope  and  phase.  The  specifics  of  imple¬ 
menting  the  model  on  a  computer  were  then  discussed.  Some  of  the 
details  of  the  development  of  the  model  were  expanded  upon  in  the  three 
appendices. 

VI .  3 . 1*  Chapter  IV 

Chapter  IV  began  the  analysis  of  the  data.  The  pro¬ 
cedure  for  sampling  the  data  to  form  ensembles  was  described,  along 
with  the  assumptions  inherent  in  this  sampling  process.  A  10  ms  rever¬ 
beration  return  was  sampled  at  320  kHz,  resulting  in  3200  samples  in 
time  for  each  of  108U  pings  and  13  receivers.  Ensembles  were  formed 
at  each  sample  time,  resulting  in  3200  ensembles  of  108U  pings  each. 
Subsequent  analysis  of  the  data  was  restricted  to  the  last  500  samples 
of  each  ensemble  because  of  a  changing  reverberation  intensity  in  the 
first  part  of  the  ensembles. 

All  3200  ensembles  of  each  channel  were  first  tested  fo 
statistical  validity  and  then  for  normality.  A  general  overview  of 
statistical  hypothesis  testing  was  given,  and  a  description  of  .w 
test  which  was  employed  was  provided.  Since  each  sample  ensenb !•'  : 
treated  as  a  random  sampling  of  the  reverberation  process,  then  -a or. 
ensemble  should  consist  of  random  variables  which  are  independer.*  on  i 
identically  distributed  (i.i.d.).  To  verify  the  i.i.i.  properties  of 
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the  ensembles,  each  ensemble  was  tested  for  randomness  (independence) 
and  homogeneity  (identical  distribution).  Independence  guarantees  that 
the  generation  of  one  random  variable  in  the  ensemble  was  not  affected 
by  the  generation  of  any  others,  while  homogeneity  implies  that  all 
the  random  variables  come  from  the  same  parent  population,  that  is, 
they  all  have  the  same  distribution.  The  runs  up  and  down  test  was 
used  for  testing  for  randomness,  while  the  Kolmogorov-Smirnov  two- 
sample  test  and  the  Wilcoxon  rank-sum  test  were  used  for  testing  for 
homogeneity.  Results  of  the  tests  were  presented  for  all  channels. 
After  determining  the  validity  of  the  ensembles ,  they  were  then  tested 
to  determine  if  their  underlying  distribution  was  gaussian.  Four  tests 
for  normality  were  applied:  Pearson's  test  of  skewness,  Pearson's  test 
of  kurtosis,  D'Agostino's  test,  and  the  Kolmogorov-Smirnov  one-sample 
test . 

VI. 3. 5  Chapter  V 

The  main  objective  of  the  study  was  finally  addressed 
in  Chapter  V.  First  the  techniques  for  the  numerical  computation  of 
the  various  moments  and  the  covariance  were  given.  It  was  shown  that 
the  covariance  from  a  bandlimited  transmission  consists  of  a  time-sum 
and  a  time-difference  component.  The  validity  of  ignoring  the  time-sum 
component  was  briefly  discussed.  The  time-difference  component  of 
the  covariance  was  expressed  in  terms  of  an  envelope  function  and  a 
phase  function.  The  technique  for  normalizing  the  covariance  was 
defined.  The  sample  mean,  variance,  skew,  and  kurtosis  were  presented 
for  a  limited  section  of  the  10  ms  reverberation  return  for  all  13 
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channels.  The  difference  component  of  the  variance  was  also  presented 


and  compared  to  the  actual  variance. 

The  last  part  of  Chapter  V  provided  a  rather  thorough 
analysis  of  the  covariance.  All  the  analysis  was  performed  with  the 
difference  component.  The  results  from  the  theoretical  model  were  aiso 
presented  for  comparison  to  the  experimental  results.  First  the  variance 
was  presented  over  the  entire  10  ms  interval.  Then  the  normalized 
envelopes  and  phases  of  the  auto-covariance  were  presented  as  a  function 
of  the  time  difference  t  at  a  time  t=75.25  ms  after  transmit.  The 
effect  of  pulse  length  on  the  temporal  duration  of  the  covariance  was 
discussed,  and  stationarity  was  examined.  The  auto-covariance  was  then 

presented  over  the  entire  10  ms  interval,  both  normalized  and  unnormal¬ 

ized.  Secondly,  the  c.ross-co'i'ariance  for  t=0  was  also  given-  The 
envelope  of  the  covariance  as  a  function  of  spatial  separation  for  both 
the  horizontal  and  vertical  arrays  was  presented  for  one  time,  and  th-n; 

for  the  entire  10  ms  interval.  The  phase  of  the  covariance  at  t=0  was 

also  given  over  the  10  ms  interval.  For  the  vertical  array,  the  effect 
of  location  of  the  elements  on  the  array  was  examined  for  both  the 
envelopes  and  phases.  Thirdly,  the  covariance  envelopes  and  phases 
were  presented  for  t^O  at  one  time  as  a  function  of  spatial  separation. 
Finally,  a  qualitative  explanation  for  the  failure  of  the  model  to 
correctly  predict  the  dependence  of  the  vertical  covariance  on  spatial 
separation  was  presented.  Suggested  changes  to  the  model  were  also 
given. 

VI . U  Summary  of  the  Results  of  the  Study 

The  previous  section  contained  a  review  of  the  work  which  was 


performed  on  this  study.  The  results  of  the  data  analysis  will  now 


be  summarized.  The  results  will  be  presented  according  to  the  objec¬ 
tives  for  the  study. 

VI . k . 1  Computation  of  Vertical  and  Horizontal  Covariance 

1.  It  was  found  that  the  normalized  envelope  of  the 
horizontal  covariance  at  t=0  decreased  to  a  value  of 
0.1  at  a  separation  of  Just  over  four  wavelengths. 

A  comparison  was  made  to  the  horizontal  covariance 
measured  by  Frazer  under  similar  circumstances.  It 
was  found  that  the  present  study  obtained  smaller 
values  of  the  covariance.  The  difference  was 
attributed  to  the  illumination  of  more  "off-axis" 
scatterers  because  of  the  wider  transmit  horizontal 
beamwidth  used  in  the  present  study. 

2.  The  horizontal  covariance  at  t=0  was  also  determined 
as  a  function  of  time  after  transmit  over  the  10  ms 
interval  from  which  reverberation  samples  were  taken. 
It  was  found  that  the  normalized  envelopes  of  the 
horizontal  covariance  were  statistically  constant 
throughout  the  10  ms  interval;  any  fluctuations 
could  be  attributed  to  the  finite  ensemble  size. 

It  was  also  observed  that  the  phases  of  the  co- 
variance  were  constant  for  separations  which  had  signi 
icantly  non-zero  covariance  envelopes,  but  for 
small  values  of  the  covariance  envelopes  the  phases 


changed  with  time.  The  phases  which  were  constant 
were  non-zero,  indicating  a  non-uniform  scatter 
density  and  resulting  in  a  covariance  which  was  less 
than  the  envelope  value. 

The  dependence  of  the  vertical  covariance  on  spatial 
separation  was  also  examined.  The  vertical  co- 
variance  maintained  a  significant  level  at  much  larger 
separations  than  the  horizontal  covariance.  Even 
at  the  greatest  separation  of  a  little  over  50 
wavelengths,  the  normalized  envelope  of  the  vertical 
covariance  still  had  a  value  greater  than  0.2. 

It  was  also  determined  that  the  normalized  envelope 
of  the  covariance  between  vertically  separated 
receivers  did  not  change  significantly  over  the  10  ms 
interval.  However  the  phase  of  the  vertical  co- 
variance  did  change  with  time,  and  the  rate  of  change 
increased  with  increasing  vertical  separation. 

The  change  was  approximately  linear,  indicating  that 
the  actual  covariance  slowly  oscillated  with  time. 

The  vertical  array  also  contained  several  pairs  of 
elements  which  had  the  same  separation  at  different 
locations  on  the  array.  It  was  determined  that  the 
envelope  at  a  specified  separation  did  not  change 
significantly  with  location  on  the  array,  but  that 
the  phase  did.  Thus  the  covariance  was  dependent  not 
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only  on  the  separation  of  the  vertical  receivers, 
but  also  on  their  position. 

6.  The  horizontal  and  vertical  covariances  were  also 
examined  as  a  function  of  time  difference  (x).  It 
was  found  that  the  covariance  envelopes  decreased 
to  zero  in  a  time  difference  approximately  equal  to 
the  extent  of  the  transmitted  pulse  after  it  was 
filtered  by  the  projector  and  receiver.  The  100  us 
pulse  was  extended  in  time  to  approximately  275  us, 
and  the  envelope  of  the  covariance  decreased  to  zero 
in  approximately  250  us.  The  envelopes  exhibited  an 
approximate  even  symmetry  about  t=0,  and  the  phases 
had  approximately  an  odd  symmetry. 

7.  The  auto-covariance  for  one  receiver  was  also  shown 
as  a  function  of  t  ani  t  over  the  10  ms  interval. 
Non-stationarity  was  observed  in  the  unnormalized 
covariance  over  the  10  ms  due  to  the  change  of  the 
envelope  with  time.  However  much  of  this  non- 
stationarity  was  removed  from  the  auto-covariance  by 
normalization . 

8.  The  variance  (auto-covariance  at  t=0)  was  also  ex¬ 
amined  as  a  function  of  time.  It  was  found  to  change 
with  time  due  to  the  vertical  directionality  of  the 
projector.  The  entire  10  ms  interval  of  reverbera¬ 
tion  resulted  from  ill ’imi nation  of  the  surface  by 
the  main  lobe  from  the  vertical  directional  it v 
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pattern  of  the  projector. 

Comparison  to  Theoretical  Model 

1.  The  model  was  found  to  agree  well  with  the  spatial 
dependence  of  the  normalized  envelope  of  the  co- 
variance  from  the  horizontal  array.  This  was  especially 
significant  since  the  same  model  also  accurately 
predicted  the  horizontal  covariance  obtained  by 
Frazer,  even  though  the  two  levels  of  covariance 

were  different. 

2.  The  model  also  correctly  predicted  that  the  envelope 
of  the  normalized  covariance  would  not  change  with 
time  over  the  10  ms  interval.  The  phase  was  pre¬ 
dicted  to  be  constant  and  near  zero,  as  expected, 
since  the  model  assumes  a  uniform  scatter  density. 

The  measured  phase  was  constant  for  significant 
levels  of  covariance  but  was  non-zero;  thus  the 
actual  measured  covariance  was  less  than  predicted 
because  of  the  non-zero  measured  phase.  The  theo¬ 
retical  model  did  not  predict  the  change  in  phase 
with  time  for  the  low  levels  of  covariance  on  the 
horizontal  array. 

3.  The  envelope  of  the  covariance  as  a  function  of 
spatial  separation  for  the  vertical  array  was 
greatly  over-estimated  by  the  model.  However,  the 


model  did  correctly  predict  the  change  in  phase 
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as  a  function  of  time  for  the  vertical  array,  and 
the  change  in  phase  as  a  function  of  vertical  separ¬ 
ation. 

1*.  The  dependence  of  the  covariance  on  time  difference 
was  also  examined  by  the  model.  Very  good  agree¬ 
ment  was  obtained  between  the  model  and  the  experi¬ 
mental  results.  The  shapes  of  the  envelopes 
of  the  auto- covariance  and  horizontal  cross¬ 
covariance  were  predicted  very  accurately.  There  was 
also  good  agreement  between  the  phases  as  a  function 
of  t.  The  envelopes  of  the  cross-covariance  of  the 
vertical  array  were  overestimated,  as  was  mentioned 
earlier,  but  the  value  of  t  for  which  the  envelopes 
went  to  zero  was  correctly  predicted. 

5.  The  variance  was  also  computed  theoretically  for 
the  entire  10  ms  interva1  on  all  13  channels.  Good 
agreement  was  obtained  with  the  experimental  results. 
This  agreement  indicated  that  the  effects  of  depth, 
grazing  angle,  and  directionality  were  being  correct¬ 
ly  taken  into  account. 

6.  An  explanation  for  the  failure  of  the  model  to 
correctly  predict  the  dependence  of  the  envelope 
of  the  covariance  on  vertical  separation  was  also 
given.  It  was  shown  that  better  agreement  could 
possibly  be  obtained  by  considering  the  scatterers 


distributed  vertically  as  well  as  horizontally  in 
the  plane  of  the  surface.  Two  methods  for  including 
a  vertical  distribution  of  the  scatters were 
presented.  One  method  considered  the  scatterers 
to  be  distributed  in  a  volume  layer  upper  bounded 
by  the  plane  of  the  surface.  The  other  method 
considered  the  scatterers  to  be  distributed  only  at 
the  surface  but  allowed  the  surface  to  vary  randomly 
in  height. 

VI .  I» . 3  Computation  of  Moments  and  Tests  for  Normality 

1.  The  mean,  variance,  skew,  and  kurtosis  were  examined, 

and  all  3200  ensembles  of  each  channel  were  tested 

for  univariate  normality.  It  was  shown  that  the 

moments  contained  an  oscillatory  component.  The 

fundamental  frequency  of  the  oscillation  appeared 

to  depend  upon  the  order  of  the  moment,  that  is, 

the  r^*1  moment  appeared  to  oscillate  fundamentally 

at  a  frequency  of  r-co  .  It  was  shown  that  the  anrli- 

o 

tude  of  the  mean  was  not  statistically  large.  It 
was  also  shown  that  the  oscillation  of  the  variance 
was  not  statistically  large  compared  to  the  difference 
component  of  the  variance.  The  oscillations  :f 
the  skew  were  shown  to  be  symr..etric  about  zero. 

Since  a  normal  distribution  has  zero  skew,  it  war- 
shown  that  deviations  from  normality  iue  to  skew 
different  from  zero  were  due  to  the  or r i 1 1  at  cry 
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component  of  the  skew.  However,  this  was  found  not 
to  be  the  case  for  the  kurtosis.  A  normal  kurtosis 
has  a  value  of  3,  but  the  measured  kurtosis  was  not 
symmetric  about  a  value  of  3.  Rather,  it  was  gen¬ 
erally  larger  than  3.  Thus  deviations  from  normal¬ 
ity  were  due  to  a  larger  than  normal  kurtosis. 

2.  Several  tests  for  normality  were  performed.  The 

ensembles  were  generally  found  to  be  non-normal  due 
to  a  slight  amount  of  skew  and  a  significantly  large 
kurtosis.  This  result  is  in  contrast  to  other 
studies  which  have  indicated  that  surface  reverber¬ 
ation  was  normally  distributed. 

^  Development  of  Measurement  Techniques  and  Data  Validation 

1.  Tests  for  randomness  and  homogeneity  were  selected 
for  power  and  computational  efficiency.  Ensembles 
of  1000  samples  were  found  to  be  random  (independent) 
but  inhomogeneous.  Much  of  the  inhomogeneity  was 
removed  by  restricting  the  ensembles  to  the  last  500 
samples.  The  inhomogeneity  was  attributed  to  a 
change  in  reverberation  intensity  during  the  collec¬ 
tion  of  the  data.  Thus  it  was  possible  to  generate 
statistically  valid  sample  ensembles. 

2.  A  technique  was  demonstrated  which  reduced  the 


measurement  errors  to  an  acceptable  level. 


VI.  5 


Conclusions 


Now  that  the  work  performed  in  this  study  has  been  reviewed  and 
the  results  summarized,  a  number  of  conclusions  can  be  drawn.  The 
conclusions  will  also  be  organized  according  to  the  objectives  of  the 
study. 


VI. 5-1  Computation  of  Vertical  and  Horizontal  Covariance 

1.  Significant  levels  of  covariance  can  extend  out  tc 
relatively  large  separations.  It  is  sometimes 
assumed  that  signal  processing  gain  can  be  achieve.: 
by  lengthening  a  receiving  array  since  the  noise 

is  incoherent  over  the  extended  length  while  the 
signal  remains  coherent.  It  has  beer,  shown  here  tha 
in  some  cases  this  assumption  can  be  very  misleading 
Thus  simplistic  assumptions  in  the  design  of  sor.ar 
systems  may  lead  to  performance  predictions  which 
are  not  achievable  in  practice.  More  accurate 
designs  and  performance  predictions  can  be  obtained 
by  considering  results  such  as  presented  in  this 
study. 

2.  It  is  also  cbvious  that  the  covariance  between  *vc 
receivers  is  greatly  affected  by  the  relative  men¬ 
tation  of  the  receivers,  and  net  just  by  their 

|  separation.  In  particular,  the  covariance  bet  veer. 

;  vertically  separated  elements  is  much  larger  *  har. 

f 


between  horizontally  separated  elements 


the  case  examined.  The  time  dependence  of  the 
covariance  is  also  greatly  affected  by  the  orienta¬ 
tion  of  the  receivers.  The  horizontal  covariance 
was  constant  with  time,  while  the  vertical  co¬ 
variance  slowly  oscillated  with  time.  The  vertical 
covariance  also  depended  upon  the  location  on  the 
array.  Both  these  effects,  the  time  dependence 
and  the  location  dependence,  can  be  equated  to  a 
dependence  upon  grazing  angle.  The  normalized 
covariance  for  the  horizontal  array  was  unaffected 
by  grazing  angle  over  the  limited  range  examined, 
while  the  normalized  covariance  for  the  vertical 
array  was  greatly  affected.  Thus  geometrical 
parameters  such  as  spatial  distribution  of  the 
sensors  and  grazing  angle  can  have  a  very  signifi¬ 
cant  effect  upon  covariance. 

It  can  also  be  concluded  from  this  study  that  sonar 
parameters  can  determine  some  of  the  characteristics 
of  the  covariance.  The  temporal  extent  of  the 
covariance  is  determined  in  part  by  the  length  of 
the  transmitted  pulse  after  it  has  been  filtered  by 
the  transmitting  and  receiving  apertures.  Thus 
the  pulse  length  and  the  frequency  responses  of  the 
transmitter  and  receivers  affect  the  temporal  extent 
of  the  covariance.  It  has  been  shown  in  other 


studies^ ’ ^  that  the  frequency  spectrum  of  the 
transmitted  signal  can  also  affect  the  temporal 
extent  of  the  covariance.  It  has  been  inferred  by 
a  comparison  of  the  results  from  this  study  and  the 
study  by  Frazer  that  the  level  of  the  covariance 
between  horizontal  receivers  can  be  affected  by  the 
projector's  horizontal  beamwidth.  The  projector's 
vertical  beamwidth  certainly  affects  the  level  of 
the  unnormalized  covariance.  Therefore,  sonar 
parameters  such  as  pulse  length,  directionality, 
aperture  response,  and  signal  spectrum  affect  the 
covariance  of  surface  reverberation. 

Based  on  the  non-zero  phase  of  the  horizontal  co- 
variance,  it  was  inferred  that  surface  scatter 
ensity  was  not  uniformly  distributed.  A  non- 
uniform  scatter  density  is  most  likely  attributable 
to  a  shadowing  of  scatterers  as  a  function  of  azi¬ 
muthal  angle  to  the  surface.  Thus  the  shadowing 
function  would  depend  upon  such  factors  as  the  ancle 
between  the  wind  direction  and  the  direction  the 
array  is  pointing,  the  wave  height  distribution,  and 
the  wave  spectrum.  It  was  also  shown  that  the 
vertical  covariance  most  likely  depends  upon  wave 
height  distribution.  Thus,  although  the  present 
study  had  very  little  control  over  environmental 
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conditions ,  several  observed  phenomena  could 
possibly  be  explained  by  environmental  effects. 

Comparison  to  Theoretical  Model 

1.  A  realistic  model  must  be  able  to  correctly  take  into 
account  all  these  various  parameters  which  affect 
the  covariance  of  reverberation.  It  can  be  concluded 
that ,  at  least  over  the  range  of  parameters  consid¬ 
ered,  the  model  used  in  this  study  is  correctly 
taking  into  account  the  various  sonar  and  geometrical 
parameters ,  but  needs  to  additionally  take  into 
account  some  environmental  parameters. 

2.  In  spite  of  this  additional  need,  a  number  of  proper¬ 
ties  of  the  covariance  can  already  be  predicted. 

The  horizontal  covariance  is  being  modeled  reason¬ 
ably  well  (with  the  exception  of  non-zero  phase). 

The  temporal  extent  of  the  covariance  and  the  change 
of  the  variance  and  covariance  with  time  is  being 
predicted.  Thus  many  aspects  of  the  covariance 
can  be  predicted  without  environmental  information. 

3.  The  present  study  has  also  served  to  reinforce 
the  basic  assumptions  and  approach  of  the  point- 
scatter  model.  The  many  properties  of  the  co- 
variance  which  could  be  predicted  by  the  model  indicate 
the  reasonableness  of  its  basic  approach.  In 


addition,  a  reasonable  explanation  of  the  failure 


of  the  model  to  predict  the  vertical  covariance 
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co  Id  be  made  from  the  model,  indicating  a  failure 
not  of  the  basic  model,  but  of  its  current  imple¬ 
mentation.  Therefore  it  can  be  concluded  that  this 
point-scatter  model  is  a  reasonable  model  of 
surface  reverberation  which  has  proven  to  be  accurate 
under  certain  conditions. 

Computation  of  Moments  and  Tests  for  Normality 

1.  From  the  comparison  of  the  difference  component 
of  the  variances  and  actual  variances,  it  can  be 
concluded  that  the  sum  term  of  the  variance  is 
measurable.  In  this  case,  it  was  small  in 
comparison  to  the  difference  term  of  the  variance, 
but  it  may  not  be  small  in  every  instance.  It  is 
not  known  to  what  extent  the  sum  term  of  the  var¬ 
iance  (or  covariance)  will  affect  the  performance 
of  signal  processing  schemes  which  depend  upon 
information  about  the  covariance  of  the  noise  field. 

2.  It  was  also  observed  that  all  the  moments  examined 
contained  an  oscillatory  component  which  had  a 
varying  amplitude.  It  was  shown  that  although  the 
mean  oscillated,  the  extent  of  its  amplitude  was 
not  statistically  different  from  zero.  Thus  from 
a  statistical  viewpoint  the  mean  can  be  considered 
to  be  zero.  That  is,  the  random  amplitude  of  the 
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oscillation  was  not  significantly  different  from 
zero  in  comparison  to  the  magnitude  of  the  variance. 
In  this  case  it  can  be  concluded  that  the  moments 
examined  were  not  con: tant  with  time,  containing 
both  a  random  magnitude  which  changed  slowly  with 
time  and  a  component  which  oscillated  rapidly  with 
time.  However,  the  magnitudes  of  these  oscillations 
were  not  large. 

3.  Although  density  estimates  were  not  made  for  these 
data,  the  data  were  tested  for  normality  and  found 
to  be  non-gaussian  due  to  a  larger  than  normal 
kurtosis.  Tests  for  normality  on  other  surface 
reverberation  data  have  shown  them  to  be  gaussian. 
However,  it  is  not  well  known  under  what  conditions 
gaussian  behavior  is  to  be  expected.  Thus  it  is 
concluded  that  surface  reverberation  can  exhibit 
non-gaussian  behavior,  although  the  reasons  for 
this  behavior  are  unknown. 

VI . 5 . U  Development  of  Measurement  Techniques  and  Data  Validation 

1.  Considering  the  validation  of  the  sample  ensembles, 
one  can  conclude  that  random  ensembles  can  be 
generated,  but  that  care  must  be  taken  to  ensure 
their  homogeneity.  In  particular,  changes  in  the 
reverberation  intensity  during  the  generation  of  the 
ensemble  data  can  lead  to  inhomogeneities.  Thus 
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specific  tests  for  homogeneity  should  be  performed. 
VI .6  Accomplishments  and  Contributions 

The  summary  and  conclusions  that  have  just  been  made  for  this 
study  lead  to  the  identification  of  several  accomplishments  and  contri¬ 
butions  which  have  resulted  from  the  study: 

1.  The  covariance  of  surface  reverberation  from  a  pulsed  source 
was  measured  simultaneously  in  both  the  horizontal  and  ver¬ 
tical  orientations  and  found  to  be  quite  different.  A 
detailed  analysis  of  these  differences  was  provided. 

2.  It  was  shown  that  the  covariance  of  surface  reverberation 
depends  upon  many  environmental,  geometrical,  and  sonar 
parameters . 

3.  The  accuracy  of  the  theoretical  model  was  verified  for  the 
horizontal  covariance.  It  was  shown  that  the  model  if: 
correctly  taking  into  account  many  of  the  parameter.:  w:.i  ch 
affect  surface  reverberation. 

1*.  A  deficiency  in  the  model  was  identified.  The  rear  >n  :’->r 
the  deficiency  was  attributed  to  the  failure  of  the  •..rr-c.' 
implementation  of  the  model  to  take  into  account  cor*  air; 
environmental  conditions.  The  ability  of  *he  model  *  :■ 
explain  the  failure  of  this  implementation,  alone  wi*h  *  he 
ability  of  the  model  to  correctly  predict  many  aspect r-  'f 
the  covariance,  has  provided  additional  confidence  *!.• 
point-scatter  approach  of  the  model. 

5.  Univariate  moment s  :>f  ‘he  rev-rlnrat ;  t.  were  q.m-r.  v 
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examined.  Oscillation  of  the  moments  was  demonstrated. 

6.  The  sum  term  of  the  variance  was  shown  to  be  measurable, 
thus  bringing  into  question  the  indiscriminate  application 
of  the  difference  component  approximation  of  the  covariance. 

7.  It  was  demonstrated  that  surface  reverberation  can  be  non¬ 
gauss  i  an. 

8.  It  was  shown  that  large  ensemble  sizes  (500)  could  be 
constructed  and  validated. 

VI .7  Recommendations  for  Further  Study 

As  is  typical  of  most  studies  of  this  type,  as  many  questions 
have  been  raised  as  have  been  answered.  Several  directions  for  further 
research  have  become  obvious  during  the  course  of  the  study.  Thus 
it  seems  reasonable  to  conclude  with  a  list  of  recommendations  for 
further  study: 

1.  An  attempt  should  be  made  to  extend  the  theoretical  model 
so  that  it  will  correctly  predict  the  dependence  of  the 
covariance  on  vertical  separation.  It  is  felt  that  the 
greatest  potential  for  doing  so  lies  in  the  proper  inclusion 
of  the  statistical  distribution  of  the  surface  wave  heights. 

2.  The  model  should  also  be  extended  to  predict  the  non-zero 
phase  of  the  horizontal  covariance.  The  addition  of  a 

shadowing  function  which  depends  upon  the  direction  of  the 
waves,  the  wave  spectrum,  and  the  wave  heights  to  effectively 
produce  a  non-uniform  scatter  density  appears  to  be  a 
reasonable  approach. 


* 
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3.  In  order  to  validate  these  extensions  to  the  model,  addi¬ 
tional  experimental  data,  including  m  >urements  of  wind 
direction,  wave  heights,  and  wave  spectra,  will  be  necessary, 

4.  There  are  several  advantages  to  conducting  experiments  in  a 
controlled  environment  such  as  a  tank.  The  effects  of 
environmental  conditions  may  be  studied  individually  and 
with  better  control. 

5.  It  would  also  be  of  interest  to  extend  the  model  to  predict 
other  moments  in  addition  to  the  covariance. 

6.  Some  unexpected  effects  were  observed  when  the  covariance 
was  low.  It  was  speculated  that  secondary  scattering  may 
be  their  cause.  Although  these  effects  appear  to  have 
little  practical  significance,  further  study  of  them  may 
be  of  interest. 

7.  Additional  data  at  different  pulse  lengths,  grazing  angles, 
and  projector  horizontal  beamwidths  would  help  to  confirm 
some  of  the  tentative  conclusions  drawn  from  this  study 
concerning  the  effects  of  these  parameters. 

8.  The  model  has  yet  to  be  validated  for  other  than  pulsed  CW 
signals.  Since  experimental  results  have  shown  that  other 
transmit  types,  such  as  linear  frequency  modulated  trans¬ 
missions,  produce  results  which  differ  from  the  results  from 
pulsed  CW  transmissions,  it  would  be  useful  to  validate 

the  model  under  these  different  conditions. 

9.  Since  it  has  been  shown  that  reverberation  is  sometimes 
gaussian  and  sometimes  non-gaussian,  it  would  be  very  useful 
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to  determine  under  what  conditions  a  departure  from  normal¬ 
ity  can  be  expected. 

10.  Improved  criteria  for  accepting  or  rejecting  an  entire  set 
of  ensembles  for  a  statistical  hypothesis  test  should  be 
developed . 

11.  An  examination  of  the  power  and  cross-power  spectra  in 
conjunction  with  the  covariance  would  be  useful. 


APPENDIX  A 


DERIVATION  OF  THE  CHARACTERISTIC  FUNCTION  AND  MOMENTS 


The  following  is  a  derivation  of  the  characteristic  fun  •* i 

and  moments  of  the  reverberation  process  for  one  and  two  receivers 

The  derivation  follows  closely  Ref.  27,  but  in  abbreviated  form. 

The  receiver  output  at  time  t.  due  tc  a  single  point  sea*  *  <•: 

at  point  A  is  denoted  by  U.(t,A.).  The  receiver  output  from  all  ; 
<J  J  J 

scatterers  which  were  illuminate!  is  lust  the  sum  of  ail  eler.er.*  a; 


scattered  returns  U.(t,A 
J  J 


M(  t " ) 

x(t.t')  =  2  U  (t,A  ) 

4_1  J  J 

where  M(t')  is  the  number  of  contributing  scatterers  at  time  t'.  What 
needs  to  be  determined  here  is  the  covariance  of  X  between  two  times 
t^,  tg.  The  procedure  will  be  to  determine  the  1 oint  density  function 
of  X(t^,’  and  X(t2'  in  terns  of  the  Foisson  distribution  cf  the  sea’terers 
and  then  Fourier  transform  the  joint  density  function  to  get  the  joint 
characteristic  function,  from  which  the  covariance  can  be  calculate! 
by  taking  the  appropriate  derivatives. 

The  joint  density  function  cf  the  two  signals  >ht,)  and  X,t  '' 
and  the  number  of  contributing  scatterers  in  the  ~ra*ial  interval 
at  the  tine  t"  is 

w(x.  ,t,  ,x2,t; 1 1 

which  can  be  written  in  terms  of  the  joint  conditional  density  function 


j 
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W(X1,t1,X2,t2|H,t',A) 


as 

W(X  ,t  ,X2,t2,N|t',A)  =  P(N|t',A)W{X  ,tltX2,t2|N,t',A) 

where  P(N|t',A)  is  the  probability  function  of  the  discrete  random  var¬ 
iable  N.  The  joint  density  function  of  X^  and  X2  alone  is  just 

00 

W(X1  ,t  ,X  ,t Jt'.A)  =  2  W(x  ,t  X  ,t  ,N|t',A) 

1  1  ^  ^  N=0  1 


=  Z  P(N|t',A)W(X  ,t  X  ,t  |N;t',A)  (A.l) 

N=0 


since  N  is  discrete. 

P(N | t " ,A )  is  the  probability  of  N  events  in  the  interval  A.  The 

basic  assumption  is  that  this  probability  is  Poisson  distributed,  i.e., 

for  each  point  X^  in  A,  the  probability  that  exactly  N  events  occur  in 

a  sufficiently  small  interval  AX  about  X  is 

J  J 

P(H|t',AXj)  =  [p(t',X  )AX  ]Nexp[-pAX  ]/MJ 

where  p(t',X  )  is  the  density  of  the  process. 

By  taking  the  Fourier  transform  of  Eq.  (A.l)  for  the  smaller 
interval  AX t ,  we  get 


F2( i?l  ,ti  ,i^2  ,  t2|t',AXj  )  = 


oo 


1 

N=0 


(A. 2) 


where  F2  is  the  Joint  characteristic  function.  From  the  definition  of 
the  characteristic  function  and  the  receiver  output  X: 


F2(i?1,t1,iC2,t2|t',A)  =  n  F2  (iC1,t1,i52,t2|t',AAJ) 

=  IT  expfpAA j  ^expi(C1U1+^2U2)-l}] 

=  exp[^pAAj  (expi(C1U1+C2U2)-l)] 

J 

In  the  limit  as  AA  -*0,  this  becomes 
J 

F2(iC1>t1,i52 ,t  2|t',A)  = 

exp[/Ap(t',X)^expi(eiU(t1,A)+e2U(t2,A))-l^  dX] 

This  is  the  joint  characteristic  function  of  the  two  receiver  outputs 
X(t^)  and  X(t^).  The  characteristic  function  of  a  single  receiver 
output  would  just  be 

^(iCi^ift'jA)  =  exp[/Ap(t%  A)^'expi(C1U1)-l^dA] 

The  expected  value  of  X(t^)  is 

<x(tl»  ■  -1  A  F1 


*  -1 V  ii,  <ex[,if;i'Jr1>d,F1|t=c 

=  -i  /AP  <iU1expiC1U1>  dAF1lc=0 


=  /AP  dA 


Likewise  the  expected  value  of  the  product  of  the  two  receiver  outputs 
/(t^)  and  X( t2 )  is 


<x(t1)x(t2)>  = 


«1«2  F2Y?2-"° 


which  turns  out  to  be 


/Ap  <V:>  dX  +  /Ap<Ui>  dX  /A‘  dX 


Thus,  the  temporal  covariance 

Kx(t'-rtp|t'}  =  (x(t1)x(t2))  -  (x(t,  ))(r.Up'} 

is  just 

/Ap(t'.X)  <'u(t1,A)  U  (t2,X))  d.x 


This  is  the  t-enporal  covariance  of  a  single  receiver.  The 

the  joint  spatial-temporal  covariance  of  the  outputs  X^\' 

(2) 

X  “  (t£)  of  two  receivers  follows  exactly  the  sane  argumen 
yield  the  similar  result: 


where  X^^(t^)  is  the  output  of  receiver  1  a*  time  t,  and 
the  output  of  receiver  2  at  time  t0.  A.  ,  is  the  region 
which  contributed  to  the  si  mu.  a*  both  receivers  and  .  .  . 
density  of  these  scatterers. 


f 

I 

! 

1 

t 

L 


APPENDIX  B 


THE  ELEMENTARY  SCATTERED  WAVEFORM 


Utilizing  the  assumptions  stated  in  Chapter  III,  the  form  for 
the  elementary  scattered  wave  U(t,A)  is  developed.  The  procedure 
amounts  to  convoluting  the  transmitting,  scattering,  and  receiving 
apertures  with  the  appropriate  signal  and  solving  the  homogeneous  wave 
equation  in  the  medium. 

Beginning  with  the  transmitter,  the  pressure  density  at  the 
point  4  on  the  surface  of  the  transmitter  (Fig.  B-l)  is  just  the  con¬ 
volution  of  the  input  electrical  signal  with  the  response  function  of 
the  transmitter.  In  the  frequency  domain  this  is 

0T(t,O  =  aTU)f  TT(f)S.n(f)elaltdf 
—  00 

where 


S.  (f)  is  the  frequency  spectrum  of  the  input  electrical 


in 


signal , 


T^,(f)  is  the  frequency  response  function  of  the  transmitter,  and 
a^,(C)  is  the  transmitter  aperture,  and  is  the  spatial  Fourier 

transform  of  the  beam  pattern.  It  is  assumed  that  the  aper¬ 


ture  is  frequency  independent  and  that  is  the  same  for 

all  points  4  on  the  transmitter. 

The  incident  pressure  at  any  point  R^,  from  the  transmitter  in 

the  homogeneous  medium  is  given  by  the  solution  of  the  wave  equation 

V2  P.  -  P.  /c2  =  -G  ( t , O 
inc  me  T 

2h  U 


The  solution  is 


P. 

inc 


u.V  -  L  o, 


(t  _L  f)  H- 
lt  c’5J  W 


where  is  given  above  and 


Applying  the  farfield  assumptions  that 

1  =  1-  =  r  _  ^T 

r  R^,  an  c  c  c 


and  substituting  in  from  above  gives 


pi„=(t’V  -  ^  I-.  Vf  ST)TT(«sln(f)«1"'1  =  'if 


iu(t-^) 


where  the  integral  over  £  has  been  incorporated  into  the  beam  pattern 
A,j,  as  the  inverse  Fourier  transform  of  a  . 

To  get  the  scattered  pressure,  it  is  first  necessary  to  consider 
the  reradiated  pressure  density  over  the  scatterer.  This  pressure 
density  is  just  the  convolution  of  the  incident  pressure  on  the  scatterer 
with  the  response  function  of  the  scatterer,  Just  as  was  done  with  the 
electrical  signal  and  the  response  function  of  the  transmitter.  It  is, 
this  time,  in  the  time  domain: 

Gscat  =  fL  h(T’t|Hr)pinc(t-T,RT)dT 

where  5 ( R )  is  the  delta  function  and  is  the  spatial  aperture  of  the 
point  scatterer,  and  h(T,t|R^)  is  the  response  function  of  the  scatterer. 
The  response  function  is  time  varying  and  thus  spreads  the  incident 
wave  in  frequency  as  well  as  time  delay. 


The  scattered  pressure  is  again  the  solution  of  the  homogeneous 


scalar  wave  equation 


VP  -  P  ,  ,  2  =  -G  ( t  ,R I R  ) 

scat  scat/c  scat  1  T 


The  solution  is 


P  ,(t,r')  =  /  G  (t-— ,R|r_)  p 

scat  JV  scat  c  T  4u 


where  r  =  r'  -  R  and  r'  is  a  vector  from  the  scatterer  to  receiver. 
Substituting  in  for  G  and  integrating  over  the  delta  function 

S  C  aLTy 

6(R)  gives 


Pscat(t,r  ^  =  ,,  ,2  .  L  lRT)AT(cVTT(f! 

(Utt)  RTr 

r " 

x  S.n(f)e1  ^t__c  r}df 


where  the  impulse  response  function  h  has  been  absorbed  into  y  by 
integrating  over  t;  that  is, 


y(f,t-^|RT)  =  C  h(x,t-^!RT)e-lwTdT 

All  that  is  left  to  do  now  is  to  convolve  this  scattered  pressure 
with  the  receiving  aperture  to  get  the  elementary  scattered  return, 
namely : 


U(t,BR)  =  /  J?.R(n)/:.TR(f-)Fs.at(:-,Rs;g"  If 

R 

where  a  Cn)  is  the  receiver  aperture,  T  ( f ' )  is  the  receiver  frequency 
R  R 

response,  and  Pscat^P»Prj^  i-s  't^le  Fourier  transform  of  ^(t,R^).  It 
is  assumed  that  the  receiver  is  in  the  farfield  of  the  scatterer  so 


2U8 


that  the  farfield  assumptions  made  previously  are  also  applied  here. 
Substituting  the  equation  for  P  (f',R^  into  the  equation  for  the 

S  C  n 

elementary  scattered  return  gives 


"“•V  ■  — 5 —  C. df'/l  dfTR<f')V'TSBWf'f'-f|V 

(i*n) 

R 


«  VjVTT(f)Si„(f)e 


iu'(t— 5  )-i«— 
c  c 


where  YCfjf'-flH^)  is  the  Fourier  transform  of  y  with  respect  to  f'-f. 

It  is  a  measure  of  the  amount  a  frequency  component  f  is  shifted  by 
f'-f. 

The  simplest  approximation  to  make  is  that  the  scatterer  does 
not  spread  the  frequencies  of  the  incident  wave.  In  this  case 

Y(f,f'-f|RT)  =  Y(f|^T)5(f'-f) 
and  the  elementary  scattered  return  reduces  to 

"‘‘•v  ■  u  X ,  r-  dfyf)vfvi(fiyv^r> 

X  s . 
m 

To  a  good  approximation,  the  beam  patterns  are  independent  of  frequency 
for  a  narrowband  signal.  If  it  is  also  assumed  that  the  point  scatterer 
attenuates  all  frequencies  of  a  narrowband  signal  equally,  then 
is  independent  of  frequency,  Y(f)  =  Y.  Thus,  the  point  scatterer  is 
rso  i a  perfect  point  reflector. 

With  these  twc  assumptions  t.he  elementary  scattered  return  is 


.  ,,  V*r, 


HR+h,r  ^ 


YA  (RjA  (R  )  ioj(  t-— — 

U(t,RR)=  --R  g  —  -  /  ^  dfT  (f)T  (f)S  (f)e  c 

(M  Rp^ 


This  integral  is  just  a  function  of  the  receiver  and  transmitter  fre¬ 
quency  response  and  the  spectrum  of  the  input  electronic  signal.  It 
represents  the  input  signal  after  it  has  been  filtered  by  the  trans¬ 
mitting  and  receiving  array,  and  is  measured  experimentally  for  the 
particular  transmitter  and  receiver  used.  As  such,  it  is  a  bandiimited 
signal  and  can  be  represented  in  the  form 

V(t  )  =  a(t  )cos[w  t  ( t  )] 
r  r  o  r  r 


where 


t  =  t  - 
r 


PR  +  \ 


is  the  retarded  time.  Thus,  the  elementary  scattered  return  is 


U(t,Rp)  = 


(U)2RRR,r 


APPENDIX  C 


PET  AT  LC  OF  THE  COMPUTATION  OF  THE  DIRECTTV  l':'v 


FUNCTIONS  AND  TRANSFORMATION  FACTOR 


The  directivity  functions  for  the  projector  unci  receiver  us 
riven  in  Chapter  ITT  need  to  be  expressed  in  terms  07’  1  he  variables  of 
integration  t.  and  .  This  requires  a  change  of  variables  from  the 
coordinate  system  centered  on  an  element  to  the  variables  of  integration 
t  and  <f>.  The  following  is  a  description  of  the  transformat  ion  rela- 
t.  i  onships . 

The  products  sincyosP^,  sinu^si  nP.^,  r.ina.^  ,  and 

aina  sinB„.  can  be  related  to  the  vectors  R  and  R  by  the  following: 
nl  nl  J  n i 


sinaTcos3T  =  ~ 


sina^sinp,^, 


sinoi .  cosB,, . 
Ri  Ri 


r.lna_ .  sinp_, . 
Ri  Ri 


VO  O  *"> 

RT,+RTy+RTs. 


R  .  =  V  R‘  +R‘  .  +H‘  . 
Ri  Rix  Riy  Ri: 


R_  ,  R_  .  and  R_  are  the  Cartesian  coordinates  of  the  vector  R  . 
Tx  Tv  T"  I 


and  R  ,  R_.  ,  R_.  are  the  Cartesian  coordinates  of  the  vector  R„. 

Rix  Riy  Riz  Ri, 

as  illustrated  in  Fig.  III-3.  Since  the  coordinates  A  of  the  point 
scatterer  are  given  in  terms  of  the  variables  t  and  <t>,  it  is  necessary 
to  relate  "R^,  and  R  to  these  variables  in  order  to  compute  these 
products . 


Consider  a  coordinate  system  on  the  surface  directly  above  the 
center  of  the  array  of  receivers  and  projector,  as  shown  in  Fig.  C-l. 

The  vector  A  is  the  vector  from  the  origin  to  the  point  scatterer, 
where  the  components  of  A,  the  unit  vector  in  the  direction  of  A,  are 
( sin<J> ,  cosb,  0).  The  vector  R^.  is  the  vector  from  the  ith  receiver  to 

—Ik 

the  point  scatterer,  and  R^  is  the  vector  from  the  projector  to  the 
point  scatterer.  The  magnitude  of  R^  is 

|RtI  =  c*tT 

The  coordinate  system  at  the  surface  can  be  related  to  the  coordinate 
system  centered  at  the  projector  or  receivers  by  a  translation  to  the 
center  of  the  array  representing  the  depth  of  the  center  of  the  array 
below  the  surface,  a  rotation  of  the  array  through  an  angle  to  produce 
the  desired  grazing  angle  with  the  surface,  and  a  translation  along  the 
surface  of  the  array  to  the  location  of  the  projector  or  receivers  withir 
the  array.  Thus  the  vectors  R^,  and  R  .  are  relates  to  the  vector  A  tv 

i  r*l 

the  following  transformations: 

=  3(  A-u)  -  vT 

4  “  S(t-r)  - 


ORIGIN  OF 
COORDINATES  AT 
WATER'S  SURFACE 


z 


FIGURE  C  l 
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where 

"u  is  the  vector  from  the  origin  on  the  surface  to  the  center  of 
the  array  at  a  depth  h  below  the  surface: 


is  the  vector  from  the  center  of  the  array  to  the  projector, 
and  depends  on  the  location  of  the  projector  on  the  array; 

*Ri  iS  the  vector  from  the  center  of  the  array  to  the  ith 

receiver,  and  depends  on  the  location  of  the  receiver  on  the 
array ; 

and 


S  is  the  rotation  matrix  representing  a  rotation  through  the 
angle  ip : 


/  1  0  0 

S  =  I  0  cosip  simp 

\  0  -sini p  cosi^ 

Since  the  unit  vector  X  is  known,  it  is  only  necessary  to  determine  |x 

in  order  to  determine  the  vector  X  and  thus  compute  R,^  and  R  .  To 

determine  |x|,  consider  first  the  vector  R^  for  X  =  0.  In  this  case 

Rrp  points  to  the  origin  of  the  surface  coordinates  and  is  given  by 

Rt  =  -(Su+vT),  for  X  =  0 

The  negative  of  this  vector  is  represented  in  Fig.  C-l  by  the  vector 

"rT,  i.e.  , 


rT  =  Su  +  vT 

Next,  consider  the  vector  R^,  for  X 


X.  This  vector  is  designated  by 


\  ;  and  is  riven  by 

>  '  =  f.f.A-nl  -  v  »  ;  '  -  :?u  -  v,,,  -  -  Tv, 

From  Fig.  C-l  it  ran  be-  seen  t  bn'  the  vert  1  in  '  u.ate 

system  centered  at  ;.)ie  er,  mat  ed.  by  '  \  [.■  jus? 

y\  **v 

'  '  =  +  r?  -  SA 

Tf  a  line  is  drawn  t’ror.  the  coordinate  origin  at  the  riwjeete r  to 
sect  the  vector  A  at  a  rirht  angle,  then  the  distance  between  t i:o 
section  and  the  coordinate  origin  at  the  surface  is  r;v-  .  by  the  i 

-•'V 

product,  of  r  and  •  The  distance  between  the  intersect  i  on  and  > 
coordinate  origin  at  the  projector  is  then  given  by 


£TZ  y.«p 

r  •  r  —  l  r  *  \  ) 

T  T  ri 


Thus  the  distance  from  the  intersection  t .»  t  !.•-  p.  iut  scat  ;  c: .  t. 


|PTi2-(?v-t,.-(rT-\')") 


and  the  magnitude  of  A  is.  t.nerefore  v;v-  n  by 


A  = 


.♦A  +  v 


-n 


’V 


"r  •  A  ’  V  ) 
T 


Since  the  unit  vector  A  is  known,  t.hi.-h  <  lie  vector  A  has  thus  been 

- -4  _ k 

mined.  Knowing  A,  the  vectors  R  and  K^.  can  be  calculated  from 

— *  — k 

given  coordinate  transformations.  The  components  of  K,;,  anu  R.  .  nr 
used  to  compute  the  proda-*'  s,  r-'tuiri'  I  by  t  he  din-civily  rih'Ul-it 
Thin  thus  allows  comrut  at.  i  on  of  1  he  )  i  rcct.  i  »  i  <  ,v  fut.  *t.  ions  to"  t  :,>■ 


nrojecl.or  and  receivers  iti  t.h>*  correct  coordinate  system. 

facto-  w  can  also 


UIi* 1 

he 


iet 

h>n 

e  t  heii 
i  >  ns  . 


The  trails  format,  i  on 


now  bo  express'*- 


in  '  orr.r 


of  t  and  $ .  The  surface  integration  can  be  expressed  in  terms  of  the 
surface  coordinates  X  =  | X |  and  $  by 
dS  =  XdXd$ 

It  is  now  necessary  to  transform  the  variable  X  to  the  variable  Rr^=c'bp . 
From  the  expression  above  for  X,  it  can  be  seen  that 


Thus 


dX 

dR„ 


R„ 


V 


( rT  •rT-  ( rT  •V ) 2 ) 


XdX  =  I  2  <S,2, 

W RT-(rT*rT-(rT*X  )  ) 


1  + 


. 

rT*A 


V 


=  V  RTdRT 


=  w  c  tdtm 
T  T 


Therefore 


dS  =  wc  t_dt_d$ 
T  T 


where 


w  =  1  + 


Vv 


Vp-r^-fryl  -)2) 


For  the  geometry  of  the  present  study,  v  was  very  close 


ver  *  no 


range  of  t, 


and  $  considered. 
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